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A. 1 Introduction 

The  primary  function  of  the  Thermal  Distribution  System 
Simulation  Code  (TDIST)  is  to  aid  in  the  design  and  perfor- 
mance analysis  of  a Total  Energy  System  (TES)  supplying 
thermal  and  electrical  energy  to  a metropolitan  area  con- 
taining a variety  of  energy  consumption  categories.  To 
achieve  this  end,  the  code  is  concerned  with  the  steady- 
state  and  time-dependent  energy  flow  modelling  of  up  to 
twenty  different  energy  consumer  types  and  of  the  thermal 
utility  system  (TUS)  serving  them  with  high  temperature 
water  (HTW)  from  a central  total  energy  power  plant  site. 
Since  the  major  costs  of  such  an  integrated  energy  supply 
system  are  particularly  sensitive  to  the  TUS  design,  emphasis 
is  placed  upon  simulating  the  time-dependent  TUS  performance. 
No  provisions  are  made  for  specifically  modelling  the  asso- 
ciated electrical  energy  distribution  network,  although  the 
community's  electrical  energy  demands  are  calculated  and  are 
included  in  the  overall  TES  energy  supply  analysis. 

Briefly,  TDIST  consists  of  two  distinct,  but  intercon- 
nected, modelling  packages.  The  first  is  a set  of  building 
energy  flow  calculation  subroutines  which  determine  the  space 
conditioning,  domestic  hot  water,  and  non-space-conditioning 
electrical  loads  for  each  of  the  specified  energy  consumer 
types,  based  upon  steady-state  demand  calculations  scheduled 
over  time  according  to  building  usage  characteristics  and 
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ambient  weather  conditions.  These  demands  are  applied  to  a 
range  of  end-use  equipment  types  available  for  each  building 
in  the  community,  the  aggregate  demands  of  which  comprise 
the  time-varying  loads  to  be  supplied  by  the  thermal  and 
electrical  energy  distribution  networks.  The  second  model- 
ling package  utilizes  the  computed  thermal  energy  demands 
and  the  thermal  energy  output  of  the  central  TES  power 
plant  as  driving  functions  for  the  TUS.  Included  in  the 
dynamic  models  of  this  system’s  components  are  descriptions 
of  the  major  piping  loops,  the  HTW  circulation  pumps,  heat 
exchangers,  thermostatically-controlled  fluid  flow  regula- 
tion circuitry,  and  a central  stratified  thermal  energy 
storage  reservoir  designed  to  smooth  the  system’s  thermal 
energy  supply  and  demand  imbalances  over  a 24-hour  duty  cycle. 

Various  conceptually  independent  segments  of  these 
models  are  developed  in  separate  program  units  to  allow  flex- 
ibility in  their  formulation  and  integration  into  the  two 
packages.  These  sub-program  units  are  discussed  in  detail  at 
a later  point  when  the  derivations  of  the  individual  system 
cqmponent  models  are  presented.  However,  a brief  overview 
of  the  general  characteristics  of  the  overall  TES  model  is 
useful  in  understanding  how  each  of  the  sub-units  is  related 
to  the  others  and  how  they  are  integrated  to  yield  the  desired 
TES  behavioral  information. 
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A. 2 Overall  TES  Simulation  Model  Description 

For  the  purposes  of  energy  demand  modelling,  it  Is 
assumed  that  all  the  buildings  in  the  community  are  aggre-' 
gated  into  a set  of  up  to  20  representative  energy  consump- 
tion categories  based  upon  similarities  in  construction 
characteristics  and  usage  patterns.  Building  units  typify- 
ing each  of  these  consumption  categories  are  further  aggre- 
gated into  from  one  to  twenty-five  Consumer  Load  Centers, 
the  distribution  of  which  reflects  the  type  and  intensity  of 
use  to  which  each  geographical  region  within  the  community 
is  put.  Figure  A.l  illustrates  a unidirectional  grouping  of 
these  load  centers  postulated  for  a typical  large  U.S.  mili- 
tary installation,  which  is  fairly  representative  of  small 
industrial  or  suburban  civilian  communities  [1]. 

The  assumed  configuration  for  the  TUS  piping  connecting 
the  load  centers  to  the  central  power  station  is  shown  in 
Fig.  A. 2.  Each  Load  Center  Heat  Exchanger  is  a single  heat 
exchanger  which  supplies  a load  center  with  water  at  a 
temperature  consistent  with  the  requirements  of  the  end-use 
equipment  in  the  buildings  which  it  serves.  All  thermal 
energy  distribution  in  the  piping  loops  shown  in  the  figure 


is  through  the  HTW  medium,  although  the  use  of  low  pressure 
steam  generators  for  certain  process  heating  applications  is 
not  precluded  in  sections  of  the  utility  system  not  included 


in  the  TDIST  models.  The  Load  Center  Heat  Exchangers  may 
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CONSUMER 


Thermal  Utility  System  Piping  Schematic 
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represent  central  distribution  units  serving  several  build- 
ings through  (unmodelled*)  tertiary  piping  loops,  or  they 
may  represent  units  serving  individual  buildings,  but  each 
heat  exchanger  specified  in  the  system  must  have  at  least 
one  energy  consumer  associated  with  it.  Up  to  five  Load 
Center  Heat  Exchangers  are  supplied  from  a single  Secondary 
Distribution  Loop.  Because  many  of  the  consumers  may  require 
water  of  the  same  temperature  as  others,  the  Load  Center 
Heat  Exchangers  are  connected  in  parallel  as  shown  in 
Fig.  A.  2,  with  their  feedwater  and  return  flov^s  joining  at 
common  supply  and  return  headers.  A set  of  flow  control 
valves  govern  the  allocation  of  the  loop  fluid  to  each  Load 
Center  Heat  Exchanger  according  to  its  relative  energy  demand 
rate.  Up  to  five  Secondary  Distribution  Loops,  in  turn,  are 
supplied  with  thermal  energy  through  the  illustrated  Secon- 
dary Heat  Exchangers  in  series  with  the  Primary  HTW  Loop 
fluid  flowstream.  In  order  to  provide  individual  control 
over  the  fluid  flowrate  through  each  of  these  heat  exchangers 
while  maintaining  the  desired  fluid  temperatures  around  the 


♦Tertiary  piping  is  excluded  from  the  TDIST  models  since  it 
would  require  a prohibitive  number  of  additional  system  nodes 
and,  because  of  its  small  diameter  and  relatively  short  pipe 
runs,  it  is  assumed  to  have  negligible  effects  upon  the  over- 
all thermal  inertia  or  operational  stability  of  the  TUS . 
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primary  loop,  bypass  loops  are  provided  around  each  Secon- 
dary Heat  Exchanger,  with  the  relative  bypass  and  heat 
exchanger  fluid  flowrates  governed  by  control  valves  set 
according  to  the  heat  exchanger  thermal  demands.  Control 

over  the  total  primary  and  secondary  loop  energy  supply 

* 

rates  is  accomplished  through  variations  in  the  various  cir- 
culation pump  pressure  settings,  adjusting  the  loop  fluid 
flowrates  while  maintaining  relatively  fixed  loop  fluid 
temperature  distributions. 

As  is  shown  in  Fig.  A. 2,  the  thermal  energy  Storage 
Reservoir  serves  as  an  interface  between  the  Primary  HTW 
Loop  and  the  system’s  thermal  energy  supply  at  the  power 
plant's  TUS  Supply  Heat  Exchanger.  Figure  A. 3 shows  a more 

detailed  view  of  this  reservoir,  while  the  schematic  Fig.  A. 4 

. 

illustrates  the  method  by  which  the  two  different  temnerature 
sections  are  separated  for  modelling  purposes.  During 
periods  of  low  thermal  demand  and  excess  thermal  supply,  the 
reservoir  is  charged,  with  the  temperature  interface  moving 
into  the  cold  water  section,  and  the  hot  water  volume  increas- 
es. When  necessary,  hot  water  is  withdrawn  from  the  tank  to 
supplement  the  power  plant's  output.  Because  modern  technolo- 
gy is  limited  to  the  production  of  pressurized  steel  cylinders 
of  the  type  shown  with  dimensions  of  approximately  20  ft.  in 
diameter  and  70  ft.  long  [2],  the  actual  reservoir  would  con- 
sist of  a set  of  these  smaller  tanks  connected  in  parallel  to 
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Determined  independently  of  Primary  Loop  flow  to  allow  the 
storage  volume  to  vary 

Figure  A. 4 

Thermal  Energy  Storage  Reservoir  Model  Schematic 


produce  a large  storage  field.  European  experience  with 
this  method  of  energy  storage  indicates  that  minimal  mixing 
between  the  two  fluid  temperature  regions  takes  place  if 
charge/discharge  cycling  times  are  on  the  order  of  4-12 
hours  [33- 

The  assumed  mode  of  power  plant  operation  is  one  of 
constant  total  thermal  power  output,  with  the  station's 
electrical  energy  generation  following  the  community's  elec- 
trical demand  schedule.  Thus,  if  the  plant  is  rated  at 
100  MWt,  and  the  community's  electrical  demand  varies  from 
10  MWe  to  40  MWe , the  thermal  power  available  to  the  TUS 
will  vary  from  90  MWt  to  60  MWt.  This  thermal  energy  is 
supplied  directly  to  the  HTW  section  of  the  storage  reser- 
voir, from  which  it  is  removed  through  the  primary  loop  as 
called  for  by  the  TUS  demand  schedule. 

A. 3 Model  Subsystem  Descriptions 

As  is  stated  above,  the  overall  energy  system  model  can 
be  broken  down  into  several  conceptually  independent  sub- 
systems whose  behavior  can  be  studied  separately  and  then 
reintegrated  into  the  full  model.  Table  A.l  lists  the  sub- 
sections of  the  program  which  are  concerned  with  each  of  the 
separate  modelling  tasks  and  describes  briefly  their  func- 
tions within  the  total  system  model;  more  detailed  descrip- 
tions and  derivations  of  the  specific  models  used  are  presen- 
ted in  the  following  sections. 


TABLE  A. 1 


TDIST  SUBROUTINES 


Subroutine 


1.  Main  (MAIN) 


2.  Input  Scaling 
(CSC ALE) 


Function 

Data  input,  program  flow  control,  time 
monitor,  thermostat  control 

* 

Scales  input  fluid  mass  flowrates  accor- 
ding to  case  run  multiplier  (not  fully 
developed) 


Consumer  Heat 
Flow  (QLOAD) 

Building  Load 
Distribution 
(DLOAD) 


5.  Steady-State 

System  (STDY ) 


6.  Initialization 
(SET) 


Calculates  space  heating  and  cooling 
loads  for  all  building  types 

Calculates  domestic  hot  water  loads  and 
distributes  building  total  energy  demands 
among  specified  mixture  of  thermal  and 
electrical  end-use  equipment 

Establishes  steady-state  flow  rates  and 
temperature  distribution  throughout 
system 

Initializes  dynamic  analysis  sections 
to  converged  steady-state  system 
condit4  ons 


Fluid  Flowrate 
(FLOWS) 


Calculates  instantaneous  fluid  flowrates 
throughout  system,  monitors  storage 
reservoir  heat  flow 


8.  System  Pressure  Calculates  fluid  frictional  pressure 

Drops  (PRESS)  losses,  sets  circulation  pump  pressures 

9.  Dynamic  Energy  Sets  instantaneous  system  temperature 

Balance  (ENRG)  distribution  based  upon  heat  flow 
analyses 

10.  Precooler  Monitor  Monitors  heat  flow  through  power  plant 

(PREC)  TU3  supply  heat  exchanger 


11.  Output  (PRV) 


Prints  system  parameters 
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A. 3-1  Energy  System  Simulation  Control  (Program  MAIN) 


To  simplify  the  calculation  procedure  in  each  program 
unit,  the  entire  model  system  as  shown  in  Fig.  A. 3 is  broken 
down  into  a set  of  nodes  whose  boundaries  are  determined  by 
such  elements  as  pipes  and  heat  exchangers  having  different 
physical  properties,  fluid  flow  and  heat  transfer  character- 
istics and  transient  thermal  behavior.  A typical  secondary 
loop  serving  three  load  center  heat  exchangers  thus  consists 
of  six  nodes:  the  three  load  center  heat  exchangers,  two 

lengths  of  distribution  piping,  and  one  secondary  heat 
exchanger.  Appropriate  boundary  conditions  are  invoked  in 
each  subroutine  to  ensure  the  physical  continuity  of  the 
piping  loop  being  analyzed,  but  all  calculations  are  performed 
on  a purely  nodal  basis.  All  nodes  are  indexed  for  simple 
identification  of  their  parameters  and  locations  within  the 
system,  and  since  heat  transfer  between  the  two  nodes  on 
opposite  sides  of  a heat  exchanger  is  vital  to  the  dynamic 
simulation  development,  each  heat  exchanger  is  also  indexed 
separately  by  both  its  hot-and  cold-side  nodes.  The  program 
input,. for  each  node  consists  of  its  index  numbers,  codes  for 
its  type  (pipe  or  heat  exchanger,  series  or  parallel  branch, 
whether  or  not  it  contains  a circulation  pump),  its  overall 
heat  transfer  coefficient,  and  estimates  of  the  initial 


fluid  mass  flowrate,  inlet  and  outlet  fluid  temperatures, 
and  fluid  density,  specific  heat  and  viscosity.  Additional 
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distribution  system  input  data  includes  the  initial  volumes 
of  the  hot  and  cold  regions  of  the  storage  reservoir,  the 
HTW  storage  temperature,  and  specifications  of  the  power 
plant  coolant  mass  flowrate,  specific  heat,  and  inlet  and 
outlet  temperatures  at  the  system  supply  heat  exchanger. 

The  ambient  earth  temperature  is  specified  for  the  determina- 
tion of  heat  losses  from  the  distribution  piping  and  heat 
conduction  to  the  earth  from  the  various  buildings.  Consumer 
sector  input  data  includes  the  exposed  areas  and  thermal 
resistances  of  all  exterior  features  (walls,  windows,  roof, 
basement)  of  each  specified  building  type  within  the  community, 
specifications  of  infiltration  crack  dimensions  and  air  flow 
coefficients,  the  building  heights,  the  wall  and  roof  surface 
material  types,  peak  ventilation  air  flowrates,  peak  internal 
non-space-conditioning  equipment  heat  generation  rate'-,  peak 
domestic  hot  water  usage  rates,  shading  coefficients  and 
solar  absorbtivities  of  the  walls  and  roofs,  and  the  building 
orientations  relative  to  North.  The  number  of  units  of  each 
building  type  served  by  each  load  center  heat  exchanger  and 
the  total  number  of  buildings  of  each  type  not  connected  to 
the  TUS  are  input  in  conjunction  with  specifications  of  the 
types  and  distribution  of  space  conditioning  equipment  to  be 
used  in  the  system.  The  desired  building  room  temperatures, 
building  usage  specifications  and  domestic  hot  water  demands 
are  provided  in  the  form  of  hourly  schedules  and  are  used  in 
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conjunction  with  hourly  air  temperatures,  wind  velocities 
and  wind  directions  for  the  specified  time  period  - in  addi- 
tion to  seasonal  solar  position  and  radiation  intensity  data 
- to  determine  the  hourly  energy  demands  for  the  buildings, 
which  are  applied  to  the  thermal  and  electrical  energy  dis- 
tribution networks  after  being  scaled  according  to  the  speci- 
fied end-use  equipment  coefficients  of  performance.  Finally, 
a fixed  power  plant  thermal  output  power  level  and  electrical 
generation  efficiency  are  input  in  conjunction  with  hourly 
non-space-conditioning  electrical  energy  demand  data  (e.g., 
loading  from  lights,  appliances,  motors,  etc.)  to  determine 
the  net  energy  available  from  the  power  plant  to  supply  the 
building-specific  space  conditioning  and  domestic  hot  water 
energy  demands.  Control  codes  indicating  the  type  of  calcu- 
lation to  be  made  (building  load  analysis,  TUS  design  initial- 
ization, steady-state  TUS  analysis,  full  dynamic  system  ana- 
lysis), the  simulation  time  step  interval,  and  the  desired 
output  frequency  determine  the  overall  program  flow  and 
information  transfer. 

. In  the  building  load  analysis  mode  of  operation,  the 
code  computes  and  prints  out  the  hourly  space  conditioning 
loads  and  tho-ir  components  for  each  of  the  specified  build- 
ing types  to  allow  verification  of  the  building  input  para- 
meters prior  to  the  use  of  these  demands  in  the  dynamic  sys- 
tem simulations.  The  first  hourly  set  of  demands  are 
scaled  by  the  specified  end-use  equipment  coefficients  of 
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performance  and  the  resulting  loads  are  aggregated  according 
to  the  given  distribution  of  buildings  and  are  displayed  at 
the  load  center  heat  exchangers  to  allow  the  calculation  of 
TIJS  system  initialization  or  design  parameters  to  be  used  as 
input  data  for  the  dynamic  simulations. 

' * 

If  only  the  steady-state  TUS  analysis  is  required  in  a 
given  run  (e.g.,  to  check  a design  for  adhering  to  desired 
limits  on  fluid  flowrates  or  temperatures),  the  model  cal^ 
culates  only  a single  set  of  consumer  loads,  applies  these 
loads  to  specified  supply  and  return  temperature  conditions 
at  the  load  center  heat  exchangers,  and  determines  the  resul- 
ting steady-state  fluid  flows  and  temperature  distributions 
throughout  the  HTW  transmission  system  based  upon  specified 
system  thermal  and  electrical  energy  supply  conditions. 

(See  the  section  describing  Subroutine  STDY  for  a more 
detailed  analysis  of  the  actual  calculation  procedures  used.) 

If,  as  will  generally  be  the  case,  it  is  desired  to 
study  the  dynamics  of  the  TUS  behavior  over  a period  of  time 
during  which  both  the  energy  supply  and  demand  conditions 
will  vary.  Program  MAIN  becomes  the  overall  control  and 
time-monitor  for  the  remainder  of  the  subprograms.  Unless 
a detailed  system  description  including  all  the  nodal  flow- 
rates  and  temperatures  for  the  specified  initial  conditions 
is  known,  the  model  calculates  the  steady-state  system  dis- 
tributions for  these  conditions  and  uses  them  as  the  input 
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parameters  to  the  dynamic  calculation  subroutines.  A 

« 

is  set  and  Program  MAIN  sweeps  the  system  at  a rate  fixed 
by  the  specified  transient  time  step  interval.  At  each 
time  increment,  the  previously  computed  building  load  dis- 
tributions are  used  to  yield  instantaneous  thermal  and 
electrical  energy  demands  through  linear  interpolation 
between  the  hourly  data  points.  These  loads  are  applied  to 
the  load  center  heat  exchangers,  and  nodal  energy  and  momen- 
tum transfer  calculations  are  performed  to  determine  the  net 
imbalances  during  the  time  step  between  the  system  thermal 
demands  and  the  power  plant  thermal  energy  output  (given  by 
the  total  station  power  rating,  less  the  electrical . energy 
produced  to  meet  the  electrical  demands). 

In  addition  to  performing  the  program  flow  control  and 
monitoring  functions.  Program  MAIN  also  contains  the  thermo- 
stat models  which  govern  the  TUS  response  to  thermal  energy 
supply  and  demand  variations.  Chosen  as  temperatures  to  be 
monitored  by  these  thermostats  are  the  consumer  fluid  supply 
and  return  temperatures  at  the  load  center  heat  exchangers 
(to  insure  uniform  operational  parameters  for  the  consumers' 
end-use  equipment)  and  the  outlet  temperatures  on  the  primary 
loop  side  of  the  secondary  heat  exchangers  (to  insure  the 
maintenance  of  sufficiently  high  fluid  temperatures  at  the 
supply  to  the  last  secondary  loop  in  the  system).  At  the 
end  of  each  time  step,  after  the  system  temperatures  have 
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been  adjusted  according  to  the  energy  transfers  occurring 
during  the  increment,  each  of  these  reference  points  is 
checked  to  determine  if  the  temperatures  have  exceeded  their 
limits.  (Consumer  water  supply  temperatures  are  allowed  to 
vary  + 5 °F  from  their  initial  design  values;  return  temp- 
eratures are  allowed  to  vary  from  the  minimum  room  tempera- 
ture in  any  of  the  buildings  served  to  5 °F  above  their 
initial  design  values.  Primary  loop  heat  exchanger  outlet 
temperatures  are  allowed  to  vary  + 5 °F  from  their  design 
settings.)  If  any  temperature  is  found  to  have  surpassed 
its  tolerances,  steps  are  taken  to  re-adjust  the  loop  fluid 
flowrates  to  correct  the  causal  energy  flow  imbalances.  To 
accomplish  this  re-adjustment,  a set  of  dummy  loop  fluid 
flowrates  are  calculated  which  satisfy  the  instantaneous 
energy  demand  conditions  within  the  desired  loop  fluid 
temperature  distributions.  Valves  in  the  secondary  distri- 
bution loops  are  reset  to  allocate  the  total  loop  dummy 
fluid  mass  flowrates  among  the  load  center  heat  exchangers 
according  to  their  relative  energy  demands.  Independent 
dummy  flowrates  are  calculated  for  each  of  the  primary-side 
secondary  heat  exchanger  nodes  based  upon  their  thermal 
demands  and  desired  inlet  and  outlet  temperatures.  The 
largest  of  these  flowrates  is  assumed  to  occur  under  condi- 
tions of  no  bypass  flow  (i.e.,  the  total  primary  loop  fluid 
mass  flowrate  is  applied  to  one  of  the  heat  exchangers), 
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and  bypass  flowrates  are  established  around  each  of  the 
other  heat  exchangers  such  that  the  dummy  heat  exchanger 
flowrates  are  maintained  at  their  calculated  values  and  con- 
tinuity of  the  total  loop  fluid  flowstream  is  insured.  (i.e.. 
Any  excess  fluid  over  that  required  to  supply  the  load  for  a 
given  heat  exchanger  is  shunted  around  that  heat  exchanger 
and  re-joins  the  total  loop  flow  downstream  of  the  heat  ex- 
changer outlet).  A series  of  subroutines  is  then  called 
which  calculate  the  loop  circulation  pump  pressures  corres- 
ponding to  these  desired  dummy  flowrates,  and  the  existing 
pump  pressures  are  adjusted  to  accelerate  or  decelerate  the 
observed  fluid  flowrates  over  a period  of  time  governed  by 
the  input  pump  response  sensitivity  factor.  Valve  models 
adjust  the  relative  heat  exchanger  flowrates  to  the  desired 
flow  split  conditions.  (See  Subroutine  PLOWS  for  further 
details  of  the  adjustment  process).  The  system  thus  responds 
to  load  variations  primarily  through  adjustment  of  the  loop 
fluid  flowrates,  allowing  constrained  fluid  temperature 
variations  to  provide  secondary  compensation  through  changes 
in* the  instantaneous  nodal  heat  storage  and  supply  rates. 

During  intervals  in  which  the  thermostat  is  not  activated. 
Program  MAIN  calls  the  dynamic  simulation  subroutines  which 
calculate  the  instantaneous  fluid  flows,  temperatures  and 
pressures  throughout  the  distribution  network.  A full  print- 
out of  these  TUS  parameters  and  the  consumer  load  distribution 
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is  given  at  regular  intervals  as  specified  by  the  output 
print  code. 

A. 3. 2 Input  Scaling  (Subroutine  CSCALE) 

Originally  conceived  as  a means  for  simplifying  the 
input  data  requirements  for  a series  of  simulations  inves- 
tigating the  effects  upon  system  operation  of  varying  only 
the  community's  thermal/electrical  energy  demand  ratio, 
this  subroutine  is  only  partially  completed.  The  use  of  an 
input  mass  flowrate  scale  factor  (CSFAC)  other  than  1.0 
will  result  in  the  specified  initial  flowrates  being  multi- 
plied by  that  factor  prior  to  their  use  in  the  TUS  initializa- 
tion calculations.  Since  neither  the  heat  exchanger  ratings 
nor  the  fluid  temperatures  are  scaled  internally,  it  Is 
necessary  to  re-punch  a significant  portion  of  the  input 
data  when  varying  the  demand  ratio  and,  to  avoid  errors,  it 
is  recommended  that  the  input  mass  flowrates  be  similarly 
repunched,  with  the  value  of  CSFAC  input  as  1.0  until  the 
internal  scaling  routine  has  been  completed. 

A. 3-3  Consumer  Heat  Flow  (Subroutine  QLOAD) 

Using  a combination  of  American  Society  of  Heating, 
Refrigeration  and  Air  Conditioning  Engineers  (ASHRAE)  recom- 
mendations [4]  and  simplifications  of  procedures  employed 
In  the  National  bureau  of  Standards  Load  Determination  code 
(NBSLD)  [5],  Subroutine  QLCAD  calculates  the  hourly  net 
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heat  losses  (gains)  from  each  building  type  in  the  community 
based  upon  a steady-state  energy  flow  analysis  including  the 
effects  of  heat  conduction  from  the  building  to  its  surroun- 
dings, infiltration  and  ventilation  air  flows,  solar  radia- 
tion absorption  and  transmission,  and  internal  heat  genera- 
tion from  lights,  appliances,  people,  etc.  Due  to  the  aggre- 
gate nature  of  the  building  units  chosen  to  typify  the  large 
numbers  of  individual  structures  classified  in  each  consumer 
grouping,  emphasis  in  the  development  of  Subroutine  QLOAD 
has  been  placed  upon  the  computation  of  building  energy 
demand  schedules  reflecting  the  major  components  of  the  total 
building  loads,  but  it  is  not  as  detailed  nor  as  dynamically 
rigorous  as  the  energy  transfer  function  analyses  performed 
by  NBSLD . A major  criterion  in  the  derivation  of  the  QLOAD 
demand  models  has  been  the  minimization  of  complex  building 
parameter  specifications  which  would  inordinately  restrict 
the  applicability  of  a given  building  unit  as  being  "typical" 
of  a fairly  diverse  mixture  of  structures.  However,  enough 
detail  is  included  in  the  calculations  that,  if  it  is  desired 
and  is  statistically  valid,  sensitivity  analyses  upon  the 


building  energy  demands  may  be  performed  by  varying  certain 
design  and  construction  specifications  on  a wall-by-wall 
basis . 


Some  basic  assumptions  made  in  the  demand  calculations 
are  that  all  buildings  specified  are  square  and  have  symmet- 


rical  walls  (i.e.,  that  each  of  a given  building's  four 
walls  is  constructed  identically  to  the  others);  that  the 
roofs  are  flat  (or  that’  the  input  data  describes  a horizontal 
projection  of  the  actual  roof);  that  building  occupancy, 
appliance  and  lighting  usage,  and  forced-air  ventilation 
requirements  are  closely  correlated  and  may. be  scheduled 
according  to  the  same  building  use  function;  and  that  a 
single  temperature  setting  is  to  be  maintained  throughout 
all  the  rooms  of  each  building  type.  All  loads  are  calcu- 
lated under  instantaneous  steady-state  heat  transfer  condi- 
tions, and  no  consideration  is  made  of  such  effects  as  the 
thermal  inertia  of  the  building  walls  or  of  energy  transfers 
among  regions  within  the  building. 

The  first  component  of  the  building  heat  loads  calculated 
in  Subroutine  QLOAD  is  that  due  to  the  combined  effects  of 
infiltration  and  ventilation  air  exchanges.  For  the  purposes 
of  this  analysis,  infiltration  is  defined  as  any  air  flow 
occurring  as  a result  of  wind  and  stack  effect  induced 
pressure  differentials  across  cracks  permeating  the  building's 
exterior  surfaces.  Ventilation  air  flows  are  those  air 
exchanges  scheduled  by  a building's  occupants,  either  by  the 
use  of  air  circulation  equipment  or  by  opening  windows. 
Infiltration  is  calculated  in  Subroutine  QLOAD  according  to 
the  carck  factor  method  described  in  the  User's  Manual  for 
NBSLD  [6].  In  general,  the  pressure  loading  on  a given  wall 
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consists  of  two  components:  that  due  to  differences  in  the 

temperature  and  density  of  the  air  between  the  base  and  top 
of  the  wall  (stack  effect  air  convection)  and  that  due  to 
wind  pressure  on  the  wall's  exterior  surface.  A simplified 
formula  for  determining  the  net  stack  effect  pressure  on  a 
wall  is  given  by  Eq.  (A.l). 

PC  » 5.3508  BH  (jig  - (S.l) 

The  conversion  factor  of  5-3508  includes  the  dimensional 
conversion  factor  of  7-644  (in.  I^O^K/ft)  multiplied  by  0.7, 
since,  according  to  ASHRAE,  the  net  effects  of  a distributed 
stack  effect  pressure  loading  may  be  approximated  by  a 
single  load  concentrated  at  10%  of  the  building's  height  [4]. 
The  theoretical  direct  normal  wind-inauced  pressure  loading 
on  a wall  is  given  by  Eq.  (A. 2), 

PV  = 0.000432(WV)2  (A. 2) 

2 

where  the  multiplier  of  0.000482  converts  (mph)  to  a pres- 
sure difference  in  inches  of  water,  relative  to  normal  atmos- 
pheric  pressure.  However,  since  the  ideal  conditions  under 
which  this  theoretical  loading  is  derived  are  seldom  observed 
in  practice,  the  pressure  computed  in  Eq.  (A. 2)  must  be 
scaled  to  account  for  the  actual  building  orientation  and  for 
empirically-determined  wall  surface  pressure  reduction  effects 
[6].  Figure  A. 5 illustrates  the  definition  of  the  relative 





wind  incidence  and  wall  normal  directions  as  they  are  used 
« * 

in  Subroutine  QLOAD.  Walls  are  numbered  clockwise  from 
the  North,  with  #1  being  the  reference  wall  for  normal  angle 
WA.  Because  the  air  flow  will,  in  general,  be  neither  normal 
nor  parallel  to  a given  wall,  the  normal  pressure  term  com- 
puted in  Eq.  (A. 2)  must  be  multiplied  by  normal  and  oblique 
air  flow  correction  factors  which  depend  upon  the  angle  of 
wind  incidence  relative  to  the  wall's  orientation  as  shown 
in  Table  A. 2.  Once  these  coefficients  have  been  determined 
for  a given  wall,  the  wind  pressure  on  that  wall  is  computed 
through 

PW  = (PTKN) (PTKO)PV  . (A. 3) 

It  should  be  noted  that  PW  may  be  wither  positive  or  nega- 
tive depending  upon  whether  the  wall  being  examined  is  on 
the  windward  or  leeward  side  of  the  building.  For  conserva- 
tism, it  is  assumed  that  the  stack  effect  and  wind  pressure 
loadings  on  a given  wall  add  to  produce  a total  pressure 
difference  of 

* " PT  = PW  + PC  . ( A . U ) 

Once  the  total  pressure  difference  across  a crack  is  known, 
the  volumetric  air  flowrate  through  the  crack  may  be  deter- 
mined through  the  use  of  infiltration  coefficients  based 
upon  ASHRAE  tabulated  air  flow  measurements  [6].  Table  A. 3 
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lists  the  air  flow  coefficients  - determined  for  several 
types  of  cracks  - which  are  used  in  the  general  infiltration 
formula 


(LI)i  = (CL)i(Ci)(PT)  1 . (A. 5) 

In  Subroutine  QLOAD,  infiltration  air  flowrates  are  computed 
in  this  manner  on  a wall-by-wall  basis  for  cracks  around 
doors,  windows  and  in  the  structural  walls,  and  the  resulting 
components  are  summed  to  produce  a net  total  infiltration 
air  leakage  for  the  building. 

For  many  types  of  small  buildings,  such  as  single  family 
homes,  infiltration  heat  losses  dominate  those  due  to  venti- 
lation. However,  for  large  commercial  buildings,  which  are 
generally  better  constructed  and  are  regulated  by  municipal 
building  codes,  ventilation  is  the  dominant  mode  of  air 
exchange.  In  Subroutine  QLOAD,  ventilation  is  scheduled 
according  to  each  building's  occupancy-related  use  factor  as 
is  showr.  in  Eq.  (A.6). 

VLEAK  = (BUF) ( AINFL)  (A.6) 

During  every  simulation  time  step,  each  building's  infiltra- 
tion and  ventilation  air  flowrates  are  compared,  and  the 
larger  is  used  to  determine  the  building's  instantaneous  air 
exchange  heat  loss  rate  through 
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TABLE  A. 3 


INFILTRATION  AIR  FLOW  COEFFICIENTS 


(from  Ref.  6) 


Double-hung  wooden  windows (locked)* 
non-weatherstripped  loose  fit 

. 6 

0 . 6)6 

average  fit 

2 

0. 66 

weatherstripped  loose  fit 

2 

0 . 66 

average  fit 

1 

0.66 

Window  frames* 

masonry  frame  with,  no  caulking 

1.2 

0 . 66 

masonry  frame  with  caulking 

0.2 

0 . 66 

wooden  frame 

1 

0.66 

Swinging  doors* 

1/2"  crack 

160 

0.5 

1/4"  crack 

80 

0.5 

1/8"  crack 

40 

0.5 

Walls** 

8"  plain  brick 

1 

0.8 

• 8"  brick  and  plaster 

0.01 

0.8 

13"  plain  brick 

0.8 

0.8 

13"  brick  and  plaster 

0.004 

0.7 

13"  brick,  furring,  lath  and  plaster 

0.03 

0.9 

frame  wall,  lath  and  plaster 

0.01 

0.55 

24"  shingles  on  1x6  boards  on 

9 

0 . 66 

14"  centers 

16"  shingles  on  1x4  boards  on 

5 

0 . 66 

5"  centers 

24"  shingles  on  shiplap 

3-6 

0.7 

16"  shingles  on  shiplap 

1.2 

0 . 66 

•Values  of  C listed  for  these  openings  are  per  linear  foot 
of  crack  length. 


••Values  of  C listed  for  the  walls  are  per  square  foot  of 
surface  area. 
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QLEAK  = (paCpa)  (ALEAK )TSB  - TAIR)  . (A. 7) 

(The  infiltration  and  ventilation  flowrates  are  not  summed, 
since  it  is  assumed  that,  to  some  extent,  one  substitutes  for 
the  other  in  determining  overall  building  comfort  levels. 

A comparison  between  the  two  components  allows,  for  example, 
ventilation  to  be  turned  off  at  night,  with  infiltration 
remaining  high  due  to  high  winds). 

In  order  to  determine  the  effects  of  solar  radiation 
absorption  and  transmission  upon  a building’s  total  space 
conditioning  energy  demands,  it  is  necessary  to  first  deter- 
mine the  direction  and  intensity  of  the  sun's  rays  striking 
the  four  walls  and  roof  of  the  building.  In  Subroutine  QLOAD, 
the  direction  cosines  of  a direct  solar  beam  are  referred  to 
the  normal  to  a horizontal  surface  as  is  shown  in  Fig.  A. 6. 

At  any  time  during  the  day,  the  solar  hour  angle  relative  to 
the  horizontal  surface  normal  is 

H = 15  (TIME  - 12)  . (A. 8) 

If  " 

| COS (H) | > |-TAN(LAT)*TAN(6) | 


the  sun  has  risen  above  the  horizon,  and  the  direct  solar 
beam  direction  cosines  are  defined  by 
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■ COS(Z)  =*  SIN(LAT) *SIN( 6 )+COS (LAT) *COS ( 6 ) *COS(H) , ( A.  9a) 
COS(W)  = COS (5 ) *SIN(H) , and  (A. 9b) 

COS(S)  = (1  - COS2(Z)  - COS2(W))0’5  , (A. 9c) 

where 


COS(S)  > 0,  if  COS (H)  > TAN (6 ) /TAN (LAT)  . 

(6  is  the  angle  of  solar  declination,  obtained  from  the 
monthly  data  in  Table  A.  4). 

The  intensity  of  the  direct  solar  beam  under  cloudless 
conditions  is  given  by 

DNSI  = ASC  * exp (-AEC/C0S ( Z) ) , (A. 10) 

where  ASC  and  AEC  are  the  Apparent  Solar  Constant  and  Atmos- 
pheric Extinction  Coefficient  obtained  from  the  monthly 
solar  incidence  data  summarized  in  Table  A. ■4.  The  angle  of 
incidence  of  this  direct  radiation  relative  to  the  horizontal 
normal  to  each  of  the  building's  four  walls  is  given  by 

COS(.Bi ) = SIN(SWAi)*COS(W)  + C0S(SWAi)*C0S(S),  (A. 11a) 

where  SWA  is  the  orientation  of  the  wall  being  studied  rela- 
tive to  South  (i.e,,  angle  WA  as  defined  in  Fig.  A. 5,  minus 
180  degrees).  The  cosine  of  the  angle  of  direct  solar  inci- 
dence upon  the  building's  (flat)  roof  is 


TABLE  A. 4 

SOLAR  INCIDENCE  DATA 
(from  Ref.  7] 


Apparent  Solar 
Constant 
(BTU/hr-ft2) 


Solar 


Atmospheric 


Sky 

Diffuse 

Factor 


Declination  Extinction 


(Degrees ) 


Coefficient 


Jan. 

21 

390 

-20.0 

0.142 

Feb  . 

21 

385 

-10.8 

0.144 

Mar. 

21 

376 

0.0 

0.156 

Apr. 

21 

360 

11.6 

0.180 

May 

21 

350 

20.0 

0.196 

June 

21 

345 

23.45 

0.205 

July 

21 

344 

20.6 

0.207 

Aug. 

21 

351 

12.3 

0.201 

Sept 

. 21 

365 

0.0 

0.177 

Oct . 

21 

378 

-10.5 

0.160 

Nov, 

21 

387 

-19.8 

0.149 

Dec  . 

21 

391  . 

-23.45 

0.142 

0.058 

0.060 

0.071 

0.097 

0.121 

0.134 

0.136 

0.122 

0.092 

0.073 

0.063 

0.057 


' 'Data  values  listed  are  for  the  Northern  Hemisphere;  for 
Southern  Hemisphere  locations,  the  values  of  AEC  and  SDF 
should  be  shifted  by  six  months. 
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= 0.45.  (A.Hb) 

The  total  Intensity  of  the  diffuse  radiation  striking  each 
of  the  building's  walls  is  thus 

DPSIi  = (Ida*Yi  + 0*5Idg)#PWI,IT1  »■  (A . l’5a ) 

and  of  that  striking  the  roof. 


\ 
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DFSIR  = Ida*FWLITr  . (A. 15b) 

Heat  gained  by  a building  from  the  sun  may  be  attributed 
to  two  effects:  the  internal  absorption  of  radiation  trans- 

mitted through  windows,  and  the  absorption  of  radiation  by 
the  walls  and  roof.  The  solar  heat  gain  through  a building's 
windows  is  obtained  through  the  calculation  of  a window 
solar  heat  gain  factor  based  upon  polynomial  expansions  of 
the  window's  solar  transmissivity  and  absorptivity  [9,10]. 

SHGF^  = DSIi*Fil  + 2*DFSI1*F12  + 0 .267  ( DSIi  *F13 

+ 2*DFSI1*FiA)  (A. 16) 

The  factors  F^-F^  are  the  transmission  and  absorption 
coefficients  obtained  from  the  expansions  for  1/8"  thick 

single  plate  glass: 


■ 


I 
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I 


. f • l T-COS^Ce^,  (A. 17a) 

Direct  transmission  Fj_i  j 1 

\ /M+2)  (A. 17b) 

Diffuse  transmission  Fi2  '^=QTj 

.,on  F = E A.CCS^e  ),  and  (A. 17c) 
Direct  absorption  ^3.3  . 0 J 1 

J 

- r A / ( 1+2)  (A.17d) 

Diffuse  absorption  j_q  j 

where  the  polynomial  coefficients  A.  and  ^ are  given  in 

Tab’e  a. 5.  The  factor  of  0.267  multiplying  the  radlatlo 

• c-q  (A  16)  arises  from  a comparison  of 
absorption  term  in  Eq.  ( 

, „ surface  thermal  resistances, 

typical  inner  and  outer  glass  surfac 

and  it  is  used  to  approximate  the  ratio  of  inward  to  outwar 

neat  flows  due  to  radiation  ahsorhed  hy  the  glass  [10]. 

The  rate  of  solar  heat  gain  through  the  windows  in  wall 

at  of  the  solar  heat  gain  factor  for  those 
given  by  the  product  of  the  f 

in  that  wall  (assumed  to  be  1/ 
windows,  the  window  area  in  that 

the  total  building  window  area  due  to  the  symmetric  1 a 

assumption),  and  an  Internal  window  shading  coefficient 

depending  upon  such  factors  as  double  glaaing,  tinted  glass, 
and  the  use  of  drapes,  blinds,  etc. 

aVIDSGj  - awd^shof^wdsc  (A’l8> 

rfMMents  are  listed  in  Table  A. 6 

Representative  shading  coefficients 


1 
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TABLE  A. 5 

POLYNOMIAL  COEFFICIENTS  USED  IN  CALCULATION  OF 
WINDOW  SOLAR  HEAT  GAIN  FACTORS 
(from  Ref.  9) 


Single  Glazing,  1/8"  sheet 


1 

ii 

Zi 

0 

0.01154 

-0.00885 

1 

0.77674 

2.71235 

2 

-3.94657 

-0.62062 

3 

8.57881 

-7.07329 

A 

-8.38135 

9.75995 

5 

3.01188 

-3.89922 

TABLE  A. 6 

WINDOW  SHADING  COEFFICIENTS 


WDSC 


TYPE  OF  GLASS 


Regular  sheet:  3/32  - 1/4 

0.95 

- 

1.00 

1/4  - 1/2 

0.88 

- 

0.95 

Grey  sheet : 1/8 

0.78 

0.80 

1/4 

0.86 

- 

0.88 

Heat-absorbing  plate:  1/4 

0.70 

— 

0.74 

1/2 

0.50 

0.57 

TYPE  OF  SHADING 

Venetian  blinds:  medium 

0.42 

_ 

0.64 

light 

0.40 

- 

0.55 

Roller  shade:  dark  opaque 

0.36 

— 

0.59 

white  opaque 

0.25 

- 

0.28 

translucent 

0.31 

- 

0.39 

Draperies 

0.35 

- 

0.80 

Louvered  sun  screens 

0.15 

- 

O 

-=T 

O 

for  a variety  of  common  glass  types  and  room  decors;  more 
extensive  tabulations  are  found  in  ASHRAE  [4],  pp . 402-408. 

The  individual-wall  solar  heat  gains  are  added  to  produce  a 
total  window  solar  heat  gain  for  the  building. 

The  effects  of  solar  radiation  absorption  and  reflection 
at  the  budding's  wall  and  roof  surfaces  may  be  approximated 
by  the  use  of  a surface  "sol-air"  temperature  in  the  building 
component's  steady-state  heat  conduction  equation  [4],  The 
sol-air  temperature  for  a surface  is  given  by 

TSAj^  = TAIR  + SABi*(DSIi  + DFSIi)*H01  , (A. 19) 

where  SAB^^  is  the  solar  absorptivity  of  the  surface  and  HC^ 
is  the  surface  convection  heat  transfer  coefficient.  Table 
A. 7 lists  several  absorptivities  for  typical  building  materi- 
als, a more  extensive  summary  of  which  is  given  in  Ta'^le  C6 
of  NBSLD  [5].  The  convection  heat  transfer  coefficient  for  a 
surface  depends  primarily  upon  the  velocity  of  the  wind  strik- 
ing the  surface  and  the  roughness  of  the  surface  material: 

HO,  = 5 , (A. 20) 

Cai*(WV)  + cbi*(wv)  + Cci 

where  the  coefficients  in  this  formula  are  functions  of  the 
surface  material  type  as  given  in  Table  A. 8.  To  determine 
the  net  heat  conduction  rate  out  of  the  building  through  its 
walls  and  roof,  including  the  effects  of  solar  heating,  the 
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TABLE  A. 7 

SOLAR  A3S0RFTIVTTIE5  OF  TYPICAL  3UILDING  MATERIALS 

(from  Ref.  5) 


MATERIAL 

SAB 

Brick : 

red  common 

0 . 68 

light  buff 

0.36 

white  glazed 

0.26 

Tiles : 

dark  clay 

0.82 

red 

0.67 

brown  concrete 

0.85 

Hoofing 

asphalt 

0. 9-l 

bituminous  felt 

0.89 

black  matte  sheet 

0.97 

blue  slate 

0.86 

Wood 

0.78 

Marble 

0.58 

White  plaster 

0.29 
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TABLE  A. 8 

COEFFICIENTS  USED  IN  CALCULATING  SURFACE 


CONVECTION 

HEAT  TRANSFER 

COEFFICIENTS 

MATERIAL 

(from  Ref.  5) 
Ca 

Cb 

C 

c 

Stucco 

0.0 

0.464 

2.04 

Brick,  rough  plaster 

0.001 

0.320 

2.20 

Concrete 

0.0 

0.330 

1.90 

Clear  pine 

-0.002 

0.315 

1.45 

Smooth  plaster 

0.0 

0.244 

1 . 80 

Glass,  white  paint 
on  pine 

-0.00125 

0.262 

1.45 
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wall  and  roof  material  heat  transfer  coefficients  are  cal- 
culated by 


AW  AWD, 

AUW  = = + ±_ 

*uwi  RTW1 + WHOi  RTWD1 


(A. 21a) 


and 


AUR 


AR 

RTR  + RHO 


(A. 21b) 


The  total  heat  conduction  rate  to  the  ambient  air  is  obtained 
by  adding  the  components  for  the  four  walls  and  the  roof: 

4 

QCOND  = I AUW. * (TSB  - TSA, ) + AUR* ( TSB  - TSAR)  .(A. 22) 
a i-1  1 1 

The  heat  conducted  from  the  building  to  the  earth  is 

QCOND  = AUG*  (TSB  - gS3 +JTMB)  ^ (A.  23) 

O ^ 

where  the  average  between  the  desired  building  temperature 
and  the  ambient  earth  temperature  is  used  to  account  for 
heating  of  the  soil  in  direct  contact  with  the  building. 

The  total  heat  conduction  rate  from  the  building  is  given 
by  the  sum  of  this  ground  loss  and  the  heat  lost  to  the  air. 

Energy  flows  due  to  the  production  of  heat  within  a 
building  are  lumped  into  a single  term  accounting  for  light- 
ing, appliances,  people,  etc.,  which  are  scheduled  according 
to  the  building’s  occupancy-related  use  factor  as 
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1 

I 

1 


l.  i 


QLIG  = (BUF)CPOL)  . (A. 24) 


Finally,  the  net  total  heat  loss  from  a building  due  to 
these  component  energy  flows  is  given  by 

QL  = QCOND  + QLEAK  - QWDSG  - QLIG  . (A. 25) 

It  should  be  noted  that,  as  defined  here,  a positive  heat  flow 
indicates  a net  heat  loss  from  the  building  to  its  surroun- 
dings . 

A . 3 • 4 Building  Load  Distribution  (Subroutine  DLOAD) 

In  Subroutine  DLOAD,  the  building  heat  losses  (gains) 
computed  in  Subroutine  QLOAD  are  applied  to  the  TUS  load 
center  heat  exchangers  and  the  electrical  supply  network 
according  to  the  specified  distributions  of  buildings  and 
end-use  equipment  in  the  community  being  studied.  When 
translating  the  heat  losses  to  utility  system  demands,  care 
is  taken  to  include  the  effects  of  variations  In  the  end-use 
equipment  coefficients  of  performance  with  changes  in  the 
building  room  and  ambient  air  temperatures.  A printout  of 
the  community  load  distribution  is  provided  at  regular  inter- 
vals during  the  simulation  as  specified  by  the  user-defined 
output  code.  Since  the  demands  in  Subroutine  QLOAD  are  cal- 
culated only  at  hourly  intervals.  Subroutine  DLOAD  computes 
hourly  load  distributions  which  are  used  as  inputs  to  linear 


. 
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interpolation  equations  in  Program  MAIN  to  produce  the  con- 
tinuous load  schedules  required  for  the  dynamic  system  analy- 
sis. (Instantaneous  distributions  may  be  calculated  by  Sub- 
routine DLOAD  at  any  time  during  the  simulation  as  called 
for  by  the  output  schedule,  but  these  instantaneous  values 
are  based  upon  linear  interpolations  between  previously  com- 
puted hourly  data  points.) 

Associated  with  each  load  center  heat  exchanger  in  the 
TUS  are  a set  of  six  numbers  for  each  building  type  in  the 
community:  the  number  of  units  of  that  building  type  served 

by  the  heat  exchanger,  and  the  numbers  of  compressive  air 
conditioners,  absorptive  air  conditioners,  heat  pumps,  hot 
water  space  heating  units  and  resistance  heaters  within  these 
buildings.  (e.g.,  load  center  heat  exchanger  #5  might  serve 
20  units  of  building  type  1,  of  which  11  have  absorptive  air 
conditioners,  9 have  compressive  air  conditioners,  7 have 
heat  pumps,  13  have  hot  water  heat,  and  0 have  resistance 
heaters.)  At  each  hourly  interval,  all  the  specified  build- 
ing types  are  examined  at  each  heat  exchanger  to  determine 
how  many,  if  any,  of  each  type  have  equipment  connected  to 
the  TUS  at  that  point.  The  total  space-conditioning  heat 
loss  from  these  buildings  is  then  calculated  by 

VI. 

= NU*  * QL±  , (A. 26) 

where  QL^  is  the  net  heat  loss  (gain)  from  one  unit  of  build- 
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ing  type  i computed  in  Subroutine  QLOAD  (or  is  an  off-hour 
load  obtained  through  linear  interpolation  between  two 
hourly  values),  and  the  superscript  indicates  the  calcula- 
tion is  performed  for  buildings  served  by  heat  exchanger  k. 

Domestic  hot  water  energy  demands  are  also  computed  in  Sub- 

*• 

routine  DLOAD  for  each  of  the  building  types  according  to 
the  user-specified  water  consumption  schedules  through 

QDHW*  = NU^  * WUF1  * DHW1  . (A. 27) 

(The  service  hot  water  demands  are  computed  here  rather  than 
in  Subroutine  QLOAD  as  a matter  of  convenience.  They  are 
added  to  the  net  space  conditioning  loads  and  cannot  be  used 
to  offset  building  heat  losses  computed  in  QLOAD) . 

After  computing  the  aggregate  space  conditioning  demand 
for  a given  building  type  at  a given  load  center  heat  exchan- 
ger, Subroutine  DLOAD  determines  whether  this  demand  is  due 
to  space  heating  or  cooling  and  applies  it  to  the  appropriate 
end-use  equipment  for  that  building  type.  If  the  load  is  due 
to  space  cooling,  it  may  be  supplied  by  either  electric  com- 
pression air  conditioners  or  by  hot  water  (or  low  pressure 

* 

steam)  liquid  absorption  chillers,  i£  the  outside  air  temp- 
erature is  greater  than  that  desired  within  the  building. 

If  the  outside  air  temperature  is  lower  than  the  desired 
roorr.  temperature  (the  cooling  load  resulting  from  internal 
heat  generation  or  solar  heating),  it  is  assumed  that  venti- 
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lation  will  be  increased  to  cool  the  building,  and  no  loads 
are  applied  to  either  type  of  air  conditioning  equipment. 
(This  criterion  may  be  employed  by  the  user  to  insure  cer- 
tain buildings  of  not  being  cooled  by  simply  specifying  their 
desired  temperatures  to  vary  exactly  with  the  known  ambient 
air  temperature  schedule).  If  the  outside  air  temperature 
exceeds  the  room  temperature,  the  nominal  coefficient  of  per- 
formance of  each  air  conditioning  unit  (defined  at  an  air 
temperature  of  90  °F  and  a room  temperature  of  75  °F)  is 
scaled  according  to  the  theoretical  Carnot  cycle  performance 
for  a refrigerator: 


CO?’ 


COP 

o 


,0.028037s 
tOAB  , ’ 
TIAB 


(A. 28) 


where  the  0.028037  factor  is  defined  by  the  ideal  Carnot 
heat  engine  efficiency  at  the  specification  temperatures: 


cops 


TOAB  , _ 90+460  , 

TIAB  " 75+%0  _1 


0.028037  • 


It  should  be  noted  that  COP  , the  nominal  coefficient  of 

„ o 

performance  for  the  given  air  conditioning  unit  at  TOAB  = 

550  °K  and  TIAB  = 535  °K,  should  not  be  the  theoretical 
Carnot  cycle  COP  under  those  conditions,  but  should  reflect 
the  manufacturer's  specifications  for  the  actual  units  in- 
stalled in  the  system.  The  Carnot  cycle  efficiency  is  merely 


used  by  Subroutine  DLOAD  as  an  approximate  indicator  of  the 


relative  behavior  of  the  actual  COP  as  the  room  and  air 


temperatures  vary  over  time.  Applying  the  COP'  computed  for 


each  of  the  air  conditioning  units  to  the  aggregate  loads  to 


be  met  by  those  units,  the  total  energy  demands  for  space 


cooling  are  determined  by 


electric  compression:  QCMP 


*Q,  /bcm: 


thermal  absorption:  QABS 


where  F 


is  the  ratio  of  electrically  supplied  space  cooling 


to  the  total  space  cooling  loads  for  building  type  i at  heat 


NCMP 


NCMP 


+ NABS 


( BCMP  and  BABS  are  the  scaled  COP's  for  the  compressive  and 


absorptive  air  conditioning  units,  respectively). 

Space  heating  demands  are  treated  analagously  to  cooling 


room  temperature  for  a given  building  type,  no  space  heating 


demand  for  that  building  type  is  applied  to  the  energy  dis 


tribution  networks.  Heat  pump  COP's  are  defined  at  conditions 


air  temperature  and  68  °F  room  temperature,  with  the 


deal  Carnot  cycle  efficiency  factor  being 
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■,  TOAE  _ , 32  + i46C  _ „ neo,an 
1 " TTa^  " 1 - 53T2F0  ' °-068182' 


Variations  in  either  temperature  from  these  set-point  con- 
ditions cause  the  COP  to  be  modified  according  to 


cop,  , COPo  (SU168182,  . 

1 “ TIAB 


(A. 31) 


The  space  heating  demands  may  be  supplied  by  heat  pumps,  hot 
water  heat.,  or  by  electric  resistance  heaters,  with  the  loads 
applied  to  the  utility  systems  being  determined  by 


electric  heat  pumps: 

QHPM^  = 

FhPi*«£/BHPM> 

(A. 32a) 

hot  water: 

QHWS^  = 

Ph»i*Ql/EHWS-  and 

(A. 32b) 

resistance : 

QRES^  = 

Fri*Qi/ERES> 

(A. 32c) 

where  the  load  fractions  for  building  type  i at  heat  exchanger 
k are  defined,  in  general,  by 


xi  NHPmJ  + NHWsi1  + NRES^ 
i i i 


(A. 33) 


All  heat  pumps  in  the  system  are  assumed  to  be  provided  with 
auxiliary  resistance  heating  coils.  At  each  hourly  interval, 
the  scaled  heat  pump  COP  (BHPM)  is  compared  with  the  COP  of 
a resistance  heating  unit  (ERE S,  assumed  to  be  1.0),  and  if 
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the  heat  pump  COP  is  less  than  that  of  the  resistance  heater, 
all  heat  pump  loads  are  transferred  to  their  resistance 
back-up  units. 

Domestic  hot  water  demands  are  assumed  to  be  supplied 
through  the  use  of  hot  water  heat  exchangers  in  buildings 
supplied  with  hot  water  heat  from  the  TUS,  or  through  the 
use  of  electric  hot  water  heaters  in  buildings  containing 
heat  pumps  or  resistance  space  heaters.  The  service  hot 
water  demands  represent  an  additional  heating  load  which  may 
be  conceptually  applied  to  the  same  type  of  equipment  as 
that  supplying  the  space  heating  requirements.  Therefore, 
for  buildings  with  domestic  hot  water  loads,  the  equipment 
energy  demands  computed  on  Eqs.  (A.32b,c)  are  modified  to 
yield 

hot  water:  QHWS^  = QHWS^( space  heating)  + 

k 1 k (A. 34a) 

Fhwi*QDHWi/EKWS’  and 

k k 

resistance:  QRES.  = QRES. (space  heating)  + 

k k k (A. 34b) 

(Fhpi  ri } *QDHWi/ERES . 

This  load  analysis  is  performed  for  every  building  type 
at  each  load  center  heat  exchanger  for  each  hour  of  the 
simulation  period.  The  resulting  thermal  demands  define  the 
distributed  hourly  load  schedules  to  be  met  by  the  TUS,  and 
the  electrical  demands  are  aggregated  to  form  one  component 
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of  the  community's  electrical  load  schedule.  In  addition 
to  the  buildings  served  by  the  TUS  and  distributed  among  its 
load  center  heat  exchangers,  there  may  be  several  buildings 
in  the  community  which  are  not  provided  with  HTW  service 
but  which  are  supplied  with  electricity  from  the  TES  power 
plant.  For  the  purposes  of  calculating  these  "external 
system"  loads  in  Subroutine  DLOAD,  all  of  the  units  of  a 
given  building  type  not  served  by  the  TUS  are  aggregated  into 
a single  grouping,  independently  of  their  actual  geographical 
locations  within  the  community.  (Since  no  detailed  analysis 
is  performed  on  the  electrical  distribution  network,  there 
is  no  need  to  specify  transmission  distances  to  these  con- 
sumers.) Proceeding  analogously  to  the  computations  described 
above  for  those  buildings  connected  to  cr.e  TUS,  the  external 
system  energy  demands  are  calculated  as  follows. 

1)  Total  space  conditioning  demand  for  all  units  of  build- 
ing type  i not  servied  by  the  TUS: 

QEXT±  = NUEXi*QL1  , (A. 35) 

where  QL.  is  the  hourly  demand  computed  in  Subroutine 

* i 

QLOAD,  or  is  an  interpolated  load  if  the  demands  are 
being  computed  during  one  of  the  simulation  hours. 

2)  Total  domestic  hot  water  energy  demand  for  all  units  of 
building  type  i not  served  by  the  TUS: 

QDEXTi  = NUEX1*WUFi*DHW1  . (A. 36) 
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3)  Total  compressive  air  conditioning  energy  demand  for 
all  units  of  building  type  i not  served  by  the  TUS : 

QCEXT1  = QEXT  /BCMP  , (A. 37) 

where  QCEXT.^  is  defined  for  TOAB  > TIAB^,  and  BCMP  is 
the  air  conditioner  COP  scaled  according  to  the  rela- 
tion in  Eq . (A. 28).  It  should  be  noted  that  all  build- 
ings in  the  external  system  requiring  cooling  are 
assumed  to  be  served  by  compressive  air  conditioning 
units.  Those  few  building  not  served  by  the  TUS  using 
absorption  chillers  will  also  have  in-house  fossil- 
fired  furnaces  to  supply  their  heating  loads  - due  to 
the  high  cost  of  installing  and  operating  dual  heating 
and  cooling  systems  - and,  for  the  purposes  of  this  TES 
analysis,  will  have  no  effect  upon  the  overall  energy 
system  design. 

4)  Total  heat  pump  energy  demand  for  all  units  of  building 
type  i not  served  by  the  TUS: 

QHEXT^  = FHPEX1 *QEXT^/BHPM  , (A. 38) 

where  all  space  heating  demands  are  defined  only  for 
TOAB  < TIAB^,  and  QHEXT^  is  zero  if  the  scaled  COP 
(BHPM,  from  Eq . (A. 31))  is  less  than  1.0. 

Total  resistance  heating  energy  demand  .for  all  units  of 
building  type  i not  served  by  the  TUS: 


5) 


56 . 

QREXT.  = FREXT,  “QEXI^/ERES  + 

(FREXT.  + FHPEXi)*QDEXT1/ERES  , (A. 39) 

and  QREXTi  includes  the  heat  pump  space  heating  demands 
if  the  heat  pump  COP  is  less  than  ERES. 

Total  hot  water  energy  demand  for  all  units  of  building 
type  i not  served  by  the  TUS : 

QHWEXi  = (1  - FHPEX.  - FREXTi)*(3EXTi  + QDEXTi)/EHWS  . 

(A. 40) 

This  demand  is  assumed  to  be  supplied  by  in-house  fossil 
fuel  furnaces  and  is  not  applied  to  the  TES  analysis 
(other  than  as  an  entry  in  the  external  system  building 
load  distribution  summary).  Although,  as  is  described 
above,  it  is  assumed  that  any  external  system  buildings 
containing  absorptive  air  conditioning  units  will  also 
contain  fossil-fired  furnaces  and  will  not  affect  the 
TES  analysis,  present  construction  practices  demonstrate 
that  it  is  very  likely  that  a building  having  an  in- 
house  fossil  furnace  will  also  have  compressive  air 
conditioning.  These  buildings  will  affect  the  total 
electrical  loads  to  be  met  by  the  TEC  power  : '.ant,  and 
they  must  therefore  be  included  in  this  d<j~  an  : a-  . y Is. 
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In  order  to  compute  the  total  system  hourly  electrical 
energy  demands,  the  electrical  demands  computed  in  Eqs . 

(A. 29a),  (A. 32a),  and  (A.3^b),  summed  over  all  the  building 
types  and  all  the  load  center  heat  exchangers,  are  added  to 
the  total  external  system  electrical  demands.  These  hourly 
space  conditioning  and  domestic  hot  water  loads  are  trans- 
mitted to  Program  MAIN,  where  they  are  added  to  the  aggre- 
gate non-space-conditioning  electrical  load  for  the  commu- 
nity to  produce  the  hourly  total  electrical  demands  to  be 
met  by  the  TES  power  plant  electrical  output. 

A. 3- 5 Steady-State  System  (Subroutine  STDY) 

Subroutine  STDY  is  an  iterative  steady-state  heat  flow 
subprogram  which  establishes  the  nodal  fluid  flowrates  and 
temperatures  throughout  the  TUC  for  given  thermal  load  and 
supply  conditions.  The  essence  of  the  routine  lies  in  a 
steady-state  heat  balance  performed  on  each  of  the  nodes 
which  requires  equality  between  the  convective  and  conductive 
heat  flows  from  the  node.  The  further  requirements  of  con- 
tinuity of  the  temperature  profiles  and  fluid  flows  across 
the  boundaries  between  consecutive  nodes  within  each  loop 
serve  to  physically  Join  the  separate  nodal  parameters  and  to 
insure  conservation  of  energy  and  momentum  around  the  loops. 
The  loop  flows  and  temperatures  resulting  from  the  heat 
balance  are  adjusted  iteratively  until  the  system  converges 
to  the  conditions  satisfying  both  nodal  and  total  system 
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energy  balance. 

« 

The  bases  of  the  nodal  heat  balance  calculation  are  the 
steady-state  convection  and  conduction  equations  given  by: 


convection:  Q = MASS*CAP* (TI  - TO)  , and 


(A. 41) 


conduction:  Q = AUF*DTLM  . 


(A. 42) 


DTLM  is  the  log-mean  temperature  difference  between  the 
fluid  in  the  node  and  the  fluid  to  which  heat  is  being 


transferred  and  is  expressed  by: 


DTLM  = 


TN  - TD 

\ipD' 


(A. 43) 


where  TN  = temperature  difference  across  the  node  inlet,  and 
TD  = temperature  difference  across  the  node  outlet. 


For  a length  of  pipe  transferring  heat  to  the  earth. 


Eq.  (A. 42)  becomes: 


Q = (AUF)  * 


TI  - TO 
, TI  - TAMB\ 
1 ''TO  - TAMB' 


(A. 44) 


since  the  earth  temperature  is  assumed  to  be  constant  along 
the  length  of  the  pipe.  Combining  Eqs.  (A. 41)  and  (A. 44) 
yields  the  following  expressions  for  the  nodal  inlet  and 


outlet  fluid  temperatures  for  a length  of  pipe: 


TI  = TO  exp  [AUF/(MASS‘:fCAP)]  + TAMB  [l  - exp[AUF/(MASS*CAP  ) ]]  , 


(A. 45a) 


TO  = TI  exp  [-AUF/( MASS# CAP ) ] + TAMB|l-  exp [-AUF/ (MASS# CAP ) ]]  , 


(A. 45b) 


All  heat  exchangers  are  assumed  to  be  of  the  single-pass 
counterflow  type#.  Designating  the  node  being  analyzed  as 
the  "1"  or  "primary"  side  of  the  heat  exchanger  and  the  oppo- 
site node  as  the  "2"  or  "secondary"  side,#*  the  conductive 
heat  flow  out  of  the  primary  side  can  be  expressed  by: 

(see  Fig.  A. 7 for  reference) 


= ( AUF) 


(TI1  - T02)  - (T01  - TI2) 
* TI  - TO- 

ln  (T01  - TI2} 


(A. 46) 


#If  it  is  desired  to  use  a different  configuration  in  the 
actual  TUS  design,  a conversion  factor  may  be  applied  to  the 
heat  transfer  coefficient  of  the  actual  heat  exchanger  to 
transform  it  to  an  equivalent  counterflow  heat  exchanger  for 
the  purposes  of  this  model.  [11] 

**In  this  context,  "primary"  and  "secondary"  are  not  neces- 
sarily synonymous  with  hot  and  cold.  The  analysis  developed 
here  is  independent  of  the  side  of  the  heat  exchanger  being 
studied;  If  the  "1"  side  is  the  cold  side,  the  heat  transfer 
out  of  the  node  will  simply  be  negative. 
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The  steady-state  convective  heat  transfer  rates  on  both 
sides  of  this  heat  exchanger  must  be  equal  to  the  conductive 
heat  flowrate  between  them  and  are  given  by: 

Q = MASS1*CAP1*(TI1  - T01),  and  (A. 47a) 

Q = MASS2*CAP2«(T02  - TI2)  . (A.*47b) 

Combining  Eqs.  (A. 46),  (A. 47a),  and  (A. 47b)  to  eliminate  TO^, 
the  primary  nodal  inlet  temperature  is  determined  by  the 
relation : 


TI, 


mass2*cap2 

MASS1*CAP1 


-1 


(to2  - ti2) 


1 - exp  AUF(MASs2»CAp2  ‘ mAssi*CAP1) 


+ T02  (A. 48) 


hd 


A closed  solution  for  all  the  nodal  fluid  mass  flowrates 
and  temperatures  can  be  obtained  only  when  information  is 
available  regarding  the  total  thermal  demands  at  each  load 
center  heat  exchanger  and  the  TUS  steady-state  heat  input. 

In  practice,  Subroutine  STDY  requires  as  input  data  the  inlet 
and  outlet  temperatures  and  fluid  flowrates  on  the  consumer 
sides  of  all  the  load  center  heat  exchangers,  the  primary 
HTW  loop  fluid  flowrate,  and  the  primary  loop  fluid  supply 
temperature  at  the  storage  reservoir.  (See  Fig.  A. 2 for  a 
schematic  representation  of  the  TUS.)  Subroutines  QLOAD  and 
DLOAD  determine  the  steady-state  loads  to  be  applied  to  each 
load  center  heat  exchanger,  and  the  required  fluid  flowrates 
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are  calculated  to  supply  these  loads  within  the  user-defined 
heat  exchanger  temperature  limits.  The  primary  loop  steady- 
state  energy  supply  data  are  specified  by  the  user  according 
to  the  desired  system  design.  Estimates  of  the  total  loop 
fluid  flowrate  and  of  the  load  center  heat  exchanger  parallel 
branch  flow  splits  in  each  of  the  secondary  distribution  loops 
may  also  be  supplied  to  Subroutine  STDY  by  the  user,  and  will 
speed  the  convergence  of  the  iterative  solutions  if  chosen 
properly . 

The  system  calculations  are  divided  into  two  sections 
corresponding  to  the  two  sets  of  heat  exchangers  at  the 
secondary  loop  and  consumer  supplies.  For  each  secondary 
loop.  Subroutine  STDY  first  calculates  a load  center  heat 
exchanger  supply  header  reference  temperature  through  the 
use  of  Eq.  (A. 48),  the  computed  load  center  heat  exchanger 
consumer-side  demand  data,  and  an  estimated  value  of  the  first 
secondary  loop  parallel  branch  mass  flowrate.  Since  all  the 
load  center  heat  exchangers  are  supplied  from  a common  point 
in  the  secondary  loop,  the  inlet  temperatures  on  their  supply 
sides"  must  all  be  equal  to  the  calculated  reference  tempera- 
ture. Therefore,  trial  inlet  temperatures  computed  through 
Eq.  (A. 48)  for  each  of  the  load  center  heat  exchangers  are 
compared  with  the  reference,  and  the  branch  flowrates  are 
adjusted  until  the  calculated  temperatures  are  within  0.001? 
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of  the  required  value.*  When  the  temperatures  have  con- 
verged, they  are  set  exactly  equal  to  the  reference  value  to 
avoid  instabilities  in  future  iterations,  and  the  nodal 
outlet  temperatures  are  calculated  by  requiring  equality 
between  the  convective  heat  flows  on  both  sides  of  the  heat 
exchangers.  Upon  completion  of  this  iterative  process  for 
each  of  the  load  center  heat  exchangers  served  by  the  secon- 
dary loop,  the  parallel  branch  fluid  flowrates  are  added  to 
produce  the  total  secondary  loop  flow,  and  the  parallel  nodal 
outlet  temperatures  are  combined  to  yield  a mass-flow-averaged 
inlet  temperature  to  the  secondary  loop  return  pipe.  The  out- 
let temperature  on  the  secondary  loop  side  of  the  secondary 
supply  heat  exchanger  is  then  calculated  from  Eq . (A . 45a ) using 
the  parallel  branch  supply  header  reference  temperature  as 
TO;  the  inlet  temperature  to  the  secondary  supply  heaL.  exchan- 
ger is  calculated  using  Eq.  (A.^5b)  and  the  mass-flow-averaged 
return  header  temperature. 

The  inlet  temperature  on  the  primary  loop  side  of  the 
first  secondary  supply  heat  exchanger  is  computed  by  using 
in  Eq . (A. 45b)  the  specified  loop  fluid  flowrate  and  the 
HTW  supply  temperature  at  the  storage  reservoir.  Since  the 


*Due  to  round-off  errors  in  the  calculations,  the  two  values 
will  never  become  exactly  equal;  a higher  degree  of  precision 
than  the  0.001$  chosen  leads  to  much  longer  times  for  final 
system  convergence. 


primary  HTW  loop  fluid  flowrate  and  supply  temperature  are 
fixed  by  the  system  design,  this  inlet  temperature  is  con- 
stant and  defines  the  HTW  reference  temperature  for  the 
first  secondary  supply  heat  exchanger.  The  loop  fluid  flow- 
rate  and  the  inlet  and  outlet  temperatures  on  the  secondary 
side  of  the  heat  exchanger  calculated  in  the  iterative  ana- 
lysis of  the  load  center  heat  exchangers  are  then  used  in 
Eq.  (A. 48)  to  yield  a trial  primary-side  supply  temperature 
which  is  compared  with  the  reference.  If  the  difference 
between  these  values  is  greater  than  the  0.001%  convergence 
limit,  the  secondary  loop  fluid  flowrate  is  adjusted  to  de- 
crease the  error.  The  new  secondary  loop  flow  is  split 
among  the  parallel  load  center  heat  exchangers  in  proportion 
to  their  loads,  a new  supply  header  reference  temperature  is 
calculated,  and  the  load  center  heat  exchanger  iteration 
process  is  begun  again. 

This  double  iteration  procedure  continues  until  all  the 
calculated  temperatures  for  the  first  secondary  distribution 
loop  fall  within  the  required  convergence  limits.  The  outlet 
temperature  on  the  primary  loop  side  of  the  secondary  supply 
heat  exchanger  is  then  set  by  requiring  the  specified  primary 
1 op  fluid  mass  flowrate  to  supply  the  total  heat  load,  includ- 
ing losses,  of  the  secondary  loop  being  examined.  A reference 
HTW  inlet  temperature  for  the  next  secondary  supply  heat 
exchanger  is  calculated  using  Eq . (A. 45b)  and  the  first  heat 


exchanger  outlet  temperature,  and  the  entire  loop  iteration 
procedure  is  started  for  the  next  secondary  distribution 
loop.  When  all  the  loops  have  converged,  the  steady-state 
heat  flows  across  each  of  the  .heat  exchangers  are  calculated 

and  the  nodal  fluid  mass  flowrates  and  temperature  data  are 

* 

printed  out  in  a TUS  steady-state  design  summary. 


A. 3. 6 Initialization  (Subroutine  SET) 


Subroutine  SET  is  an  intermediate  program  flow  control 
subroutine  which  is  called  to  control  the  building  load 
analysis  section  of  TDIST,  to  initialize  the  dynamic  system 
simulation  parameters  to  the  converged  steady-state  condi- 
tions computed  in  Subroutine  STDY,  or  to  control  the  program 
flow  .when  the  Program  MAIN  thermostat  model  requires  adjust- 
ment of  the  loop  fluid  flowrates. 

During  system  initialization  or  steady-state  design  ana- 
lyses (KODE  = 0),  Subroutine  SET  assumes  control  from  Program 
MAIN  and  causes  the  hourly  (or  single  time  instant)  building 
energy  demands  and  the  resulting  TUS  load  distributions  to 
be  computed  by  Subroutines  QLOAD  and  DLOAD.  If  only  an  analy 
sis  of  the  building  loads  is  desired  (IDES  = 1),  Subroutine 
SET  causes  the  building  energy  demand  schedules  and  the 
desired  load  distribution  to  be  printed  out  and  then  returns 
control  to  Program  MAIN,  from  whence  the  design  run  is  termin- 


If  it  is  desired  to  use  the  loads  computed  at  the  begin- 
ning hour  of  the  simulation  period  as  input  data  to  the 
steady-state  TUS  calculations  performed  in  Subroutine  3TDY 
(IDES  = 0),  Subroutine  SET  uses  the  lead  distribution  from 
Subroutine  DLOAD  and  the  user-specified  load  center  heat 
exchanger  supply  and  return  temperatures  to  compute  the 
fluid  mass  flowrates  on  the  consumer  sides  of  all  the  load 
center  heat  exchangers.  To  prevent  destruction  of  the  exis- 
ting system  temperature  and  flowrate  data  in  the  event  of 
non-convergence  in  Subroutine  STDY,  a set  of  dummy  flowrates 
and  temperatures  are  initialized  to  the  existing  system  condi- 
tions and  are  transmitted  to  Subroutine  STDY  as  the  input 
data  required  by  the  heat  balance  analyses.  Upon  return  of 
the  converged  steady-state  TUS  parameters.  Subroutine  SET 
transmits  the  loop  fluid  flowrates  to  Subroutine  PRESS,  in 
which  the  loop  fluid  frictional  pressure  losses  and  required 
circulation  pump  pressure  settings  are  computed.  Initializa- 
tion of  the  TUS  nodes  to  the  computed  steady-state  conditions 
begins  with  the  resetting  of  the  actual  nodal  temperatures 
and  fluid  mass  flowrates  to  the  values  returned  from  Sub- 


routine STDY.  The  temperature  of  the  cold  water  section  of 
the  thermal  energy  storage  reservoir  (see  Fig.  k.h)  is  set 
equal  to  the  fluid  temperature  at  the  outlet  of  the  last 
pipe  node  in  the  primary  HTW  loop;  the  hot  water  section 
temperature  is  fixed  at  the  desired  system  HTW  supply  tempera- 
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ture  specified  by  the  user,  and  the  fluid  mass  flowrate  on 
the  TUS  side  of  the  power  plant  thermal  energy  supply  heat 
exchanger  is  calculated  so  as  to  absorb  the  net  rate  of 
thermal  energy  production  by  the  power  plant  after  meeting 
the  steady-state  system  electrical  demands.  (The  interfacing 
heat  exchanger  between  the  power  plant  and  the  TUS  is  not 
modelled  in  detail  by  TDIST  and  is  simply  assumed  to  have  an 
"extremely  large"  heat  transfer  coefficient.  It  is  included 
as  a separate  item  in  the  system  energy  flow  analysis  to 
emphasize  its  existence  and  to  provide  a mathematical  means 
for  transferring  energy  between  two  physically  different 
fluid  flowstreams ) . Following  the  computation  of  this  TUS 
supply  fluid  flowrate,  Subroutine  PREC,  which  monitors  the 
TUS  supply  heat  exchanger  to  insure  that  the  Second  Law  of 
Thermodynamics  is  always  obeyed  by  its  combined  energy  flows, 
is  called  to  verify  the  validity  of  the  system  supply  condi- 
tions. Subroutine  SET  then  calculates  measures  of  the  energy 
st  red  in  each  of  the  reservoir  sections,  based  upon  their 

i 1 volumes  and  fluid  temperatures.  (In  the  present  form 
of  the  code,  these  are  not  true  measures  of  the  actual  BTUs 
of  energy  stored,  since  the  computations  are  performed  in 
degrees  Fahrenheit,  referred  tr  0 °F  rather  than  0 °K. 

As  a result,  the  printout  data  summarizing  the  reservoir 
energy  storage  should  only  be  used  as  an  indicator  of  the 
relative  supply/demand  variations  of  the  system  over  time. 


— 
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Since  the  reservoir  volume  output  data  is  a correct  measure 
of  the  actual  amount  of  water  stored  in  each  reservoir  sec- 
tion at  any  given  instant,  this  data  may  be  used,  if  neces- 
sary, to  calculate  the  true  energy  storage  variations  over 
time.)  Steady-state  TUS  nodal  temperature  and  mass  flowrate 
data  and  the  distribution  of  the  steady-state  building  energy 
demands  are  printed  out  in  the  steady-state  system  summary. 
Finally,  the  pressure  settings  of  all  the  loop  circulation 
pumps  are  initialized  to  the  steady-state  loop  fluid  fric- 
tional losses,  and  program  control  is  returned  to  Program 
MAIN  for  either  the  beginning  of  the  dynamic  simulation  or 
for  termination  of  the  run  if  only  the  steady-state  system 
analysis  output  is  desired. 

If  Subroutine  SET  is  called  during  the  dynamic  system 
analysis  (KODE  = 1),  its  function  is  to  transfer  the  desired 
nodal  fluid  mass  flowrates  computed  by  the  thermostat  model 
in  Program  MAIN  to  Subroutine  PRESS,  where  they  are  trans- 
lated into  a set  of  corresponding  loop  pressure  losses. 

These  pressures  are  then  shuttled  to  Subroutine  FLOWS,  in 
which  they  determine  the  desired  settings  of  the  loop  circu- 
lation pumps  to  compensate  for  the  out-of-balance  energy 
flows  detected  by  the  thermostat. 


A . 3 . 7 Fluid  Flowrate  (Subroutine  FLOWS) 

As  was  mentioned  previously,  the  primary  mechanism  em- 
ployed in  the  TUS  model  to  adjust  for  changes  in  the  system 
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loads  and  supply  conditions  is  the  variation  of  the  loop 
fluid  flowrates.  Subroutine  PLOWS  contains  the  model  for 
this  flowrate  adjustment  mechanism  and  also  serves  as  the 
calling  routine  for  the  dynamic  energy  balance  and  fluid 
pressure  loss  calculation  subroutines.  Since  the  consumer 
tertiary  distribution  piping  is  assumed  to  have  negligible 
fluid  and  thermal  inertia,  the  fluid  flows  in  the  consumer 
sides  of  the  load  center  heat  exchangers  are  allowed  to  vary 
instantaneously  with  the  thermal  demands.  However,  since 
the  secondary  distribution  and  primary  HTW  loops  contain  a 
considerable  volume  of  fluid,  circulation  pump  models  are 
provided  for  each  of  these  loops  which  determine  their 
steady-state  fluid  flowrates  and  the  rates  of  fluid  accelera- 
tion during  transients. 

At  the  end  of  .the  initiating  steady-state  calculations, 
the  fluid  friction  pressure  losses  corresponding  to  the 
converged  flows  are  calculated  for  each  of  the  loops  (see 
Subroutine  PRESS),  and  the  pump  pressures  are  set  to  main- 
tain these  flowrates.  When  the  thermostat  is  activated,  a 
set  of  desired  fluid  flowrates  is  calculated  in  Program  MAIN 
by  requiring  instantaneous  energy  balance  across  each  of  the 
heat  exchangers  in  the  system.  The  circulation  pump  pressures 
corresponding  to  these  desired  flows  are  calculated  in  Sub- 
routine PRESS  and  are  supplied  to  Subroutine  FLOWS,  where  they 
are  adjusted  to  speed  system  response  before  being  applied 
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to  the  actual  pumps.  If  the  thermal  loads  increase,  setting 
the  pump  pressures  to  the  values  necessary  merely  to  main- 
tain the  projected  loop  fluid  flowrates  will  cause  the 
fluids  to  accelerate  relatively  slowly  and  to  asymptotically 
approach  the  desired  conditions  over  a fairly  long  period 
of  time.  Since  it  is  desired  to  have  the  system  respond  as 
quickly  as  possible  to  changes  in  the  thermal  demands,  the 
fluid  acceleration  rates  are  increased  by  setting  the  actual 
pumps  at  pressures  somewhat  higher  than  the  required  steady- 
state  values.  The  factor  by  which  the  pressures  are  adjusted 
is  proportional  to  the  difference  between  the  existing  and 
the  required  pressures  and,  therefore,  varies  with  the  magni- 
tude of  the  change  in  the  thermal  demands.  By  adjusting  the 
proportionality  constant  in  this  correction  factor,  it  is 
possible  to  achieve  any  desired  degree  of  load  following  in 
the  system  fluid  flowrates.  However,  very  short  reaction 
times  (i.e.,  nearly  instantaneous  adjustment  of  the  flows 
to  exactly  match  the  load  conditions)  tend  to  cause  instabi- 
lities in  the  system  because  the  time  increments  over  which 
they,  .occur  must  be  of  finite  size  and  the  thermal  inertia 
of  the  system  does  not  allow  instantaneous  adjustment  of 
the  fluid  temperatures.  In  the  current  model,  the  pump 
pressure  adjustment  equation  used  is: 


PEX  = PP  + PSENS*(PP  - PS) 


f 


(A.  h9 ) 
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where  PP  is  the  desired  pressure  setting  computed  in  Sub- 
routine PRESS  from  the  desired  loop  flowrates,  PS  is  the 
existing  pump  pressure  setting,  and  a value  for  the  pump  • 
response  sensitivity  factor  PSENS  of  2.0  has  been  found  to 
be  capable  of  maintaining  realistic  pressure  variations  and 
system  stability  while  allowing  fairly  rapid  response  to 
demand  changes.  When  the  fluids  have  accelerated  to  within 
1%  of  their  projected  flowrates,  the  pump  pressures  (PEX) 
are  reduced  to  their  required  values  (PP)  to  prevent  exceed- 
ing the  projected  conditions  and  possible  instabilities  in 
the  resulting  heat  flows. 

If  the  thermostat  is  not  activated.  Subroutine  FLOWS 
controls  the  dynamic  energy  and  momentum  transfer  calcula- 
tions which  yield  the  instantaneous  fluid  temperatures  and 
flowrates  throughout  the  TUS.  In  these  dynamic  calculations, 
it  is  assumed  that  system  conditions  remain  constant  during 
each  time  interval,  and  new  values  for  the  fluid  temperatures, 
pressures  and  flowrates  are  calculated  at  the  end  of  a time 
step,  based  upon  the  instantaneous  system  thermal  energy 
supply  and  demand  conditions,  and  upon  the  overall  energy 
flows  and  fluid  momentum  transfers  during  the  preceding 
interval.  Since  this  is  a step-wise  approximation  to  the 
continuous  variation  of  these  parameters,  it  is  important  that 
the  specified  simulation  time  step  site  be  small  compared 
with  the  time  constants  of  the  externally-varying  heat  de- 
mand and  thermal  energy  supply  functions  to  maintain  system 
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stability . 

The  instantaneous  nodal  fluid  temperatures  are  deter- 
mined in  dynamic  energy  flow  calculations  in  Subroutine 
ERRS,  which  is  called  from  Subroutine  FLOWS  before  any  of 
the  nodal  fluid  flowrates  are  changed.  The  existing  fluid 
flows  in  the  secondary  distribution  and  primary  HTW  loops 
are  then  supplied  tc  Subroutine  PRESS,  where  the  total  fluid 
friction  pressure  losses,  including  the  effects  of  flow  regu- 
lation valves,  are  calculated  for  each  of  the  loops.  These 
pressure  losses  are  returned  to  Subroutine  FLOWS  and  are 
compared  with  the  circulation  pump  pressure  settings  to 
determine  the  acceleration  or  deceleration  of  the  fluid  dur- 
ing the  preceding  time  interval.  The  instantaneous  loop 
flows  at  the  end  of  the  time  step  are  cb rained  by  adjusting 
the  flowrates  at  the  beginning  of  the  interval  according  to 
the  calculated  acceleration  rate  through: 


DM  = MASS  + DTIME*FCTR*(PEX  - PP ) , 


(A. 50) 


'where  PP  is  used  here  as  the  instantaneous  fluid  friction 
pressure  loss  calculated  in  Subroutine  PRESS  for  the  beginning- 
of-interval  fluid  mass  flowrate  MASS,  and 


FCTR  = 


*gc(DIA)I 2 

(long) 


(A. 51) 


Q p 

with  gc  = i). 1^7x10°  lbm-f t/lbf-hr  . The  relative  flows  in 
the  parallel  branches  of  the  secondary  distribution  loops  are 
determined  by  the  new  total  loop  flowrates  given  by  Eq.  (A. 50) 
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and  the  settings  of  the  branch  flow  regulation  valves,  which 
divide  the  loop  flows  according  to  the  individual  load 
center  heat  exchanger  thermal  demands.  The  relative  bypass 
and  heat  exchanger  flowrates  at  each  secondary  supply  heat 
exchanger  in  the  primary  HTW  loop  are  similarly  determined 
according  to  the  thermal  demands  at  those  heat  exchangers. 

Subroutine  FLOWS  also  performs  the  energy  transfer  calcu- 
lations for  the  thermal  energy  storage  reservoir.  As  was 
mentioned  previously,  the  reservoir  is  modelled  physically 
as  two  separate  storage  volumes  in  series  with  the  primary 
loop  supply  and  return  nodes  at  the  TUS  supply  heat  exchanger 
(see  Fig.  A. 2).  The  initial  hot  water  storage  section  temper- 
ature is  specified  in  the  input  data  as  the  desired  HTW 
supply  temperature  for  the  primary  loop.  Since  large  varia- 
tions in  the  primary  loop  inlet  temperature,  coupled  with  the 
continuous  variation  of  the  fluid  flowrate  to  compensate  for 
load  changes,  tend  to  cause  instabilities  in  the  system  heat 
flows  during  periods  of  peak  or  minimum  demand,  this  tempera- 
ture is  held  constant  throughout  the  simulation.  The  cold 
water  storage  temperature  is  initialized  to  the  primary  loop 
steady-state  fluid  return  temperature.  The  fluid  flowrate 
between  the  two  reservoir  sections  is  then  calculated  by  re- 
quiring equality  between  the  convective  heat  transfer  rates 
on  both  sides  of  the  TUS  supply  heat  exchanger.  The  energy 
transferred  to  the  water  in  the  heat  exchanger  from  the 


power  plant  coolant  flow  is  equal  to  the  difference  between 
the  constant  power  plant  total  thermal  energy  output  and 
the  instantaneous  energy  demand  for  electrical  generation. 
Since  the  supply  heat  exchanger  outlet  temperature  is  con- 
stant and  the  inlet  temperature  varies  relatively  slowly  - 
as  determined  by  the  variations  in  the  primary  loop  return 
fluid  temperature,  buffered  by  the  large  volume  of  fluid  in 
the  cold  water  section  of  the  reservoir  - the  fluid  flowrate 
between  the  two  reservoir  sections  effectively  follows  the 
system  thermal  energy  supply.  At  the  end  of  each  simulation 
time  step,  the  water  volume  in  each  of  the  storage  sections 
is  adjusted  according  to  the  net  fluid  transfer  during  the 
interval.  (e.g.,  if  the  primary  loop  flowrate  is  greater 
than  that  in  the  supply  heat  exchanger,  there  will  be  a 
loss  of  fluid  from  the  hot  water  section  and  an  increase  in 
the  volume  of  the  cold  water  section. ) If  either  of  the 
sections  is  emptied  during  a given  time  step,  the  primary 
loop  fluid  flowrate  is  constrained  to  be  equal  to  the  flow- 
rate  through  the  thermal  energy  supply  heat  exchanger  to 
satisfy  conservation  of  mass  and  momentum  around  the  continu- 
ous piping  loop.  The  energy  contained  in  each  of  the  reser- 
voir sections  varies  with  both  the  volume  and  the  average 
temperature  of  the  water  being  stored.  At  the  end  of  each 


time  step,  this  energy  storage  is  adjusted  by  comparing  the 
thermal  energy  flows  into  and  out  of  the  section  during  the 
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interval.  (In  this  calculation,  it  is  assumed  that  heat 
conduction  losses  from  the  reservoir  will  be  small  compared 
with  the  convective  heat  flows  from  the  precooler  and  the 
primary  loop,  and  conduction  effects  are  ignored.  It  should 
also  be  noted  that,  in  the  present  version  of  the  code,  these 
energy  terms  are  only  relative  measures  of  the  actual  energy- 
stored  in  each  section,  since  they  are  computed  relative  to 
0 °F  rather  than  0 °K.  Since  the  volume  terms  depend  only 
upon  the  mass  flowrates  into  and  out  of  each  section,  they 
are  rigorously  correct  indicators  of  the  actual  reservoir 
behavior).  Thus,  during  periods  in  which  the  power  plant 
thermal  energy  supply  exceeds  the  system  demands,  the  volume 
of  water  in  the  hot  water  storage  section  increases.  During 
periods  of  higher  demand,  this  stored  energy  supplements  the 
power  plant  supply,  and  the  water  volume  decreases. 

A. 3- 8 System  Pressure  Drops  (Subroutine  PRESS) 

Subroutine  PRESS  calculates  the  total  fluid  friction 
pressure  losses  for  the  secondary  distribution  and  primary 
HTW  loops,  sets  the  circulation  pump  pressures  to  compensate 
for  these  friction  losses,  and  controls  the  flow  regulation 
valves  in  the  secondary  loop  parallel  supply  branches  to 
divide  the  loop  flows  according  to  the  relative  load  center 
heat  exchanger  thermal  demands.  The  subroutine  is  called 
during  the  initial  steady-state  system  calculations  to 
establish  the  initial  pump  pressures  and  valve  settings. 


when  the  thermostat  is  activated  to  determine  the  desired 
pressures  and  flow  splits  for  the  projected  system  flow- 
rates, and  at  the  end  of  each  simulation  time  step  to  calcu- 
late the  instantaneous  fluid  pressure  losses  in  each  loop. 
The  basis  of  the  nodal  fluid  friction  pressure  loss  calcula- 
tion is  the  Darcy  formula  [12]: 


DP  = F* ( LONG/DI A ) * ( 


MASS* 


2*G*RH0*MTB2*AC2  ' 


(A. 52) 


where  the  familiar  form  of  the  equation  has  been  modified 
slightly  by  expressing  the  fluid  velocity  in  terms  of  the 
mass  flowrate  (MASS),  the  flowchannel  cross-section  (AC), 
and  the  number  of  parallel  flowohannels  (NTB).  The  factor  F 
is  the  Darcy-V/eisbach  friction  factor 


F 


0.184 


RE 


0.2 


(A. 53) 


RE  is  the  fluid  Reynolds  number 
. ..  RE  = tt*NT3*DIA*MU  » 


(A. 54) 


with  the  fluid  velocity  again  expressed  in  terms  of  the  mass 

g 

flowrate,  and  the  conversion  constant  G = 4.147x10  lbm-ft/ 
lbm-hr2 . 

In  a simple  series  loop,  such  as  the  primary  HTW  loop, 
the  total  frictional  pressure  loss  is  simply  the  sum  of  the 


nodal  losses  determined  through  Eq.  (A. 52).  Since  the  circu- 
lation pump  for  the  loop  must  compensate  exactly  for  these 
losses  in  steady-state,  the  pump  pressure  necessary  to  main- 
tain the  given  loop  fluid  flowrate  is  set  equal  to  the  cal- 
culated total  loop  frictional  pressure  loss. 

In  calculating  the  friction  losses  for  a secondary  dis- 
tribution loop,  consideration  must  be  given  to  the  fact  that 
the  total  fluid  flow  is  split  among  several  parallel  branches 
serving  the  load  center  heat  exchangers.  Since  all  the 
branches  in  a given  loop  have  common  supply  and  return  head- 
ers, the  fluid  flowing  through  each  branch  must  experience 
the  same  total  pressure  drop.  The  flow  regulation  valve 
model  employed  in  Subroutine  PRESS  uses  this  fact  to  intro- 
duce external  pressure  losses  into  selected  branches  in  order 
to  maintain  different  flowrates  in  similar  flowchanne]  geomet- 
ries. The  desired  nodal  fluid  flowrates  are  supplied  to  Sub- 
routine PRESS  from  either  the  steady-state  system  calculation 
(Subroutine  STDY)  or  the  thermostat-controlled  system  projec- 
tion (Program  MAIN).  Since  all  the  calculations  which  estab- 
lish these  flowrates  are  based  upon  instantaneous  energy 
balances  across  each  of  the  heat  exchangers  in  the  system,  the 
flows  are  divided  among  the  parallel  branches  according  to 
the  relative  thermal  loads  which  they  must  supply.  In  each 
secondary  distribution  loop,  the  nodal  friction  pressure  loss- 
es are  calculated  through  Eq . (A. 52),  and  the  largest 
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parallel  branch  pressure  drop  is  taken  as  the  reference  for 
the  remainder  of  the  branches.  It  is  assumed  that  the  valve 
in  this  reference  branch  is  fully  open,  and  the  valves  in 
the  remaining  branches  are  set  by  adding  a valve  friction 

t 

term  to  each  nodal  friction  factor  to  make  all  the  parallel 
branch  pressure  losses  equal  to  the  reference.  This  paral- 
lel branch  pressure  drop  is  then  added  to  the  series  node 
frictional  losses  to  determine  the  total  fluid  pressure  drop 
for  the  loop  and  to  establish  the  circulation  pump  pressure 
necessary  to  maintain  the  desired  flowrates. 

During  normal  simulation  periods,  the  instantaneous 
flowrates  determined  in  Subroutine  FLOWS  are  applied  to  Sub- 
routine PRESS,  where  the  loop  fluid  friction  pressure  losses 
corresponding  to  these  flows  are  calculated.  These  frictional 
losses  are  then  compared  with  the  loop  circulation  pump 
pressure  settings  to  determine  the  rates  of  fluid  accelera- 
tion or  deceleration  as  indicated  in  the  section  describing 
Subroutine  FLOWS.  The  parallel  branch  flow  regulation  valves 
are  set  at  the  beginning  of  the  simulation  and  every  time  the 
thermostat  is  activated  to  adjust  the  pump  pressures.  During 
normal  dynamic  calculations,  however,  their  settings  (i.e., 
the  added  terms  in  the  parallel  branch  friction  factors) 
remain  constant  as  determined  by  the  desired  flow  conditions, 
while  the  fluid  friction  losses  for  each  of  the  nodes  change 
as  the  loop  flowrate  varies.  The  valves  thus  establish  the 
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relative  parallel  branch  fluid  flow  splits  based  upon  the 
desired  system  energy  supply  conditions  and  maintain  these 
flow  ratios  until  they  are  reset  by  the  thermostat.  (No 
models  of  the  bypass  valves  in  the  primary  loop  are  contained 
in  Subroutine  PRESS,  since  those  valves  are  set  through  the 
combined  effects  of  the  thermostat  flow  split  calculations  in 
Program  MAIN  and  the  pump  pressure  resetting  section  of  Sub- 
routine FLOWS.) 


A. 3-9  Dynamic  Energy  Balance  (Subroutine  ENRG) 

The  dynamic  energy  flow  calculations  in  Subroutine  ENRG 
form  the  basis  for  the  entire  time-dependent  TUS  simulation 
since  they  determine  the  rate  at  which  the  distribution  sys- 
tem follows  the  thermal  energy  supply  and  demand  variations. 
The  system  fluid  flowrates  and  nodal  temperatures  are  assumed 
to  be  constant  during  each  time  interval  according  to  the 
net  imbalances  in  the  heat  flowrates  and  momentum  transfers 
computed  during  the  time  step.  Subroutine  ENRG  calculates 
the  nodal  inlet  and  outlet  fluid  temperatures  throughout 
the  TUS  based  upon  an  analysis  of  the  nodal  thermal  energy 
flows  during  each  time  step.  It  is  called  from  Subroutine 
FLOWS  before  any  other  system  parameters  are  adjusted  because, 
unless  the  thermostat  is  activated,  the  fluid  flowrates  vary 
independently  of  the  fluid  temperatures,  while  the  energy 
flows  and  temperature  variations  depend  strongly  upon  the 
fluid  mass  flowrates  existing  during  the  interval  being 
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s cudied . 

All  energy  transfer  calculations  in  Subroutine  ENRG  have 
seen  derived  under  the  assumption  that  the  fluid  in  a given 
node  travels  only  a relatively  small  fraction  (e.g.,  10-20*) 
of  the  length  of  the  node  during  a single  time  step.  (This 
criterion  should  be  the  dominant  factor  in  the  determination 
of  the  time  step  size  specified  for  a given  simulation.) 
During  transients,  the  fluid  inlet  and  outlet  temperatures 
on  either  side  of  a heat  exchanger  may  be  perturbed  so  as  tc 
cause  the  temperature  profiles  to  deviate  from  their  steady- 
state  shapes.  In  order  to  follow  such  perturbations  along 
the  entire  length  of  a node,  it  is  rigorously  necessary  to 
incorporate  into  the  analysis  r.n  incremental  flow  indexing 
scheme  which  traces  the  time  history  of  each  element  of 
fluid  as  it  travels  through  each  node  in  its  respective 
piping  loop.  However,  the  effects  of  these  continuous  varia- 
tions in  the  fluid  temperature  profiles  may  be  modelled  ade- 
quately for  most  situations  encountered  under  normal  system 
operating  conditions  by  concentrating  only  upon  the  incre- 
mental energy  flow  variations  at  the  inlet  and  outlet  of 
each  node.  Subroutine  ENRG  therefore  examines  only  the 
• ermal  energy  transferred  into  and  out  of  the  incremental 
. jme  of  fluid  displaced  during  the  time  interval  being 
• d In  Pig.  A. 3,  it  is  assumed  that  a particle  of 
■ • x = L - f.  at  the  beginning  of  a time  step  (t  = 0) 


Dynamic  Nodal  Energy  Transfer  Parameters 


TODP  TO 


NTB  channels/ node , 
single  channel  cross-section  AC 
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flows  to  x = L at  the  end  of  the  interval  (t  = DTIME);  thus 

l = v* DTIME  , (A. 55) 

where  v is  the  nodal  fluid  velocity.  Expressing  the  fluid 
velocity  in  terms  of  the  nodal  mass  flowrate,  the  fraction 
of  the  length  of  the  node  travelled  by  the  fluid  in  DTIME 
is 


£ = nr  - MASS*DTIME  , „ 

L " NTB*FtHO*AC*LONG  ' 

However,  since  (NTB*RHO*AC*LONG ) is  just  the  total  fluid  mass 
contained  within  the  node  (TMAS),  Eq . (A. 56)  can  be  rewritten 
as 


DL 


Mass*  DTIME 
TMAS 


(A. 57) 


During  the  interval  DTIME,  the  fluid  in  the  n 1e  trans- 
fers heat  to  (or  receives  heat  from)  its  surroundings.  There- 
fore, an  element  of  fluid  located  at  position  x = L-  £ at  the 

beginning  of  a time  step  will  experience  a temperature  change 
* ✓ ► 

during  the  interval  due  to  the  net  heat  loss  or  gain  it 
experiences  while  traversing  the  distance  £.  Denoting  the 
temperature  of  the  fluid  at  t = 0 , x = L-£,  as  TODP,  the 


temperature  of  the  fluid  at  t = DTIME,  x = L,  as  TOP,  and  the 
temperature  change  of  the  fluid  as  it  passes  from  x = L-£  to 
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nodal  fluid  outlet  temperature  at  t = TDIME  is 


TOP  = TODP  - DTC  . (A. 58) 

At  any  given  instant,  the  only  parameters  specified  for 
each  node  are  the  fluid  mass  flowrate,  the  Inlet  and  oulet 
temperatures,  and  the  external  conditions  affecting  the 
heat  transfer  rate  from  the  node  (e.g.,  the  ambient  earth 
temperature  for  a pipe  node  or  the  fluid  mass  flowrate  and 
inlet  and  outlet  temperatures  on  the  opposite  side  of  a 
heat  exchanger) . The  detailed  nodal  temperature  profiles  are 
not  calculated,  necessitating  an  approximation  to  be  made  in 
the  determination  of  the  fluid  temperature  TODP.  An  exponen- 
tial temperature  difference  profile  is  assumed  for  each  node, 
having  the  general  form: 


where 


Td(x)  = TN  exp (-ax ) , (A. 59) 

T^(x)  = temperature  difference  between  the  fluid 
in  the  node  and  its  surroundings  at  posi- 
tion x, 

TN  = temperature  difference  between  the  fluid  in 
the  node  and  its  surroundings  at  the  node 
inlet,  and 

a = temperature  difference  decay  constant. 


t 


Since  the  nodal  outlet  temperature  difference,  TD,  is  also 


know 


(A. 59)  can  be  solved  to  yield  a value  for  a in 


r* , tj  q . 

terns  cf  the  instantaneous  values  of  TN  and  TD: 

~aL  = ln  > (A. 60) 

and  the  temperature  difference  profile  becomes* 

Td(x)  = TN  exp[ln(g)*£]  . (A.6l) 

It  Is  also  assumed  that  the  slope  of  the  fluid  temperature 
profile  at  any  point  along  the  length  of  a node  is  directly 
proportional  to  the  heat  transfer  rate  from  the  node  at  that 
point  and,  hence,  to  the  difference  between  the  nodal  fluid 

*It  should  be  noted  that  this  exponential  difference  profile 
will  yield  the  well-known  log-mean  temperature  difference  for 
steady-state  heat  conduction  in  a counterflow  heat  exchanger: 


fL 

TrdU  T , (x  )dx 
■*  o 


= AU  AT 


£mc 


This  profile  is  thus  valid  for  steady-state  heat  transfer  cal- 
culations and  can  be  extended  to  transient  calculations  of  the 
type  made  in  the  current  model  with  little  error  in  the 
resulting  heat  flows. 
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and  external  surroundings  temperatures  at  that  point.  If 
the  average  slope  of  the  fluid  temperature  profile  for  the 
entire  node 

TS(av)  = TI Z 'T'Q  (A. 62) 


is  the  result  of  heat  conduction  from  the  node  at  the  aver- 
age log-mean  temperature  difference  between  the  fluid  and  its 
surroundings 


DTLM 


TD  - TN 

/TDn 
In  ( 


(A. 63) 


then  the  ratio  of  the  average  slope  of  the  temperature  profile 
in  the  region  between  x = L - i and  x = L to  the  average  slope 
for  the  entire  node  is  given  by  the  ratio  of  the  log-mean 
temperature  difference  over  that  region  to  the-  average  nodal 
difference.  (i.e.,  The  ratio  of  the  profile  slopes  is  equal 
to  the  ratio  of  the  heat  conduction  rates.)  Since  the  log- 
mean  temperature  difference  over  the  region  of  the  node  from 
x = L - l to  x = L is  given  by: 


DTLMX 


TD  - TDX 


r TD  . 
in  (TDX) 


(A. 64) 


with  TDX  being  calculated  through  Eq.  (A. 61)  at  y = L-£,  the 
average  slope  of  the  nodal  fluid  temperature  profile  over 
this  region  of  the  node  can  be  approximated  by 


► i 


T (x)  = T (av'1 
1SU)  XS^a  ' DTLM 


(A. 65) 


Knowing  the  average  slope  of  the  temperature  profile  and  the 
nodal  fluid  outlet  temperature,  it  is  possible  to  estimate 
the  temperature  of  the  fluid  at  x = L - l to  be: 


TODP  = TO  + DL*  (TI  - TO)  * ( ^lTT ) * 


( A . 66 ) 


In  order  to  calculate  the  temperature  change  of  the 
fluid  due  to  heat  conduction  effects  as  it  flows  from  x = L - l 
to  x = L,  it  is  necessary  to  determine  the  average  heat  trans- 
fer rate  from  the  fluid  to  its  surroundings  over  this  region. 
Integrating  the  assumed  exponential  temperature  difference 
profile  from  x = L-H  to  x = L yields  the  result  that  the 
average  conduction  heat  transfer  rate  from  the  fluid  is 


• _ TD  - TDX 

Qc  ~ AUF  777TD  — 
In  %) 


(A. 67) 


where  AUF  is  the  overall  heat  transfer  coefficient  for  the 
tdtal'node.  If  it  is  assumed  that  the  heat  lost  from  the 
incremental  volume  of  fluid  in  node  length  i during  DTIME 
results  in  a change  in  the  average  fluid  temperature  of  mag- 
nitude DTC,  energy  conservation  requries 


QC*DTIME  = NTB*RHO*AC*£*CAP*DTC  . 


(A. 68) 
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Using  Eq.  (A. 56)  to  express  l in  terms  of  the  fluid  mass 
flowrate  and  simplifying  the  resulting  expression  yields  the 
average  temperature  change  of  the  fluid  due  to  conduction 
heat  losses : 


DTC 


MASS* CAP 


(A. 69) 


The  expressions  for  the  nodal  fluid  outlet  temperature 
assuming  no  heat  loss  (Eq.  (A. 66))  and  for  the  temperature 
change  due  to  heat  conduction  effects  (Eq.  (A. 69))  are  com- 
bined through  Eq.  (A. 58)  to  yield  a new  value  for  the  nodal 
fluid  outlet  temperature  at  the  end  of  the  time  step  being 
studied.  As  indicated  by  the  preceding  discussion,  the  cal- 
culations leading  to  this  temperature  essentially  consist  of 
a comparison  between  the  convective  and  conductive  heat  flow- 
rates  from  the  fluid  in  the  final  volume  increment  of  the 
node.  If  the  nodal  thermal  energy  demand  increases  over  its 
steady-state  value,  DTC  will  increase  in  magnitude  and  the 
outlet  temperature  will  decrease  as  heat  is  removed  from  the 
fluid  at  a rate  faster  than  it  is  supplied  to  the  node  by 
the  steady-state  fluid  flowrate.  If  the  demand  decreases, 
the  outlet  temperature  will  increase  as  the  fluid  is  swept 
through  the  node  faster  than  heat  is  removed  from  it  through 
conduction. 

The  basic  nodal  analysis  outlined  above  for  heat  exchan- 


88. 


ger  nodes  is  also  employed  in  the  adjustment  of  pipe  node 
outlet  temperatures,  with  the  only  modifications  being  that 
the  temperature  "profile"  of  the  "opposite  node"  is  the  con- 
stant earth  temperature  TAM3  and  that  - due  to  the  generally 
small  differences  between  the  pipe  node  inlet  and  outlet 
temperatures  - all  exponential  temperature  relations  are 
approximated  by  linear  functions. 

Following  computation  of  all  the  revised  nodal  outlet 
temperatures,  the  nodal  inlet  temperatures  are  set  by  requir- 
ing continuity  of  the  fluid  temperature  profiles  across  the 
boundaries  between  adjacent  nodes.  For  series  nodes,  the 
inlet  temperature  of  a node  is  simply  set  equal  to  the  new 
outlet  temperature  of  the  preceding  node  in  the  loop.  The 
inlet  temperatures  to  the  pipe  nodes  following  the  primary 
loop  secondary  supply  heat  exchangers  are  given  by  the  heat 
exchanger  and  bypass  fluid  mass-flowrate-average  of  the  pre- 
ceding node’s  inlet  and  outlet  temperatures.  The  inlet  temp- 
erature to  the  return  pipe  in  a secondary  distribution  loop 
is  the  mass-f lowrate-averagea  outlet  temperature  from  the 
parallel  load  center  heat  exchanger  branches.  The  inlet 
temperatures  to  the  consumer-side  nodes  of  the  load  center 
heat  exchangers  are  calculated  in  Program  MAIN  from  the 
heat  exchanger  fluid  flowrates,  the  fluid  outlet  temperature 
computed  through  the  dynamic  heat  exchanger  analysis  in  this 
subroutine,  and  the  thermal  demands  obtained  from  linear 
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interpolation  of  the  hourly  loads  computed  in  Subroutine 
DLOAD . 

During  the  computation  of  the  new  nodal  temperatures. 
Subroutine  ENRG  institutes  a fairly  complicated  system  of 
temperature  monitoring  and  adjustment  calculations  to  insure 
that  the  approximations  made  do  not  violate  physical  heat 
transfer  laws  or  lead  to  mathematical  relations  exceeding 
the  computational  capabilities  of  the  computer.  The  heat 
exchanger  inlet  and  outlet  temperature  differences,  TN  and 
TD,  are  allowed  to  have  minimum  magnitudes  of  10-^  °F  to 
insure  the  finiteness  of  calculations  involving  logarithms 
of  their  ratio.  Similarly,  if  the  ratio  TN/'TD  is  within 
10~6  of  1.0  (i.e.,  if  the  fluid  temperature  profiles  are 
virtually  parallel  along  the  entire  length  of  the  heat  exchan- 
ger), the  energy  flow  calculations  described  above  ar°  per- 
formed using  linear  approximations  to  all  the  exponential 
relations,  thereby  avoiding  the  calculation  of  logarithms  of 
0,  dividing  by  0,  etc.  Although  TN  and  TD  are  allowed  to 
have  opposite  signs  (indicating  a crossing  of  the  temperature 
profiles  due  to  instantaneous  transient  conditions  at  either 
end  of  the  heat  exchanger),  .if  this  condition  is  detected, 
linear  approximations  to  the  exponential  relations  are 
employed  to  return  the  parameters  to  the  unidirectional  ener- 
gy flow  conditions  satisfying  the  Second  Law  of  Thermodynamics. 
Finally,  if  the  computed  nodal  outlet  temperatures  are 
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found  to  violate  the  Second  Law  of  Thermodynamics  (o.g.,  if 
the'  outlet  fluid  temperature  on  the  supply  side  of  a heat 
exchanger  is  lower  than  .the  inlet  fluid  temperature  on  the 
demand  side),  finitely  small  temperature  differences  of 
magnitude  10-^  °F  and  of  the  correct  sign  are  assigned 
between  the  errant  values,  and  the  system  is  allowed  to 
correct  itself  through  the  continuing  transient  temperature 
adjustments  made  during  the  succeeding  time  intervals. 

A. 3-10  Precooler  Monitor  (Subroutine  PREC) 

The  precise  identification  assigned  to  the  TUS  supply 
heat  exchanger  will  depend  upon  the  method  and  location  in 
which  heat  is  extracted  from  the  power  plant's  coolant  flow- 
stream.  However,  for  many  plant  types,  this  heat  exchanger 
could  take  the  form  of  a coolant  precooler  inserted  into  the 
turbine  exhaust  flowstream  upstream  of  the  fluid  compressors 
[13] • Although  this  heat  exchanger  is  not  modelled  by  TDIST 
in  the  detail  afforded  to  the  TUS  distribution  loop  units, 
its  temperature  and  fluid  flowrate  parameters  are  monitored 
during  each  simulation  time  step  to  insure  that  the  power 
plant-TUS  energy  interchange  is  at  least  obeying  the  physical 
constraints  of  the  Second  Law  of  Thermodynamics.  The  heat 
transfer  coefficient  of  this  heat  exchanger  is  assumed  to  be 
"very  large",  such  that  any  variations  in  the  fluid  tempera- 
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tures  on  one  side  are  immediately  reflected  in  a temperature 
variation  on  the  opposite  side.  In  the  current  version  of 


TDIST,  Subroutine  PREC  simply  adjusts  the  power  plant  coolant 
flowstream  inlet  and  outlet  temperatures  and  the  coolant 
flowrate  to  insure  that  the  coolant  temperatures  are  always 
greater  than  (or  in  the  limiting  case,  equal  to)  the  water 
temperatures  on  the  TUS  side  of  the  heat  exchanger  and  that 
the  rate  of  heat  transfer  matches  the  power  plant ' s net  rate 
of  thermal  energy  production. 

A. 3. 11  Output  (Subroutine  PRV) 

Subroutine  PRV  provides  the  primary  output  function  for 
the  code  and  is  called  at  regular  intervals  specified  by  the 
user-defined  print  frequency  code.  Output  data  includes  the 
nodal  fluid  mass  flowrates  and  inlet  and  outlet  temperatures, 
the  instantaneous  heat  flowrates  through  each  heat  exchanger, 
information  on  the  volume  of  water  stored  in  both  sections 
of  the  reservoir,  the  temperatures  and  flowrates  on  both 
sides  of  the  TUS  supply  heat  exchanger,  and  the  instantaneous 
loop  circulation  pump  pressure  settings. 

A. 4 Simulation  Run  Options 

Figure  A, 9 is  a simplified  flowchart  for  TDIST  illustrat- 
ing the  major  program  components  and  decision  control  branches. 
Depending  upon  the  type  and  extent  of  the  system  analysis 
desired,  the  code  provides  the  flexibility  of  three  general 
run  options:  (1)  the  calculation  and  distribution  of  the 

individual  building  thermal  and  electrical  energy  demands, 
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(2)  the  iterative  steady-state  convergence  of  the  TUS  nodal 
temperatures  and  fluid  flowrates,  and  (3)  the  full  time- 
dependent  TUS  energy  flow  analysis.  The  following  sub-sec- 
tions describe  briefly  the  computation  format  and  parameter 
specifications  for  each  of  these  options.  Detailed  input 
data  requirements  and  output  listings  for  the  code  are- des- 
cribed in  sections  A. 5 and  A. 6. 

Two  user-defined  program  flow  control  codes  determine 
the  simulation  option(s)  chosen  for  a given  computer  run  and 
set  the  paths  for  the  inter-subroutine  data  transfer  and  in- 
formation outputs.  The  primary  run  control  code  is  KODE, 
which  determines  the  type  of  analysis  to  be  performed. 

KODE  = 0 allows  either  the  calculation  and  output  of  the 
hourly  building  energy  demands  and  the  distribution  of  the 
initial  loads  according  to  the  TUS  load  center  heat  exchanger 
locations  or,  depending  upon  the  value  of  the  second  control 
code  (IDES),  it  causes  the  TUS  steady-state  analysis  to  be 
performed.  KODE  = 1 indicates  that  the  hourly  building  use 
and  weather  data  specifications  should  be  read  in  and, 
unless  it  is  over-ridden  by  IDES,  causes  the  dynamic  system 
simulation  to  be  performed.  KODE  = 2 indicates  that  the 
run  is  to  be  terminated.  Whenever  the  value  KODE  = 0 is 
input  to  the  program,  the  design  run  option  code  IDES  must 
be  specified.  IDES  = 0 indicates  that  the  run  to  be  made 
is  simply  the  TUS  steady-state  analysis,  based  upon  a computed 
single  time  instant  building  load  distribution.  IDES  = 1 
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allows  the  calculation  of  the  building  loads  and  their  dis- 
tribution among  the  specified  load  center  heat  exchangers, 
but  terminates  the  run  prior  to  the  actual  application  of  the 
demands  to  the  heat  exchangers.  The  hourly  building  loads 
for  the  entire  simulation  period  are  computed  and  printed 
out  if  the  control  code  sequence  KODE  = 0,  IDES  = 1,  KODE  = 1 
is  input.  Table  A. 9 summarizes  these  basic  control  functions. 

A. 4.1  Building  Energy.  Demands  (KODE  = 0,  IDES  = 1) 

This  run  option  allows  the  calculation  of  the  magnitudes 
and  distribution  of  the  building  energy  demands  without  apply- 
ing the  loads  to  the  TUS  load  center  heat  exchangers  or  per- 
forming any  TUS-related  analyses.  When  the  sequence  KODE  = 0, 
IDES  - 1 is  input  to  the  program,  single  time  instant  values 
are  read  in  for  each  building  type's  desired  room  temperature, 
equipment  use  factor  and  domestic  hot  water  use  factor.  Also 
input  are  the  ambient  air  temperature,  wind  velocity,  wind 
direction  and  non-space-conditioning  electrical  demand  factor 
for  the  community  at  the  time  instant.  Subroutine  QLOAD  is 
called,  and  the  computed  building  space  conditioning  loads 
are  printed  out.  These  loads  are  supplied  to  Subroutine 
DLOAD,  in  which  they  are  scaled  by  the  appropriate  end-use 
equipment  coefficients  of  performance  and  are  aggregated  into 
load  center  thermal  energy  demands  and  a total  electrical 
demand  for  the  community.  The  distribution  of  the  building 
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TABLE  A. 9 

TDIST  RUM  OPTION  CONTROL  CODE  SPECIFICATIONS 


Code 

Value 

Resulting  Program  Performance 

KODE 

IDES 

0 

0 

Building  demands  calculated  at  one  time 
instant  and  distributed  to  TUS : TUS 
steady-state  analysis  performed  and 
system  summary  printed 

KODE 

IDES 

0 

1 

Building  demands  calculated  at  one  time 
instant  and  distributed  to  TUS:  build- 
ing loads  and  distribution  printed 

KODE 

IDES 

KODE 

0 

1 

1 

( sequence ) 

Building  demands  calculated  at  one  time 
instant  and  distributed  to  TUS:  build- 
ing loads  and  distribution  printed; 
hourly  buildings  loads  calculated  and 
printed 

KODE 

1 

Dynamic  TUS  simulation  performed;  TUS 
parameter  summary  and  load  distribution 
printed  at  user-specified  output 
intervals 

KODE 

2 

Run  terminated 

i 
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loads  is  printed  out  by  building  type  at  each  load  center 
heat  exchanger,  and  a summary  is  provided  of  the  aggregate 
thermal  energy  demands  at  each  load  center  heat  exchanger. 
When  control  is  returned  to  Program  MAIN,  a succeeding  value 
of  KODE  must  be  input.  This  value  of  KODE  should  not,  in 
general,  be  equal  to  0,  unless  it  is  desired  to  repeat  the 
same  calculation  sequence  or  to  perform  a TUS  steady-state 
analysis.  If  KODE  = 2,  the  run  is  terminated.  If  KODE  = 1, 
up  to  24  hourly  data  values  for  the  building  temperatures 
and  use  factors,  the  hourly  weather  data,  and  the  hourly  non- 
space-conditioning demand  factors  are  read  in,  are  printed 
out  in  a set  of  Input  data  summary  tables,  and  are  supplied 
to  Subroutine  QLOAD.  An  hour-by-hour  space  conditioning  load 
calculation  is  performed  for  each  building  type,  and  the 
resulting  hourly  load  data  is  printed  out.  Following  this 
printout,  the  run  is  terminated. 

A. 4. 2 TUS  Steady-State  Analysis  (KODE  = 0,  IDES  = 0) 

This  run  option  initializes  the  TUS  to  steady-state 

energy  and  momentum  flow  conditions  based  upon  a set  of  com- 
* ✓ - 

puted  building  energy  demands  and  a specified  energy  supply 
rate  at  the  TUS  supply  heat'  exchanger.  When  the  sequence 
KODE  = 0,  IDES  = 0 is  input,  single  time  instant  values  are 
read  in  for  each  building  type's  desired  room  temperature, 
equipment  use  factor  and  domestic  hot  water  use  factor.  Also 


input  are  the  ambient  air  temperature,  wind  velocity,  wind 
direction  and  non-space-conditioning  electrical  demand  factor 
for  the  community  at  that  time  instant.  Subroutine  QLOAD  is 
called,  and  the  computed  building  space  conditioning  loads 
are  printed  out.  These  loads  are  supplied  to  Subroutine 
DLOAD,  in  which  they  are  scaled  by  the  appropriate  end-use 
equipment  coefficients  of  performance  and  are  aggregated 
into  load  center  thermal  energy  demands  and  a total  electrical 
demand  for  the  community.  Prom  these  loads  and  the  user-speci- 
fied supply  and  return  water  temperatures  at  each  load  center 
heat  exchanger  are  calculated  the  steady-state  fluid  mass 
flowrates  in  all  the  consumer  energy  distribution  piping 
nodes.  The  total  space  conditioning,  domestic  hot  water, 
and  non-space-conditioning  electrical  loads  are  subtracted 
from  the  specified  power  plant  total  thermal,  power  output  to 
obtain  the  net  rate  of  heat  supply  to  the  TUS.  Subroutine 
STDY  is  called,  in  which  iterative  steady-state  heat  balance 
calculations  are  employed  to  determine  a converged  set  of 
nodal  fluid  mass  flowrates  and  temperatures  which  satisfy 
the  combined  energy  supply  and  demand  conditions  faced  by 
the  TUS.  The  resulting  loop  fluid  flowrates  are  transmitted 
to  Subroutine  PRESS,  and  the  corresponding  steady-state  loop 
circulation  pump  pressure  settings  are  determined.  Although 
the  primary  loop  heat  exchangers  are  analyzed  in  steady-state 
with  no  bypass  flows  around  them,  the  parallel  branch  flow 
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regulation  valves  in  the  secondary  loops  are  set  according 
to  the  relative  thermal  energy  demands  at  each  load  center 
heat  exchanger.  The  fluid  temperatures  and  flowrates  on 
both  sides  of  the  power  plant's  TUS  supply  heat  exchanger 
are  determined,  and  the  temperatures  of  the  water  contained 
in  the  hot  and  cold  sections  of  the  storage  reservoir  are 
initialized  to  the  steady-state  supply  and  return  tempera- 
tures in  the  primary  HTW  loop.  The  steady-state  TUS  para- 
meters are  printed  out  in  a nodal  energy  flow  summary, 
followed  by  a printout  of  the  distribution  of  the  building 
space  conditioning  energy  demands  by  building  type  at  each 
load  center  heat  exchanger.  When  control  is  returned  to 
Program  MAIN,  a succeeding  value  of  KODE  must  be  specified. 

If  KODE  = 0,  the  user  has  the  option  of  performing  a second 
steady-state  system  analysis  or  of  calling  the  building  load 
analysis  run  option  described  in  Section  A. 4.1.  If  KODE  = 1, 
the  computed  steady-state  nodal  parameters  are  used  to  ini- 
tialize the  TUS  for  the  dynamic  energy  flow  simulation  des- 
cribed in  Section  A. 4. 3.  If  KODE  = 2,  the  run  is  terminated. 

^ »•  » | 

A. 4. 3 Dynamic  System  Analysis  (KODE  = 1) 

Unless  it  is  suppressed  by  the  condition  IDES  = 1 (see 
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section  A. 4.1),  an  input  value  of  KODE  = 1 causes  the  dynamic 
energy  system  analysis  to  be  performed  for  a similation  period 
of  up  to  24  hours  in  length.  Although  it  is  theoretically 


possible  to  specify  a consistent  set  of  TUS  input  data 
which  will  allow  the  dynamic  system  analysis  to  be  success- 
fully started  the  first  time  the  program  requires  KODE  to 
be  input,  the  stability  of  the  time-dependent  simulation  is 
so  dependent  upon  the  system  initialization  parameters  that 
it  is  recommended  that  the  steady-state  TUS  initialization 
run  option  (KODE  = 0,  IDES  = 0)  always  be  used  as  the  front 
end  of  the  dynamic  system  analysis.  In  addition  to  reducing 
the  amount  of  hand  calculations  required  of  the  user  and 
minimizing  the  chances  for  TUS  instabilities  developing  due 
to  poorly  specified  initial  parameter  values,  use  of  the 
internal  initialization  option  allows  the  user  to  run  low 
cost  design  calculations  to  insure  system  convergence  prior 
to  gambling  the  larger  sums  of  time  and  money  associated 
with  the  full  dynamic  analysis  on  a set  of  untested  input 
data.  Therefore,  in  the  following,  it  is  assumed  that  KODE = 1 
is  read  into  the  computer  upon  return  of  control  to  Program 
MAIN  from  the  steady-state  analysis  calculations  described 
in  Section  A. 4. 2. 

Input  data  for  each  building  type  consists  of  up  to  24 
hourly  desired  room  temperatures,  domestic  hot  water  usage 
factors,  and  occupancy-related  building  use  factors  which 
are  used  to  schedule  internal  equipment  heat  generation  and 
forced  air  ventilation.  Weather  data  inputs  include  hourly 
ambient  air  temperatures,  wind  velocities  and  wind  directions. 


An  hourly  community  non-space-conditioning  electrical  energy 
demand  schedule  is  also  required  for  the  determination  of 
the  components  of  the  power  plant  load  not  directly  associated 
with  providing  heating,  cooling,  or  hot  water  to  each  person 
in  the  community.  These  hourly  schedules  are  printed  out  in 
input  data  summary  tables  and  are  supplied  to  Subroutine 
QLOAD.  Hourly  net  space  heating  and  cooling  energy  demands 
are  determined  for  each  specified  building  type,  and  the 
building  load  schedules  are  printed  out.  The  individual 
building  demands  are  also  supplied  to  Subroutine  DLOAD, 
which  computes  and  stores  the  resulting  set  of  hourly  thermal 
demands  at  each  load  center  heat  exchanger  and  the  hourly 
total  electrical  energy  demand  for  the  community.  Initiali- 
zation calculations  in  Subroutine  SET  transfer  the  steady- 
state  nodal  parameters  calculated  in  Subroutines  STDY  and 
PRESS  to  the  appropriate  nodes  in  the  TUS  model,  and  the 
primary  loop  and  load  center  heat  exchanger  consumer-side 
thermostat  reference  temperature  limits  are  set  in  Program 
MAIN  based  upon  the  converged  steady-state  fluid  temperatures 
and  the  design  values  specified  by  the  user.  An  internal 
time  clock  is  set,  and  the  dynamic  analysis  begins  at  the 
hour  of  the  day  at  which  the  steady-state  system  parameters 
are  specified. 

During  each  simulation  time  step,  linear  interpolations 
between  the  stored  hourly  data  values  for  the  air  temperature , 
system  electrical  loads,  and  load  center  heat  exchanger  ther- 


mal  demands  provide  the  instantaneous  energy  demand  condi- 
tions applied  to  the  TUS  and  the  power  plant.  The  rate  of 
thermal  energy  supply  to  the  TUS  is  determined  by  subtracting 
the  instantaneous  total  electrical  energy  demand  from  the 
specified  power  plant  constant  total  thermal  power  output. 
Subroutine  ENRG  is  called  to  determine  the  instantaneous 
fluid  temperature  distribution  throughout  the  TUS,  based 
upon  the  assumption  of  constant  energy  flow  conditions 
during  the  time  step,  with  temperature  adjustments  made  at 
the  end  of  the  interval.  Subroutine  PRESS  is  called  to  deter- 
mine the  fluid  friction  pressure  losses  in  each  of  the  piping 
loops  for  the  fluid  flowrates  existing  during  the  time  step, 
and  Subroutine  FLOWS  adjusts  these  flowrates  at  the  end  of 
the  interval,  based  upon  a comparison  between  the  actual 
circulation  pump  pressure  settings  and  the  computed  pressure 
losses.  The  instantaneous  fluid  flowrate  in  the  TUS  side 
of  the  power  plant  thermal  supply  heat  exchanger  is  calculated 
so  as  to  transfer  the  available  thermal  energy  to  the  TUS 
storage  reservoir,  and  the  hot  and  cold  water  storage  volumes 
are  adjusted  through  a comparison  between  the  inter-section 
flowrate  through  the  supply  heat  exchanger  and  the  flowrate 
through  the  primary  HTW  loop.  When  the  program  control  is 
returned  to  Program  MAIN,  the  new  fluid  temperatures  are 
compared  with  the  reference  values  at  the  load  center  heat 
exchangers  and  at  the  outlets  of  the  primary  loop  heat  exchan- 
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gers.  If  none  of  the  temperature  limits  has  been  exceeded, 
the  time  clock  is  incremented  and  the  linear  interpolation 
procedure  is  begun  to  determine  the  demands  for  the  next 
time  step.  However,  if  any  of  the  temperature  limits  has 
been  exceeded,  the  thermostat-controlled  fluid  flowrate 
correction  process  is  begun  by  calculating  for  the  entire 
TUS  a set  of  desired  fluid  mass  flowrates  which  will  exactly 
satisfy  the  given  instantaneous  thermal  demands  at:  (1)  the 

temperature  limits  for  the  load  center  heat  exchangers,  and 
C2 } the  desired  fluid  reference  temperatures  at  the  primary 
loop  heat  exchangers.  The  computed  parallel  branch  flowrates 
are  added  to  form  desired  total  loop  flowrates  for  each  of 
the  secondary  distribution  loops,  and  these  flowrates  are 
used  in  the  calculation  of  the  required  flow  splits  among 
the  primary  loop  heat  exchangers.  The  primary  loop  nodal 
fluid  flowrates  computed  in  this  manner  for  each  of  the 
secondary  supply  heat  exchangers  will,  in  general,  be  unique- 
ly defined  for  each  unit.  To  preserve  continuity  of  the 
total  primary  loop  flowstream,  therefore,  the  heat  exchanger 
with'-the  largest  computed  flowrate  is  assumed  to  have  no 
bypass,  and  bypass  flows  are  established  around  each  of  the 


other  units  based  upon  the  difference  between  this  maximum 
and  the  desired  unit  flows.  The  desired  heat  exchanger/ 
bypass  flowrate  ratio  is  then  used  to  scale  the  existing 
nodal  flowrates  to  effectively  set  the  required  flow  splits 
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(but  not  the  magnitudes  of  the  flows)  at  the  end  of  the 
time  step  being  studied.  The  computed  desired  fluid  flow- 
rates  are  transmitted  .to  Subroutines  PRESS  and  PLOWS,  in 
which  a set  of  corresponding  desired  circulation  pump  pres- 
sure settings,  scaled  by  the  user-defined  pump  response 
sensitivity  factor,  are  calculated  and  are  applied  to  the 
system  pumps.  The  existing  system  fluid  flowrates  are  then 
re-calculated  for  the  given  time  step  based  upon  these  re- 
vised pressure  settings.  (Waiting  until  the  end  of  the 
following  time  step  to  update  the  flows  could  result  in  an 
overall  inability  to  adequately  follow  large  load  variations 
occurring  during  a single  time  interval.)  Following  the  re- 
setting of  the  pump  pressures  and  the  re-calculation  of 
the  fluid  flowrates,  the  thermostat  model  returns  normal 
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program  flow  control  to  Program  MAIN,  the  time  clock  is 
incremented,  and  the  next  time  step's  calculations  are  begun. 

Output  from  the  program  during  the  dynamic  system  ana- 
lysis is  controlled  by  a user-defined  printout  frequency 
code.  At  regular  intervals  (as  frequently  as  every  time 
step  or  as  infrequently  as  once  each  hour)  an  instantaneous 
TUS  summary  is  printed  which  includes  all  the  nodal  fluid 
flowrates,  temperatures,  and  energy  demands,  the  total 
system  thermal  and  electrical  energy  supply  and  demand  condi- 
tions, the  instantaneous  circulation  pump  pressure  settings, 
and  the  thermal  energy  storage  reservoir  fluid  volume  and 
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to  yield  successively  more  detailed  specifications  of  the 
community’s  demands  and  of  the  energy  system  design  para- 
meters. Some  general  guidelines  which  have  been  found  to 
be  helpful  in  the  development  of  such  an  energy  system  model 
are : 

(1)  Identify  a set  of  up  to  20  different  building 
classifications  which  span  the  range  of  architec- 
tural and  occupancy  characteristics  observed  in 
the  community. 

(2)  For  each  classification,  define  a representative 
building  unit  which  exhibits  the  typical  structural 
and  usage  characteristics  of  its  class.  This 
building  may  be  either  an  existing  unit  whose  para- 
meters are  deemed  typical,  or  it  may  be  a repre- 
sentative model  composed  of  the  aggregate  character- 
istics of  several  individual  units. 

(3)  Distribute  these  representative  building  units 
throughout  the  community  according  to  the  observed 
distribution  of  actual  buildings.  The  total  number 
of  representative  units  specified  at  each  location 
in  the  community  should  be  chosen  to  approximate 
the  total  floor  area  (or  some  other  equivalent 
energy  demand  index)  of  the  actual  buildings  at 
that  location. 

(4) .  Aggregate  groups  of  buildings  into  energy  demand 
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load  centers,  based  upon  their  number  density  and 
similarity  of  use.  In  sizing  these  load  centers, 
it  should  be  remembered  that  each  one  connected 
to  the  TUS  will  be  served  by  a single  heat  exchan- 
ger. 

(5)  Design  a proposed  TUS  piping  layout  which  will 
supply  HTW  to  the  desired  load  centers.  (Not 
necessarily  every  load  center  will  be  served  by 
the  TUS,  since  the  economically  optimum  TES  may 
require  some  buildings  to  be  supplied  all-electri- 
cally.)  Each  load  center  should  be  served  by  a 
single  heat  exchanger,  up  to  five  of  which  should 
be  piped  together  in  parallel  to  form  a single 
secondary  distribution  loop.  A maximum  of  five 
secondary  loops  is  allowed  in  the  system. 

(6)  Once  this  rough  system  design  has  been  specified, 
estimates  of  pipe  sizes  and  lengths,  pumps  and  heat 
exchanger  locations,  and  pipe  insulation  parameters 
should  be  coded  into  the  TDIST  TUS  input  data  for- 
mats (see  Section  A. 5).  The  distribution  (number 
of  units  served  by  each  load  center  heat  exchanger 
and  total  number  of  units  not  served  by  the  TUS) 
and  specifications  of  each  representative  building 
type  should  be  coded  into  the  TDIST  building  input 
data  formats. 
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(7)  Specify  the  hourly  weather  data  to  be  used  during 
the  simulation  period. 

(8)  Use  the  building  design  analysis  run  option 
(KODE  = 0,  IDES  = 1,  KODE  = 1)  to  provide  the 
space  conditioning  load  profile  for  each  building 
type  over  the  simulation  period  and  to  aggregate 
the  total  energy  demands  to  be  supplied  by  the  TUS 
at  the  beginning  of  the  period. 

(9)  Examine  the  computed  load  profiles  to  identify  any 
questionable  or  invalid  building  data  specifica- 
tions (i.e.,  if  the  load  profile  is  not  "reason- 
able", identify  the  controlling  load  component 
and  the  building  parameters  determining  that  com- 
ponent). Any  parameter  revisions  should  be  veri- 
fied by  re-running  the  code  and  generating  a correc- 
ted set  of  loads. 

(10)  If  necessary,  due  to  significant  imbalances  in  the 
thermal  load  distribution  within  a given  secondary 
loop,  the  load  centers  should  be  redefined  or 
rearranged  to  yield  reasonably  uniform  loads  among 
the  heat  exchangers  in  each  loop,  and  the  code 
should  be  re-run  to  generate  the  revised  load 
distribution . 


(11)  Once  the  building  loads  have  been  verified,  use  the 
given  demand  distribution  and  specifications  of  the 


desired  TUS  fluid  temperatures  to  size  the  TUS 
heat  exchangers.  In  general,  it  is  advisable  to 
run  the  code  twice  during  this  sizing  exercise, 
unless  the  peak  hourly  demands  occur  at  the  begin- 
ning of  the  simulation  period.  One  run  (KODE  = 0, 
IDES  = 1,  KODE  = 1,  as  summarized  in  steps  8-10) 
is  used  to  generate  the  hourly  building  demands 
and  the  TUS  loads  to  be  used  in  specifying  the  re- 
quired initial  system  temperatures  and  flowrates; 
the  second  run  (KODE  = 0,  IDES  = 1,  KODE  = 2)  is 
used  to  determine  the  single  time  instant  TUS  peak 
demands  to  be  used  in  sizing  the  heat  exchangers. 

(12)  After  specifying  the  heat  exchanger  designs,  perform 
the  steady-state  analysis  (KODE  = 0,  IDES  = 0)  of 
the  system  at  the  beginning  of  the  simulation 
period  to  insure  that  the  system  design  is  self- 
consistent  and  that  the  system  will  converge  to  the 
desired  fluid  flowrates  and  temperatures. 

(13)  If  necessary,  make  design  revisions  (such  as  chang- 
ing heat  exchanger  ratings,  the  initial  primary 
HTW  loop  fluid  mass  flowrate,  or  the  sizes  of  the 
distribution  loop  piping)  and  re-run  the  steady- 
state  analysis  until  the  desired  fluid  flowrate, 
pressure  loss,  and  fluid  temperature  conditions  are 


achieved . 
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(14)  Perforin  the  time-dependent  TES  simulation  by  call- 
ing the  steady-state  analysis  option,  followed 
by  the  dynamic  system  analysis  (KODE  = 0,  IDES  = 0), 
KODE  = 1) . 

Expansion  of  the  run  options  in  this  manner  from  the  build- 
ing load  analysis  to  the  full  dynamic  system  simulation  allows 
the  user  to  clearly  define  the  community  loads  (which  are, 
of  course,  independent  of  the  assumed  method  of  energy  sup- 
ply) prior  to  performing  any  TES  analyses  other  than  choos- 
ing a proposed  piping  layout.  Once  the  proper  building  models 
have  been  specified  and  the  appropriate  TUS  configuration  has 
been  chosen,  the  resulting  heat  exchanger  designs  and  pro- 
posed fluid  flowrate  and  temperature  specif ications  may  be 
fully  tested,  and  any  required  system  design  revisions  may 
be  made  prior  to  performing  the  dynamic  system  analysis. 
Following  full  specification  of  the  system  design,  the  dynamic 
analysis  run  option  provides  information  as  to  how  the  parti- 
cular design  chosen  behaves  over  time.  Although  possibly 
seeming  over-conservative  in  its  approach  to  system  design, 
this  methodology  allows  a logical  system  development  and, 
based  upon  the  computation  time  estimates  outlined  in  the 
following  section,  provides  a cost-effective  sequence  of 
system  analysis. 
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A. *4. 5 Computer  Execution  Time  Specifications 

TDIST  is  written  in  FORTRAN  IV  and  is  directly  compa- 
tible with  the  IBM  System  370  computer  using  the  FORTRAN  H- 
level  compiler.  Program  compilation  requires  the  allocation 
of  288  K bytes  of  working  storage;  execution  of  the  compiled 
load  module  requires  128  K bytes.  The  source  deck  consists 
of  approximately  1800  cards.  The  compiled  load  module  may 
be  stored  on  a total  of  8 tracks  of  disk  space. 

When  using  a pre-compiled  load  module  called  from  re- 
serve storage,  the  user  may  expect  execution  times  for  the 
various  program  options  to  approximate  the  values  listed  in 
Table  A. 10.  In  each  case,  the  time  shown  is  the  total 
elapsed  time  for  the  run,  including  program  set-up  and  block- 
ing of  the  output  data  for  subsequent  printing;  all  times  are 
in  CPU-minutes.  The  lower  limit  of  the  range  corresponds  to 
the  analysis  of  a TUS  containing  two  secondary  distribution 
loops,  one  serving  three  and  the  other  serving  five  load 
center  heat  exchangers.  The  upper  limit  corresponds  to  the 
analysis  of  a TUS  having  four  secondary  loops,  serving  a 
total- of  18  load  center  heat  exchangers.  A total  of  11 
building  types  are  specified  for  each  system.  Steady-state 
system  input  data  includes  estimates  of  the  initial  secondary 
loop  fluid  flowrates,  which  serve  to  reduce  the  total  conver- 
gence time  required  for  each  system.  The  time  step  size  for 
the  dynamic  analyses  is  0.001  hour,  and  each  simulation 
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TABLE  A. 10 

REPRESENTATIVE  OPTION  EXECUTION  TIMES 


Execution  Time 

Program  Option  (CPU  - minutes) 


Single-time  building  load  analysis 
(KODE=0,  IDES=1,  KODE=2 ) 

■» 

0.096  - 0.112 

2^-hour  building  load  analysis 
(KODE=0,  IDES=1 , K0DE=1 ) 

0.110  - 0.123 

Steady-state  system  analysis 
(K0DE=0 , IDES=0 ) 

0.105  - 0.108 

24-hour  dynamic  simulation 
(DTIME  = 0.001) 

(KODE=0,  IDES-0,  K0DE=1) 

4.374  - 8.691 

Program  compilation  and  load  module 
storage 

0.692 

covers  a period  of  2b  hours.  The  dynamic  simulations  are 
initialized  internally  by  first  calling  the  steady-state 
system  analysis  option.  Data  output  for  the  dynamic  analyses 
occurs  at  every  half  hour. 

A. 5 TDIST  Input  Data 

All  input  to  TDIST  is  from  cards  supplied  by  the  user. 
Since  fixed-format  input  statements  are  used  throughout  the 
program,  it  is  important  to  follow  the  data  format  instruc- 
tions outlined  below.  The  data  which  must  be  specified  for 
each  program  run*  consists  of  four  general  categories; 

(1)  specifications  of  the  physical  parameters,  heat  transfer 
coefficients,  and  average  fluid  conditions  in  each  TUS  node, 

(2)  specifications  of  the  heat  transfer  characteristics,  time- 
dependent  usage,  and  distribution  of  each  building  type, 

(3)  hourly  weather  data,  and  (4)  program  run  control  and 
initialization  codes.  To  some  extent,  these  categories  are 
segregated  within  the  input  data  deck  for  ease  of  parameter 
modification.  However,  the  segregation  is  by  no  means 
complete,  and  care  must  be  taken  to  formulate  the  deck  in  the 

to  , „ 

exact  order  presented.  Input  diagnostics  exist,  but  are  mini- 
mal and  are  only  available  during  runs  in  which  Subroutine 

*Input  data  flexibility  is  extremely  limited.  Depending  upon 
the  type  of  run  to  be  performed,  some  data  values  may  be  left 
unspecified,  but  their  locations  in  the  data  deck  must  be 
filled  with  blanks  or  dummy  values  to  facilitate  correct 
input  of  all  the  required  information  (see  Section  A. 5. 30 


The  primary  HTW  loop  nodes  include  the  pipes  supplying  heat 
from  and  returning  heat  to  the  thermal  energy  storage  reser 
voir,  but  do  not  include  the  reservoir  sections  or  the  TUS 
supply  heat  exchanger.  Pipe  runs  between  heat  exchangers 
within  the  loop  must  be  included,  no  matter  how  short,  to 
result  in  effective  calculation  of  losses.  A maximum  of  13 
nodes  is  allowed,  including  a maximum  of  5 heat  exchangers. 

Variables:  NH  - number  of  primary  HTW  loop  nodes 

Format:  c.c.  1-3  13  NH 

3 . Secondary  Loop  Number  Card 

This  card  inputs  the  number  of  secondary  distribu- 
tion loops  in  the  system.  Since  each  of  these  loops  is 
supplied  from  a heat  exchanger  in  the  primary  HTW  loop,  the 
number  of  secondary  loops  must  correspond  to  the  number  of 
heat  exchanger  nodes  in  the  primary  HTW  loop.  A maximum  of 
5 secondary  distribution  loops  is  allowed. 

Variables:  LS  - number  of  secondary  distribution 

loops 

Format:  c.c.  1-3  13  LS 

• 

4 . Secondary  Loop  and  Consumer  Branch  Node  Number 
Card 

This  card  inputs  the  number  of  nodes  in  a secondary 
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distribution  loop  and  the  number  of  load  center  heat  exchan- 
gers supplied  from  the  loop.  The  secondary  loop  nodes 
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5 . Node  Dimension  and  Location  Code  Card 
This  card  inputs  the  node  identification  and  loca- 
tion codes  for  any  of  the  nodes  in  the  system  and  the  physi- 
cal dimensions  and  heat  transfer  coefficient  for  the  pipe  or 
heat  exchanger  represented.  One  card  must  be  supplied  for 
each  of  the  nodes  specified  in  the  system. 

Basic  node  identification  is  accomplished  through  the 
use  of  three  indicators,  generally  designated  in  the  program 
as  K,  L,  and  N.  K identifies  the  major  subsection  of  the 
distribution  system  in  which  the  node  appears:  1 = primary 

MTW  loop,  2 = secondary  distribution  loops,  3 = consumer 
sides  of  load  center  heat  exchangers.  L indicates  the  speci- 
fic loop  in  which  the  node  appears:  for  the  primary  HTW 

loop  L = 1;  each  of  the  secondary  distribution  loops  is 
numbered  in  the  order  in  which  it  is  supplied  from  the  pri- 
mary loop,  L ranges  from  1-5;  the  load  center  heat  exchangers 
supplied  from  a given  secondary  loop  are  coded  with  the  L 
number  of  that  loop.  N is  the  specific  node  number  within  a 
loop:  in  the  primary  HTW  loop,  node  1 is  the  outlet  pipe 

fi^om  the  hot  water  section  of  the  storage  reservoir  and 
the  last  node  is  the  return  pipe  to  the  co_d  water  section 
from  the  final  secondary  supply  heat  exchanger  (maximum  of 
N = 13  nodes);  in  each  secondary  distribution  loop,  node  1 
is  the  outlet  pipe  from  the  primary/secondary  heat  exchanger 
supplying  the  loop  and  the  last  node  is  that  heat  exchanger; 
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parallel  branch  nodes  are  numbered  consecutively  according 
to  their  schematic  representation  (see  Pig.  A. 10)  (maximum  of 
N = 13  nodes  per  loop).  The  load  center  heat  exchangers 
served  by  each  secondary  loop  are  numbered  consecutively 
according  to  their  schematic  representation  (maximum  of  5 
heat  exchangers  per  loop).  Following  this  convention,  a 
nodal  parameter  with  the  designation  X(2,3,8)  would  be  asso- 
ciated with  node  number  8 in  the  third  secondary  distribution 
loop;  Y(3,2,4)  would  correspond  to  the  fourth  load  center 
heat  exchanger  served  by  the  second  secondary  loop;  Z(l,l,l) 
would  apply  to  the  first  primary  HTW  loop  node. 

The  identification  codes  on  this  card  must  include  the 
N for  the  node  being  specified  (K  and  L are  supplied  inter- 
nally according  to  the  order  in  which  the  nodal  data  is  sup- 
plied), an  indicator  for  whether  the  node  is  a series  flow  or 
parallel  branch  node  (all  primary  loop  nodes  are  series,  the 
parallel  secondary  loop  branches  and  load  center  heat  exchan- 
gers are  parallel),  an  indicator  for  whether  or  not  the  node 
is  a heat  exchanger,  and  an  indicator  for  whether  or  not  the 
node  contains  the  circulation  pump  for  the  loop  (one  pump 
node  per  loop  is  allowed). 

Dimensional  data  which  must  be  supplied  includes  the 
nodal  flow  channel  diameter  (in  feet)  and  the  physical  length 
of  the  node  (in  feet).  For  pipe  nodes,  the 'flow  channel 
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Secondary  Loop  and  Load  Center  Heat  Exchanger 
Node  Numbering  Conventions 


diameter  is  the  specified  I.D.  of  the  pipe;  for  heau  exchan- 
gers, the  flow  channel  diameter  to  be  specified  is  the  I.D. 
of  a single  tube  or  channel  between  the  tubes  (not  the 
overall  heat  exchanger  shell  diameter).  The  length  of  the 
node  is  either  the  total  pipe  run  or  the  lateral  dimension 
of  the  heat  exchanger  tube  bundle. 

The  heat  transfer  coefficient  specified  must  be  the 
overall  heat  transfer  coefficient  for  the  node  given  in  units 
of  BTU/hr-°F  (i.e.,  the  "UA"  value  for  the  node).  For  a 
pipe  node,  this  should  include  the  heat  transfer  character- 
istics of  the  fluid  film,  the  pipe  material,  any  insulation, 
and  the  thermal  resistance  of  the  soil  in  which  the  pipe  is 
buried.  The  UA  value  specified  for  a heat  exchanger  must  be  the 
total  value  for  the  heat  exchanger  assuming  a single-pass 
counterflow  geometry.  Since  this  value  is  assumed  to  remain 
constant  throughout  the  simulation  and  is  assumed  to  be  uni- 
form along  the  length  of  the  node,  it  should  be  specified  at 
the  design  average  fluid  flowrate  and  temperature  conditions 
for  the  node . 

Variables:  N - node  number  code  (1-13  for  primary  loop, 

1-13  for  secondary  loop,  1-5  for  load 
center  heat  exchangers) 

ISER(K,L,N)  - series/parallel  node  indicator  (0  = 
series  node,  1 = parallel  node) 


c . c . 

1-3 

13 

c . c . 

4-6 

13 

c . c . 

7-9 

13 

c . c . 

10-12 

13 

c . c . 

13-19 

F7 . 4 

c . c . 

20-30 

Ell. 4 

c . c . 

31-41 

Ell.  4 

ISER(K,L,N) 

IEX(K,L,N) 

IPMP ( K , L , N ) 

DIA(K,L,N) 

LONG(K,L,N) 

AUF(K,L,N) 


6 . Nodal  Fluid  Parameter  Information  Card 
This  card  inputs  the  initial  fluid  mass  flowrate, 
nodal  fluid  inlet  and  outlet  temperatures,  and  average  fluid 
specific  heat,  density  and  viscosity  for  any  of  the  nodes  in 
the  system.  One  card  must  be  supplied  for  each  node. 
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The  nodal  fluid  mass  flowrate  must  be  specified  in 
units  of  lbm/hr.  Nodal  values  supplied  should  correspond  to 
the  design  values  for  the  system  initial  conditions  and  should 
be  consistent  with  the  flow  channel  geometry  within  each  of 
the  loops  (e.g.,  the  series  primary  loop  nodes  should  each 
have  the  same  mass  flowrate  to  preserve  continuity  of  the 
flow).  In  all  cases,  the  primary  loop  flowrate  must  be  given 
as  accurately  as  possible,  since  it  serves  as  basic  input  to 
the  steady-state  convergence  subroutine.  If  the  steady-state 
subroutine  is  to  be  used  for  system  initialization,  the  flow- 
rates  for  the  secondary  loop  and  load  center  heat  exchanger 
nodes  need  not  be  specified  precisely,  since  they  are  computed 
internally  during  the  system  convergence  calculations.  How- 
ever, estimates  of  these  flowrates  are  used  to  initialize 
the  steady-state  calculations,  and  the  data  provided  should 
be  approximately  equal  to  the  expected  values  to  speed  the 
system  convergence. 

If  the  steady-state  analysis  is  used  to  initialize  the 
dynamic  system  calculations,  the  only  temperatures  required  as 
input  data  are  the  inlet  temperature  to  the  first  node  in 
the  primary  HTW  loop  (which  becomes  the  constant  TUS  fluid 
supply  temperature)  and  the' supply  and  return  temperatures 
on  the  consumer  sides  of  the  load  center  heat  exchangers. 

These  consumer  temperatures  should  be  set  according  to  speci- 
fications for  the  end-use  equipment  served,  and  the  actual 


•I 


123- 

temperatures  are  allowed  to  vary  from  these  design  values 
during  the  simulation  according  to  the  internally-defined 
thermostat  control  limits. 

The  specified  nodal  fluid  specific  heat,  density  and 

viscosity  should  correspond  to  design  values  calculated 

from  the  average  fluid  temperature  conditions  for  the  node. 

The  specific  heat  must  be  given  in  units  of  BTU/lbm-°F,  the 

•3 

density  in  lbm/ft  , and  the  viscosity  in  lbm/ft-hr.  These 
values  are  assumed  to  be  constant  throughout  the  simulation 
and  uniform  along  the  length  of  the  node;  discontinuities  at 
nodal  boundaries  are  allowed  to  account  for  variations  in 
the  average  nodal  fluid  temperature  around  each  of  the  loops. 
Since  these  parameters  are  used  in  the  determination  of  the 
specific  nodal  heat  and  fluid  momentum  transfer  rates,  the 
values  specified  should  be  as  precise  as  possible  wituin  the 
nodal  averaging  constraints  outlined  above. 

Variables : 

MASS(K,L,N)  - nodal  fluid  mass  flowrate  (in  lbm/hr) 
TI(K,L,N)  - nodal  fluid  inlet  temperature  (in  °F) 
TO(K,L,N)  - nodal  fluid  outlet  temperature  (in  °F) 
CAP(K,L,N)  - nodal  fluid  specific  heat  (in  BTU/lbm-°F) 
RHO ( K , L , N ) - nodal  fluid  density  (in  lbm/ft3) 

MU(K,L,N)  - nodal  fluid  viscosity  (in  lbm/ft-hr) 
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Format : 

c.c. 

1-11 

Ell. 4 

MASS (K,L,N) 

c.c. 

12-18 

F7.2 

TI(K,L,N) 

c.c. 

19-25 

F7.2 

TO  ( K , L , N ) 

c.c. 

26-32 

F7.4 

CAP(K,L,N) 

c.c. 

33-39 

F7.4 

RH0(K,L,N) 

c.c. 

40-46 

F7.4 

MU(K,L,N) 

7. 

Heat 

Exchanger  Number 

Card 

This  card  inputs  the  number  of  heat  exchangers  in 
the  HTW  distribution  system.  Included  should  be  the  secon- 
dary supply  and  the  load  cer ter  heat  exchangers,  but  not  the 
TUS  supply  heat  exchanger.  A maximum  of  31  heat  exchangers 

is  allowed. 

Variables : 

NEXCH  - number  of  heat  exchangers  in  the  distribution 
system 

Format:  c.c.  1-3  13  NEXCH 

8 . Heat  Exchanger  Data  Card 

This  card  inputs  the  heat  exchanger  code  number, 
the  identification  codes  for  its  supply  ar.d  demand  side  nodes, 
and  the  number  of  parallel  flow  tubes  it  contains. 

The  system  heat  exchangers  are  numbered  consecutively 
from  the  first  heat  exchanger  in  the  primary  HTW  loop  to  the 
last  load  center  heat  exchanger.  The  code  number  for  .the 
heat  exchanger  given  on  this  card  should  correspond  to  the 


j 

4 


rLi* 


I 
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values  for  IEX(K,L,N)  given  on  the  two  Node  Dimension  and 
Location  Cards  containing  its  dimensions  and  heat  transfer 
characteristics.  A maximum  of  31  heat  exchangers  are  allowed 
in  the  system. 

For  the  purposes  of  heat  exchanger  identification,  the 
supply-side  nodes  for  the  primary/secondary  heat  exchangers 
are  the  primary  HTW  loop  nodes  serving  them;  the  supply-side 
nodes  for  the  load  center  heat  exchangers  are  the  parallel 
secondary  loop  nodes  serving  them. 

The  number  of  tubes  in  a heat  exchanger  should  be  con- 
sistent with  the  specified  heat  transfer  coefficient,  physi- 
cal dimensions,  single-pass  counterflow  geometry,  and  design 
fluid  flow  conditions  for  the  heat  exchanger. 

Variables : 

NEX  - heat  exchanger  code  number 

Kl(NEX), Ll(NEX) , Nl(NEX)  - heat  exchanger  supply-side 
identification  codes  [NEX  = IEX(K1,  LI,  Nl)] 

K2(NEX),  L2(NEX),  N2(NEX)  - heat  exchanger  demand- 

i 

side  node  identification  codes  [NEX  = IEX(K2,L2  , N2)  ] 

NTB(NEX)  - number  of  parallel  flow  tubes  in  heat  ex- 


_ J 


changer 


! 
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Format : 


c.c. 

1-3 

13 

NEX 

c.c. 

9-11 

13 

Kl(NEX) 

c.c. 

12-14 

13 

LI (NEX) 

c.c. 

15-17 

13 

Nl(NEX) 

c.c. 

23-25 

13 

K2 (NEX) 

c.c. 

26-28 

13 

L2 (NEX) 

c.c. 

29-31 

13 

N2 (NEX) 

c.c. 

37-41 

15 

NTB(NEX) 

9 . Circulation  Pump  Number  Card 
This  card  Inputs  the  number  of  circulation  pumps 
in  the  distribution  system.  To  simplify  the  loop  pump  pres- 
sure requirement  analysis  and  fluid  acceleration  calculations, 
only  one  pump  per  loop  is  allowed.  The  pumps  must  be  located 
in  pipe  nodes  as  indicated  by  the  IPMP(K,L,N)  code  on  the 
Node  Dimension  and  Location  Cards.  A maximum  of  6 circula- 
tion pumps  are  allowed  in  the  system. 

Variables : 

NPMP  - number  of  circulation  pumps  in  the  distribution 
system 


Format:  c.c.  1-3  13  NPMP 


10 . Storage  Reservoir  Card 

This  card  inputs  the  initial  hot  and  cold  water 
storage  reservoir  section  volumes  and  the  initial  hot  water 
storage  temperature.  The  volumes  must  be  specified  in  units 


of  ft  and  the  temperature  in  °F. 

Since  variations  in  the  hot  water  section  volume  during 
the  simulation  determine  the  required  storage  capacity  for 
system  load  smoothing,  the  exact  values  of  the  initial  vol- 
umes are  not  critical  to  the  system  heat  flow  calculations. 
However,  the  volumes  should  be  specified  large  enough  to 
preclude  emptying  of  either  section  during  the  time  period 
under  investigation.  The  reservoir  sections  are  assumed  to 
be  located  between  the  TUS  supply  heat  exchanger  and  the 
primary  loop  supply  and  return  pipes;  no  node  identification 
codes  are  required  for  the  reservoir. 

The  initial  hot  water  storage  temperature  should  corres- 
pond to  the  specified  fluid  inlet  temperature  for  the  first 
primary  loop  node.  This  temperature  is  held  constant  through- 
out the  simulation  unless  the  hot  water  reservoir  section  is 
emptied . 

Variables : 

HVO  - hot  water  storage  section  initial  volume  (in  ft^) 
CVO  - cold  water  storage  section  initial  volume  (in  ft^) 
THW  - hot  water  storage  temperature  (in  °F) 


Format : 

c . c . 

1-11 

Ell.  4 

HVO 

c . c . 

12-22 

Ell. 4 

CVO 

c . c . 

23-29 

F7.2 

THW 

11 . TUS  Supply  Heat  Exchanger  Data  Card 
This  card  inputs  the  power  plant  coolant  inlet 
and  outlet  temperatures,  the  initial  fluid  mass  flowrate, 
and  the  fluid  specific  heat  on  the  supply-side  of  the  TUS 
supply  heat  exchanger.  The  temperatures  must  be  specified 
in  °F,  the  flowrate  in  lbm/hr,  and  the  specific  heat  in 
BTU/lbm-°F. 

The  supply  side  inlet  and  outlet  temperatures  are  held 
constant  throughout  the  simulation  unless  the  heat  flow  across 
the  heat  exchanger  demands  that  either  one  or  both  be  modi- 
fied to  comply  with  the  Second  Law  of  Thermodynamics.  The 
inlet  temperature  should  be  specified  high  enough  to  provide 
heat  transfer  to  the  primary  HTW  loop  inlet  fluid;  the  outlet 
temperature  should  be  as  low  as  possible  for  efficient  power 
plant  operation,  yet  high  enough  to  provide  heat  transfer  to 
the  primary  loop  return  fluid. 

The  exact  value  for  the  coolant  mass  flowrate  required 
to  transfer  the  available  heat  within  the  specified  tempera- 
ture  limits  is  computed  internally  during  the  initialization 
of  the  system.  However,  a non-zero  value  for  this  parameter 
should  be  supplied  to  the  program  to  avoid  the  possibility  of 
computational  overflow  errors  during  this  portion  of  the 
simulation. 

The  specific  heat  of  the  fluid  should  correspond  to  the 


value  calculated  for  the  average  fluid  temperature  in  the 
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heat  exchanger;  the  specific  heat  remains  unchanged  during 
the  simulation. 

Variables : 

THG  - coolant  inlet  temperature  (°F) 

TCG  - coolant  outlet  temperature  (°F) 

MGAS  - coolant  mass  flowrate  (lbm/hr) 

CAPG  - coolant  specific  heat  (BTU/lbm-°F) 


Formate : 

c.c. 1-7 

F7.2 

THG 

c.c. 8-14 

F7.2 

TCG 

c.c. 15-25 

Ell. 4 

MGAS 

c.c. 26-32 

F7  • 4 

CAPG 

12 . Ambient  Earth  Temperature  Card 

This  card  inputs  the  ambient  earth  temperature  at 
the  average  burial  depth  for  the  distribution  piping.  It  is 
useu  in  the  calculation  of  heat  losses  from  buried  pipes  and 
conduction  losses  to  the  earth  from  buildings.  The  earth 
tamperature  must  be  given  in  °F. 

Variables : 

TAMB  - ambient  earth  temperature  (in  °F) 

Format:  c.c.  1-7  F7.4  TAMB 

13-  Power  Plant  Thermal/Electric  Energy  Conversl 
Efficiency  Card 

This  card  inputs  the  overall  power  plant  thermal/' 
electric  energy  conversion  efficiency  as  the  ratio  of  the 
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electrical  energy  output  from  the  generators  to  the  thermal 
energy  content  of  the  fuel  burned.  It  is  used  in  determining 


the  total  energy  consumption  for  electricity  production  in 
conjunction  with  the  heat  available  to  the  thermal  system. 

The  efficiency  should  be  given  as  a decimal  fraction. 
Variables : 

EFP  - power  plant  overall  thermal/electric  energy 
conversion  efficiency  (decimal  fraction) 

Format:  c.c.  1-7  F7 . 4 EFF 

14 . Pump  Response  Sensitivity  Factor  Card 

This  card  inputs  the  circulation  pump  response 
sensitivity  factor  for  all  the  pumps  in  the  system.  This 
factor  is  applied  to  the  difference  between  the  projected 
pressure  settings  and  the  existing  pump  pressures  to  provide 
an  accelerated  response  to  the  loop  fluid  flowrates  when  the 
thermostat  is  activated.  A factor  of  F results  in  the  pumps 
being  set  at  pressures  equal  to  the  required  pressure  plus 
F times  the  difference  between  the  required  and  existing 
pressures.  Therefore,  a factor  of  2.0  results  in  settings 
of: 

P = 3P  - op 

set  J required  existing 

Higher  factors  result  in  faster  response  times;  very  fast 

response  times  should  be  avoided  since  they  can  result  in 

overshooting  of  the  desired  flow  conditions  and  instabilities 

in  subsequent  heat  flow  calculations.  The  factor  is  dropped 

from  the  pump  settings  when  the  loop  flows  have  accelerated 


A 


: 1 

t 1 


to  within  1%  of  the  desired  values. 

Variables : 

PSENS  - pump  response  sensitivity  factor 


Format:  c.c.1-7 


F7.4 


PSENS 


15 • Time  Increment  and  Print  Interval  Data  Card 
This  card  inputs  the  number  of  time  increments  per 
hour  to  be  used  in  the  transient  calculations,  the  number  of 
increments  between  system  information  printout,  and  the  size 
of  one  time  step.  The  number  of  increments  per  hour  and 
the  number  of  increments  between  printouts  should  be  given 
as  integer  numbers;  the  number  of  increments  per  hour  must 
be  a whole  number  multiple  of  the  number  of  increments  per 
printout.  A minimum  of  1 printout  per  hour  is  allowed  (i.e., 
the  number  of  steps  between  printouts  = the  number  of  steps 
per  hour).  The  time  step  si.e  should  be  given  as  a decimal 
fraction  of  an  hour. 

A small  time  step  size  results  in  increased  accuracy  in 
the  dynamic  calculations  by  allowing  a closer  modelling  of 
the  continuously  varying  system  conditions.  However,  short 
time  steps  also  result  in  considerable  increases  in  the 
total  computational  time  for  any  given  simulation.  Conversely, 
longer  time  steps  reduce  both  computation  time  and  modeling 
accuracy.  For  each  simulation,  a compromise  value  must  be 
chosen  which  reflects  the  more  important  of  the  computational 


parameters  desired  for  the  given  run. 


L 


_ 


Although  it  is  possible  to  print  the  entire  set  of  sys- 
tem parameters  for  every  time  step  in  the  simulation,  it  is 
advisable  both  for  run  time  and  data  reduction  purposes  to 
print  out  this  data  only  at  a frequency  high  enough  to  give 
meaningful  output.  If  the  system  excitation  conditions  vary 
greatly  from  hour  to  hour,  a print  interval  of  once  every 
quarter  hour  may  be  necessary  to  yield  the  desired  detail  in 
the  output;  for  slowly  varying  conditions,  an  output  listing 
once  every  hour  could  be  sufficient. 

Variables : 

NSTEP  - number  of  time  steps  per  hour 

NPRNT  - number  of  time  steps  between  printouts 


DTI  ME 

- single  time 

step  size 

(as  a dec: 

tion 

of  an 

hour ) 

Format : 

c . c . 

1-5 

15 

NSTEP 

c . c . 

6-10 

15 

NPRNT 

c . c . 

11-18 

F8 . 5 

DTIME 

1 6.  Building  Type  Number  Card 

This  card  inputs  the  number  of  different  building 
types  to  be  served  by  the  total  energy  system.  For  the  pur- 
poses of  determining  this  number,  a distinct  building  type 
is  defined  as  one  having  representative  physical  dimensions, 
occupancy,  room  temperature  requirements,  and  materials  heat 
transfer  characteristics.  Buildings  used  for  similar  purposes 
and  having  similar  physical  characteristics  should  be  placed 
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in  a single  type  class,,  even  though  there  may  be  difference.: 
in  their  detailed  construction  or  layout.  A maximum  of  20 
different  building  types  is  allowed. 

Variables : 

NTYPE  - number  of  building  types  served  by  system 
Format:  c.c.  1-3  13  NTYPE 

17  • Bulding  Material  Heat  Transfer  Parameter  Cl-::  . 

This  card  inputs  the  exposed  areas  and  thermal 

resistances  of  the  walls,  windows,  roof,  and  ground-contact 

material  for  the  building  type  specified  by  its  building 

2 

code  number.  The  areas  must  be  given  in  ft  and  the  thermal 
resistances  in  (BTU/hr-ft  -°F)  . One  card  must  be  input  f 

each  building  type  as  specified  on  the  Building  Type  Number 
Card . 

The  building  type  code  number  is  simply  an  integer 
which  identifies  the  particular  building  type  being  specific" 

Types  should  be  numbered  consecutively  for  any  given  Simula 
tion;  a maximum  of  20  building  types  is  allowed. 

The  areas  specified  should  be  the  total  exposed  surfa 
for  each  comoonent  of  the  building.  TDIST  assumes  all  buil  - 

j 

ings  to  be  square  and  symmetrical  and  it  will  use  25%  of  th 
total  specified  areas  of  the  walls  and  windows  when  calcula- 
tions are  performed  on  a oer-wall  basis.  The  ground-contact 
area  of  the  building  should  be  the  total  basement  surface 
area  in  contact  with  the  surrounding  soil  or  the  surface  e, 


of  the  building's  base  slab,  if  no  basement  exists.  The 


roof  area,  if  the  roo 


is  not  flat,  should  be  the  horizontal 


projection  of  the  total  pitched-roof  area.  The  parameter 


values  should  be  typical  for  the  similar  buildings  within 


It  is  recommended  that 


buildings  within  the  community  being  studied 


fied  according  to  type  and  then  that  a representative  build 


ing  be  chosen  as  the  norm  for  each  grouping.  Physical  dimen 


sions  and  construction  materials  specified  for  each  building 


type  should  be  consistent  with  actual  data  for  these  norms 


Care  should  be  taken  to  in 


those  sections  of  the  building  exposed  to  the  external  ele 


ments,  since  the  heat  transfer  calculations  assume  losses 


from  the  building  only  to  the  surrounding  air  and  the  earth 


The  heat  transfer  characteristics  for  compound  materials  must 


be  correctly  combined  to  produce  a single  coefficient  for  the 


particular  section  of  the  building  specified.  Wall  and  roof 


thermal  resistances  should  include  the  inner  surface  film 


convection  resistance  and  the  composite  thermal  resistance 


of  the  materials  from  which  they  are  constructed;  the  code 


computes  the  outer  surface  convection  resistance  based  upon 


a specified  surface  roughness  coefficient  and  the  instantane 


ous  wind  velocity,  and  adds  this  surface  term  to  the  user 


supplied  resistance  to  produce  a total  heat  transfer  coeffi 


cient  for  the  component  being  analyzed.  If  walls  are  con 


structed  differently  throughout  the  building,  the  various 


section  values  should  be  combined  to  produce  a single  average 
wall  coefficient  for  the  building.  When  specifying  these 
parameters,  it  should  be  remembered  that  all  factors  influen- 
cing the  net  transfer  of  heat  from  the  interior  of  the  build- 
ing to  its  surroundings  must  be  included  but  that  heat  trans- 
fer within  the  building  itself  is  neglected. 

Variables : 

I - building  type  code  number 

AW (I)  - building  exposed  wall  area  (in  ft  ) 

RTW(I)  - average  wall  thermal  resistance  (in  hr-ft£-°F 
BTU) 

AWD(I)  - building  window  area  (in  ft“) 

RTWD(I)  - average  window  thermal  resistance  (in  hr-ft" 
°F/BTU ) 

AR(I)  - building  roof  area  (in  ft") 

RTR(I)  - average  roof  thermal  resistance  (in  hr-ft^- 
°F/BTU) 

AF(I)  - building  exposed  floor/basement  area  (in  ft2) 

RTF(I)  - ground-contact  material  thermal  resistance 


(in  br-ft2-°F/BTU) 

Format:  c.c.  1-2  12  I 

c.c.  3-10  F8 . 2 AW( I ) 

c.c.  11-18  F3 . 3 RTW(I) 

c.c.  19-26  F8.2  AW,D  ( I ) 

c.c.  27-3*1  F8 . 3 RTWD(I) 
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c .c . 

35-42 

F3 . 2 

AH  ( I ) 

c . c . 

43-50 

F8.3 

P.TR(I) 

c . c . 

51-58 

FS.2 

AF  ( I ) 

c . c . 

59-66 

F8 . 3 

RTF ( I ) 

l8 . Building  Air  Infiltration  Factor  Card 
This  card  inputs  the  infiltration  air  flow  coeffi- 
cients C and  N and  the  infiltration  crack  dimensions  for  the 
windows  and  doors  of  the  building  type  specified  by  its  code 
number.  Also  input  are  the  coefficients  C and  N for  the 
building's  walls.  One  card  must  be  supplied  for  each  build- 
ing type  as  specified  on  the  Building  Type  Number  Card. 

The  coefficients  C and  N (obtained  from  Table  A. 3)  are 
used  to  translate  the  combined  effects  of  wind  and  stack 
effect  induced  pressure  differences  into  air  flowrates  in 
CFM  per  linear  foot  of  crack  (or,  in  the  case  of  the  walls, 
CFM  per  square  foot  of  surface  area).  Since  the  infiltration 
air  flowrate  for  the  building  is  directly  proportional  to  C, 
extreme  care  should  be  exercised  in  selecting  values  from 
the  extremely  wide  range  of  available  specifications  which 
accurately  reflect  the  representative  building's  construction 
characteristics.  The  value  of  C input  for  the  windows  should 
be  the  sum  of  the  window  and  frame  coefficients;  N should  be 
‘•.he  average  of  the  window  and  frame  coefficients. 

Because  air  generally  flows  into  a building  on  one  side 
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and  out  of  the  building  on  the  opposite  side,  the  crac  : 
lengths  specified  around  doors  and  windows  should  be 
of  the  total  cracks  measured  for  the  building  to  avoid 
double  volume  counting  in  the  flow  calculations. 
Variables  : 

I - building  type  code  number 
DCL(I)  - door  crack  length  (in  ft.) 

DINC(I)  - door  crack  flow  coefficient  C 

DINN(I)  - door  crack  flow  coefficient  N 

WDCL(I)  - window  crack  length  (in  ft.) 

WDINC(I)  - window  crack  flow  coefficient  C 
WDINN(I)  - window  crack  flow  coefficient  N 
WINC(I)  - wall  crack  flow  coefficient  C 

WINN(I)  - wall  crack  flow  coefficient  N 

Format:  c.c.  1-2  12  I 


c.c. 

3-10 

F8 . 3 

DCL(I) 

c.c. 

11-18 

F8 . 3 

DINC(I) 

c.c. 

19-26 

F8 . 3 

DINN(I) 

c.c. 

27-34 

F8 . 3 

WDCL(I) 

c.c. 

35-42 

F8 . 3 

WDINC(I) 

c.c. 

43-50 

F8 . 3 

WDINN(I) 

c.c. 

51-58 

F8 . 3 

WINC(I) 

c.c. 

59-66 

F8 . 3 

WINN  ( I ) 
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19 . Building  Miscellaneous  Data  Card 

This  card  inputs  codes  for  the  surface  materials 
of  the  walls  and  roof,  the  height,  and  the  nominal  maximum 
ventilation  air  flowrate  for  the  building  specified  by  its 
type  code  number.  The  building  height  must  be  specified  in 
ft.  and  the  air  flowrate  in  CFM.  One  card  must  be  input  for 
each  building  type  as  specified  on  the  Building  Type  Number 
Card . 

The  wall  and  roof  surface  material  type  codes  are  used 
as  indices  for  internally-stored  surface  roughness  coeffi- 
cients used  in  the  calculation  of  the  surface  convection 
thermal  resistances.  The  codes  should  be  input  as  integer 
values  according  to  the  following  classifications: 

1 - scucco  4 - clear  pine 

2 - brick,  rough  plaster  5 - smooth  plaster 

3 - concrete  6 - glass,  white  paint 

on  pine 

When  selecting  the  appropriate  code  for  the  wall  or  roof 
material  in  the  building  being  specified,  it  should  be  remem- 
bered that  the  important  parameter  indicated  by  these  codes 
is  the  relative  surface  roughness  and  not  the  exact  type  of 
surface  material. 

The  peak  ventilation  air  flowrate  is  used  as  an  index 
for  the  scheduling  of  ventilation  whicK  is  compared  with  the 
computed  infiltration  air  flowrates  to  obtain  the  heat  loss 
due  to  direct  air  exchanges  between  the  building's  interior 
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and  exterior.  This  flowrate  need  not  be  limited  to  only 
buildings  containing  forced-air  ventilation  equipment,  since 
the  code  considers  the  opening  of  windows  to  cause  a venti- 
lation rather  than  an  infiltration  air  exchange.  The  peak 
value  input  on  this  card  should  generally  be  obtained  from 
either  building  code  specifications  of  the  required  ventila- 
tion rates  for  various  structural  and  occupancy  classifica- 
tions, or  it  should  be  based  upon  measured  air  exchange  data 
(such  as  a typical  number  of  air  changer  per  hour)  for  the 
specified  building. 

Variables : 

I - building  type  code  number 
MTW(I)  - wall  surface  material  type  code 
MTR(I)  - roof  surface  material  type  code 
BH(I)  - building  height  (in  ft.) 

AINFL(I)  - peak  ventilation  air  flowrate  (in  CFM) 


Format:  c.c 

. 1-2 

12 

I 

c . c 

. 5-6 

12 

MTW(I) 

c . c 

. 9-10 

12 

MTR(I) 

c.c 

. 13-19 

F7.2 

BH(I ) 

c . c 

. 20-29 

F10.2 

AINFL(I) 

20.  Building  Solar  Data,  Internal  Heating,  and 

Dome 

Stic  Hot 

Water  Data 

Card 

This  card 

inputs 

the  building 

orientation,  internal 

window  shading  coefficient,  wall  and  roof  shading  coeffi- 
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cienas,  wall  and  roof  material  solar  absorptivities , the 
nominal  maximum,  rate  of  non-space-conditioning  heat  produc- 
tion, and  the  nominal  maximum  rate  of  service  hot  water 
consumption  for  the  building  specified  by  its  type  number 
code.  One  card  is  required  for  each  building  type  as  speci- 
fied on  the  Building  Type  Number  Card. 

The  orientation  of  the  building  is  measured  in  degrees 
clockwise  from  North,  and  is  referenced  to  that  wall  whose 
normal  is  closest  to  North.  Thus,  if  a building  has  one  wa . 1 
facing  north,  its  orientation  angle  would  oe  0°;  one  wall 
facing  Northeast  would  result  in  an  orientation  of  4o°.  If 
the  specified  building  represents  an  aggregation  of  several 
randomly-oriented  units,  the  reference  wall  angle  should 
approximate  any  observed  trends  in  the  actual  orientations. 
(It  should  be  noted  that  the  orientation  must  always  lie  with 
in  the  range  0°  £ WA < 90°,  since  the  building  is  assumed  to 
be  square).  For  the  purposes  of  computing  loads  on  a per- 
wall  basis,  the  reference  wall  is  designated  as  #1,  and  the 
walls  are  numbered  clockwise  from  the  North. 

The  internal  window  shading  coefficient  input  on  this 
card  is,  in  a sense,  a "reciprocal  shading  coefficient”  since 
it  actually  indicates  a measure  of  the  amount  of  solar  energy- 
incident  upon  a window  which  is  transmitted  and  converted  to 
heat  within  the  building.  The  coefficient  is  input  as  a 
decimal  fraction,  with.  0 indicating  full  shading  and  1.0 
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indicating  no  shading.  Included  in  this  value  should  be  the 
effects  of  such  fixtures  as  double  glazing,  tinted  glass, 
Venetian  blinds,  drapes,  shades,  etc.  (i.e.,  anything  in- 
side the  exposed  surface  of  the  window  which  would  modify 
the  solar  heating  obtained  from  clear  1/8"  single  pane  glass). 
Table  A. 6 contains  shading  coefficients  for  a range  of 'typi- 
cal window  configurations,  and  References  4 and  5 contain 
more  complete  listings. 

Exterior  surface  shading  coefficients  (actually,  reci- 
procal shading  coefficients  or  lighted  fractions)  are  input 
on  this  card  for  each  of  the  four  walls  of  the  building  and 
its  roof.  As  is  discussed  above,  the  walls  are  numbered 
consecutively  from  1 to  4,  clockwise  from  the  North,  and  the 
roof  Is  assigned  an  equivalent  wall  index  of  5.  The  exterior 
shading  coefficient  for  each  of  these  surfaces  should  provide 
an  estimate  of  the  fraction  of  that  surface,  under  direct 
sunlight  conditions,  which  is  actually  lit  by  solar  radiation. 
Thus,  if  the  building  is  immediately  adjacent  to  a larger 
building  abutting  its  southern  wall,  the  coefficient  for 
wall  #3  would  be  0;  conversely,  if  the  North  wall  is  complete- 
ly free  of  shading,  even  though  the  sun  may  never  actually 
shine  directly  upon  it  due  to  its  position  in  the  sky,  the 
coefficient  for  that  wall  should  be  specified  as  1.0  to  indi- 
cate that  any  radiation  available  to  the  wall  will  reach  it 
unimpeded.  When  aggregating  several  randomly-oriented  build- 
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ings  into  a single  representative  unit,  estimates  should  be 
made  of  the  net  effects  of  exterior  shading  upon  these 
buildings  and  the  aggregate  data  should  be  used  to  derive 
effective  wall  and  roof  shading  coefficients. 

The  wall  and  roof  material  solar  absorptivities  are 
measures  of  the  fraction  of  the  solar  radiation  incident  upon 
those  surfaces  which  is  absorbed  and  converted  to  heat  within 
the  building.  In  general,  the  values  specified  will  depend 
upon  both  the  color  and  the  composition  of  the  surface.  A 
range  of  typical  absorptivities  is  listed  in  Table  A. 7,  which 
has  been  extracted  from  more  comprehensive  summaries  in 
References  4 and  5 • 

Also  input  on  this  card  is  a value  for  the  nominal  maxi- 
mum rate  of  internal  building  heat  generation  from  all  sources 
not  directly  associated  with  space  conditioning  (e.g.,  lights, 
motors,  appliances,  people).  This  peak  value  should  be  input 
in  kW,  and  it  serves  as  an  index  for  the  scheduling  of  non- 
space-conditioning heat  generation  as  a function  of  building 
occupancy  and  intensity  of  use. 

Finally,  a measure  of  the  nominal  maximum  rate  of  service 

b . , 

hot  water  usage  (e.g.,  showers,  wash  basins,  cooking,  etc.)  is 
required  to  be  input  in  units  of  BTU/hr.  This  peak,  as  are 
those  for  ventilation  and  internal  heat  generation,  is  used 
as  an  index  for  an  hourly  hot  water  heating  demand  schedule 
for  the  building.  Since  most  tabulations  and  field  measure- 
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ments of  domestic  hot  water  consumption  are  presented  in 
units  of  gallons/hour,  converting  these  rates  to  the  data 
value  required  on  this  card  necessitates  estimates  to  be 
made  of  both  the  municipal  water  supply  temperature  to  the 
building  and  of  the  water  temperature  at  the  heater  outlet 
to  determine  the  energy  demand  per  unit  flow. 

Variables : 

I - building  type  code  number 

WA ( I ) - building  orientation  (in  degrees  from  North) 
WDSC(I)  - window  internal  shading  coefficient  (decimal 
fraction  of  incident  radiation  which  is 
transmitted  as  heat) 

FWLIT(I,NN)  - fraction  of  surface  NN  open  to  direct 
30lar  radiation  [NN=l-4:  walls, 

NN  = 5 : roof] 

WSAB(I)  - wall  material  solar  absorptivity  (decimal 
fraction  of  incident  radiation  absorbed) 
RSAB(I)  - roof  material  solar  absorptivity  (decimal 
fraction  of  incident  radiation  absorbed) 
POL(I)  - nominal  peak  internal  heat  generation  rate 
from  non-space-conditioning  sources  (in  kW) 
DHW(I)  - nominal  peak  domestic  hot  water  consumption 
rate  (in  BTU/hr) 

* 

! 
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Format : 


c . c . 

1-2 

12 

I 

c . c . 

3-9 

F7.2 

WA  ( I ) 

c . c . 

10-14 

F5.2 

WDSC(I) 

c . c . 

15-19 

F5.2 

FWLIT( I , 1 ) 

c . c . 

20-24 

F5.2 

FWLIT(I,2) 

c . c . 

25-29 

F5.2 

FWLIT(I ,3) 

c . c . 

30-34 

F5.2 

FWLIT(I,4) 

c . c . 

35-39 

F5.2 

FWLIT(I , 5 ) 

c . c . 

40-44 

F5.2 

WSAB(I) 

c . c . 

45-49 

F5.2 

RSAB(I) 

c . c . 

50-58 

F9.2 

POL(I) 

c .c . 

59-69 

Ell. 4 

DHW(I) 

21 . Building  External  Supply  Data  Card 
This  card  inputs  the  number  of  units  of  a building 
type  which  are  not  served  by  the  TUS  and,  of  that  number. 


the  fraction  which  have  electric  resistance  heaters  and  the 
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fraction  having  heat  pumps.  One  card  must  be  supplied  for 
each  of  the  building  types  specified  on  the  Building  Type 
Number  Card . 

Those  buildings  not  supplied  with  HTW  are  considered  to 

r * 

be  part  of  the  TES  "external  system"  and  have  their  energy 
needs  supplied  electrically.  Because  no  detailed  electrical 
distribution  system  analysis  is  performed  by  the  code,  only 
the  total  number  of  buildings  and  their  mode  of  energy 
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supply  are  required  for  the  calculation  of  the  external  sys- 
tem electrical  demand.  It  is  assumed  that  all  the  specified 
external  system  buildings  have  compressive  air  conditioners. 
However,  the  required  resistance  and  heat  pump  fractions 
need  not  sum  to  1.0,  since  the  code  allows  a fraction  of  the 
external  buildings  to  be  heated  from  in-house  fossil-fired 
furnaces  which  present  no  load  to  the  TES. 

Variables : 

I - building  type  code  number 

NUEX(I)  - number  of  units  not  connected  to  TUS 
FREXT(I)  - fraction  of  electrically-supplied  units 

having  resistance  heat  (decimal  fraction) 
FHPEX(I)  - fraction  of  electrically-supplied  units 
having  heat  pumps  (decimal  fraction) 


Format : 

c.c.  1-2 

12 

I 

c.c.  3-7 

15 

NUEX(I) 

c.c.  8-12 

F5.2 

FREXT(I) 

c.c.  13-17 

F5.2 

FHPEX(I) 

22 . Building  Type  and  End-Use  Equipment  Distribu- 
tion Card 

This  card  inputs  the  number  of  buildings  of  a given 
type  served  by  a load  center  heat  exchanger.  Also  given  are 
the  total  number  of  compressive  and  absorptive  air  condi- 
tioning units  and  the  number  of  electrical  resistance  heating 
units,  heat  pumps  and  hot  water  heating  units  contained  in 
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these  buildings.  Since  the  program  reads  these  cards 

according  to  building  type  within  each  heat  exchanger  area,  , 

one  card  must  be  supplied  for  each  building  type  for  each 
load  center  heat  exchanger,  even  if  a particular  heat  ex- 
changer  does  not  supply  every  type  (e.g.,  if  there  are  4 
secondary  loops  with  5 heat  exchangers  per  loop  and  a total 

_ 

of  6 buildings  types  in  the  community,  a total  of  4x5x6  = 120 
Building  Type  and  End-Use  Equipment  Distribution  Cards  must 
be  supplied).  The  order  in  which  these  cards  should  be 
placed  in  the  data  deck  is:  for  the  first  load  center  heat 

exchanger,  NTYPE  cards  should  be  input  consecutively  accord- 
ing to  the  building  code  numbers  defined  by  the  Building 
Type  Number  Card;  this  sequence  is  then  repeated  in  succes- 
sion for  every  load  center  heat  exchanger. 

Buildings  should  be  apportioned  to  the  various  load 
center  heat  exchangers  so  as  to  provide  a fairly  uniform 
load  distribution  among  the  heat  exchangers  in  a given  loop 
and  among  all  the  loops  in  the  system.  If  no  buildings  of  a 
given  type  are  served  by  a given  heat  exchanger,  "0"  may 
bq  specified  for  the  number  of  buildings  on  the  Building 
Type  and  End-Use  Equipment  Distribution  Card,  or  the  card 
may  be  left  blank. 

End-Use  equipment  should  be  distributed  throughout  the 
system  so  as  to  provide  an  optimum  thermal/electric  energy 
utilization  efficiency  for  the  community  as  a whole.  The 


J 


147. 


data  given  on  this  card  for  the  numbers  of  individual  heat- 
ing and  cooling  units  serving  a particular  building  type  in 
a given  area  of  the  community  is  used  by  the  program  to 
distribute  the  space  conditioning  loads  among  the  various 
equipment  categories  provided.  The  specific  numbers  given 
can  be  the  number  of  buildings  served  by  a particular  equip- 
ment type  or  the  total  number  of  individual  units  of  that 
type  in  the  area  being  specified.  However,  the  numbers 
should  be  specified  consistently  throughout  the  system  and 
must  present  a true  indication  of  the  relative  equipment 
distribution  in  each  consumer  sector. 

Variables : 

NU(IX,I)  - number  of  building  units  of  type  I served 
by  heat  exchanger  IX 

NCMP(IX,I)  - number  of  compressive  air  conditioning 
units  in  buildings  of  type  I for  heat 
exchanger  IX 

NABS(IX,I)  - number  of  absorptive  air  conditioning 
units  in  buildings  of  type  I for  heat 
exchanger  IX 

NRES(IX,I)  - number  of  electrical  resistance  heating 
units  in  buildings  of  type  I for  heat 
exchanger  IX 

NHPM(IX,I)  - number  of  heat  pumps  in  buildings  of  type 


I for  heat  exchanger  IX 
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NHWS(IX,I)  - number  of  hot  water  heating  units  in 

buildings  of  type  I for  heat  exchanger  IX 


Format : 

c . c . 

1-5 

15 

NU ( IX , I ) 

c . c . 

6-10 

15 

NCMP ( IX , I ) 

c . c , 

11-15 

15 

NABS (IX, I) 

c . c . 

16-20 

15 

NRES (IX, I ) 

c . c . 

21-25 

15 

NHPM(IX,I) 

c . c . 

26-30 

15 

NHWS ( IX , I ) 

23. 

End- 

Use  Equipme 

nt  Efficiency  Data  Card 

This 

card 

inputs  the 

eff 

iciencies  (or  coefficients 

of  performance) 

of 

the  compressive 

and  absorptive  air  condi- 

tioning  units. 

the 

electrical 

resistance  heating  units,  heat 

pumps,  and  hot  water  heating  units  in  the  system.  For  the 
compressive  air  conditioning  units,  heat  pumps  and  electrical 
resistance  heating  units,  this  efficiency  should  be  expressed 
in  terms  of  the  rated  BTU/hr  heating  or  cooling  capacity  per 
BUT/hr  equivalent  of  the  required  input  electrical  power. 

For  the  absorptive  air  conditioning  and  hot  water  heating 
units,  the  efficiency  should  be  expressed  in  terms  of  the 
rated  BTU/hr  heating  or  cooling  capacity  per  BTU/hr  of  the 
input  thermal  power.  The  specified  efficiencies  should  be 
consistent  with  the  operating  characteristics  of  the  various 
equipment  under  system  design  supply  temperature  and  room 
space  conditioning  requirements.  The  nominal  heat  pump  COP 
should  be  specified  at  an  ambient  air  temperature  of  32  °F 
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and  a room  temperature  of  68  °F;  nominal  compressive  and 
absorptive  air  conditioning  COPs  are  defined  at  an  air 
temperature  of  90  °P  and  a room  temperature  of  75  °F.  A 
non-zero  value  for  the  efficiency  of  each  equipment  type 
should  be  given,  even  though  a given  type  may  not  appear 
anywhere  in  the  system. 

Variables : 

ECMP  - compressive  air  conditioning  unit  coefficient 
of  performance  (dimensionless) 

EABS  - absorptive  air  conditioning  unit  coefficient 
of  performance  (dimensionless) 

ERES  - electrical  resistance  heating  unit  efficiency 
(dimensionless ) 

EHPM  - heat  pump  coefficient  of  performance  (dimen- 
sionless) 

EHWS  - hot  water  heating  unit  efficiency  (dimensionless) 
Format:  c.c.  1-5  F5.2  ECMP 

c.c.  6-10  F5.2  EABS 

c.c.  11-15  F5-2  ERES 

c.c.  16-20  F5.2  EHPM 

c.c.  21-25  F5.2  EHWS 


the  Nodal  Fluid  Parameter  Information  Cards.  In  general,  a 
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value  for  this  factor  other  than  1.0  should  not  be  specified 
unless  the  only  input  parameters  to  be  varied  between  two 
given  simulations  are  the  nodal  mass  flowrates,  which  are 
scaled  uniformly  throughout  the  TUS.  In  this  instance, 
this  factor  should  be  the  ratio  between  the  desired  flowrates 
and  the  originally  specified  values,  and  the  Nodal  Fluid 
Parameter  Information  Cards  need  not  be  re-punched. 

Variables: 

CSFAC  - case  input  mass  flowrate  scale  factor 
Format:  c.c.  1-7  F7.4  CSFAC 

25 • Solar  Data  Card 

On  this  card  are  input  the  latitude  of  the  communi- 
ty being  studied,  and  the  solar  declination  angle,  the  appar- 
ent solar  constant,  the  atmospheric  extinction  coefficient 
and  the  sky  diffuse  factor  for  the  day  on  which  the  simula- 
tion is  performed.  Also  input  is  the  hour  of  the  day  at 
which  the  initial  system  loads  are  to  be  calculated  and  the 
steady-state  system  analysis  is  to  be  performed. 

The  latitude  should  be  input  in  (positive)  degrees 
north  or  south  of  the  equator. 

The  solar  declination  angle,  apparent  solar  constant, 
atmospheric  extinction  coefficient  and  sky  diffuse  factor  are 
observed  insolation  parameters  which  depend  upon  the  time  of 
year  as  summarized  In  Table  A. 4.  The  solar  declination  angle 
should  be  input  in  degrees.  The  apparent  solar  constant. 
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which  measures  the  intensity  of  a direct  solar  beam  normally 
incident  upon  a horizontal  surface  under  clear  skies,  is 

p 

input  in  BTU/hr-ft  . The  atmospheric  extinction  coefficient 
indicates  the  absorption  in  the  atmosphere  of  direct  solar 
radiation,  and  the  sky  diffuse  factor  measures  the  disper- 
sion of  the  solar  beam  as  it  travels  through  the  atmosphere. 

The  time  of  day  input  on  this  card  is  the  hour  at  which 
the  initial  system  calculations  are  to  be  performed  and  is 
the  beginning  time  of  the  dynamic  analysis.  Since  all  time- 
dependent  data  is  read  sequentially  on  an  hourly  basis  with- 
out specifying  the  precise  time  of  day  with  which  it  is 
associated,  this  beginning  hour  indicator  initializes  the 
time  clock  and  is  used  to  assign  solar  time  values  to  each 
hourly  data  point.  A continuous  24-hour  time  clock  is 
employed  by  the  model,  with  midnight  being  assigned  an  hour 
indicator  of  0 and  noon  an  indicator  of  12;  a simulation 
covering  the  12-hour  period  between  2 PM  and  2 AM,  for  exam- 
ple, would  have  14  as  its  specified  beginning  hour  and  would 
terminate  at  hour  2. 

Variables : 

ALAT  - latitude  of  location  (in  degrees) 

DECL  - solar  declination  angle  (in  degrees) 

ASC  - apparent  solar  constant  (in  BTU/hr-ft^) 

AEC  - atmospheric  extinction  coefficient 
SDF  - sky  diffuse  factor 


152. 


♦ 


Format : c . c . 

1-5 

F5.2 

ALAT 

c . c . 

6-11 

F6.2 

DECL 

c . c . 

12-16 

F5-0 

ASC 

c . c . 

17-22 

F6 . 3 

AEC 

c . c . 

23-28 

F6 . 3 

SDF 

c . c . 

29-31 

13 

IBHR 

26 . Power  Plant  Rating  and  Peak  Non-Space-Condi- 
tioning Electrical  Demand  Data  Card 
This  card  inputs  the  power  plant  thermal  rating 
and  the  peak  non-space-cond^tioning  electrical  load  for  the 
community.  The  power  plant  output  must  be  given  in  MWt  and 
the  electrical  load  in  MWe. 

The  specified  power  plant  thermal  output  is  held  con- 
stant throughout  the  simulation  and  provides  the  total  power 
input  to  both  the  electrical  and  thermal  energy  distribution 
systems.  This  output  power  level  should  be  set  approximately 
equal  to  the  average  total  energy  demand  for  the  community 
over  the  period  of  the  simulation.  Higher  outputs  will 
result  in  a net  storage  of  energy  in  the  reservoir  during 
the  simulation;  lower  outputs  will  result  in  a net  removal 
of  heat  from  the  reservoir  and  possible  inability  to  meet 
the  peak  energy  demands  of  the  community.  The  program  uses 
the  power  plant  output  as  primary  input  to  the  dynamic  system 
calculations;  the  heat  available  to  the  thermal  distribution 
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system  is  given  by  the  difference  between  this  constant  out- 
put and  the  instantaneous  energy  demand  for  electrical  power 
production. 

The  electrical  load  specified  on  this  card  should  be 
the  peak  electrical  power  demand,  excluding  space  condi- 
tioning, for  the  community  over  the  period  of  the  simulation. 
To  eliminate  space  conditioning  effects,  the  electrical  load 
data  specified  should  be  representative  of  a minimum  demand 
period  during  the  Autumn  or  Spring,  when  neither  air  condi- 
tioning nor  space  heating  is  required. 

Variables : 

RTBL  - power  plant  constant  total  power  output  (in  MWt) 
ELBL  - peak  community  electrical  load  excluding  space 
conditioning  (in  MWe) 


Format : 

c.c.  1-10 

F10.0 

RTBL 

c.c.  11-20 

F10. 0 

EL3L 

27 . Primary  Program  Flow  Control  Code  Card 
This  card  inputs  the  primary  program  flow  control 
code  KODE,  which  identifies  the  general  type  of  program  exe- 
cution to  be  performed  according  the  the  following  conventions 
KODE=0:  Perform  building  load  analysis  (if  IDES=1) 

or  perform  TUS  steady-state  convergence 
(if  IDESs0 ) , 

K0DE=1:  Perform  hourly  building  lead  summary  (if 

sequence  KODE=0,  TDES=1  input  previously' 
or  perform  dynamic  system  analysis,  and 
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K0DE=2:  Terminate  program  execution. 

In  limited  cases,  KODE  may  be  input  more  than  once  for  a 
given  computer  run  to  transfer  control  of  program  execution 
from  one  program  option  to  another.  Examples  of  such  sequen- 
tial operations  are: 

KODE=0,  IDES=1 , K0DE=1:  Computes,  adds  and  distributes 

building  loads  at  one  time  instant  and  computes  up 
to  24  hourly  building  space  conditioning  loads, 

KODE=0 , IDES=0,  K0DE=1:  Steady-state  system  analysis 

used  to  internally  initialize  dynamic  simulation 
TUS  parameters,  and 

Nesting  of  several  combinations  of  KODE=0  and  IDES=Q,1 
are  also  possible  to  analyze  several  instantaneous 
system  loading  conditions  within  the  same  TUS 
configuration  without  performing  any  dynamic  energy 
system  analyses. 

In  general,  it  is  recommended  that  program  execution  be  ter- 
minated following  the  completion  of  a dynamic  system  simula- 
tion, since  the  retention  of  TUS  parameter  values  from  the 
dynamic  calculations  could  invalidate  certain  portions  of 
directly  continued  steady-state  or  building  load  analyses. 

Variables : 

KODE  - primary  program  flow  control  code  (0,1,2) 

Format:  c.c.  1-3  13 
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28 . Building  Load  Analysis  Control  Code  Card 
The  building  load  analysis  control  code  IDES  must 
be  input  immediately  following  each  specification  of  KODE=0. 

As  summarized  in  the  description  of  the  Primary  Program  Flow 
Control  Code  Card  and  in  Table  A.9S 

* 

IDES=0:  Perform  steady-state  TUS  convergence  and 

analysis,  or 

IDES=1:  Perform  building  load  analysis  and  print  net 

building  space  conditioning  demands;  aggre- 
gate and  print  initial  building  loads  and 
distribution,  and  print  initial  thermal  demands 
at  each  load  center  heat  exchanger. 

IDES  should  not  be  specified  following  the  input  of  K0DE=1 
unless  the  system  is  re-initialized  by  reading  KODE=0. 

Variables : 

IDES  - building  load  analysis  control  code  (0,1) 

Format:  c.c.  1-3  13  IDES 

29 • Initial  Air  Temperature  Data  Card 
This  card  inputs  the  ambient  air  temperature  to  be 
used  in  determining  the  initial  building  space  conditioning 
energy  demands.  If  a dynamic  simulation  is  to  be  performed 
by  using  the  steady-state  system  analysis  to  initialize  the 
TUS  parameters,  this  temperature  must  be  the  initial  hourly 
data  point  for  the  period  being  studied  (i.e.,  the  temperature 
at  time  t=IBHR).  The  temperature  must  be  specified  in  °F . 


1 


I 


JL 


TAIR( 25 ) - air  temperature  at  time  t=IBHR  (in  °F) 

(the  index  25  is  reserved  by  the  program 
for  initial  conditions  data) 

Format:  c.c.  1-5  F5.2  TAIR(25) 

30 . Initial  Wind  Velocity  Data  Card 

Input  on  this  card  is  the  initial  wind  velocity, 
in  mph,  for  the  location  at  time  t = IBHR. 

Variables : 

WV('25)  - wind  velocity  at  time  t = IBHR  (in  mnh) 
Format:  c.c.  1-5  F5 . 2 WV(25) 

31 . Initial  Wind  Direction  Data  Card 

On  this  card  is  input  the  direction  from  which  the 
wind  is  blowing  at  time  t = IBHR.  The  direction  is  measured 
in  degrees,  clockwise  from  North  (e.g.,  a wind  incident  from 
the  Southwest  would  have  a direction  angle  of  225°). 
Variables : 

WDIR(25)  - wind  incidence  direction  at  time  t = IBHR 
„ (in  degrees  from  North) 

Format : c.c.  1-5  F5.2  WDIR(25) 

32 . Initial  Electric  Load  Factor  Data  Card 
This  card  inputs  the  electric  load  factor  which 

multiplies  the  peak  community  demand  specified  on  the  Power 
Plant  Rating  and  Peak  Non-Space-Conditioning  Electrical 


Demand  Data  Card  to  determine  the  non-space  conditioning 
electric  load  during  the  steady-state  calculations.  If  a 
dynamic  simulation  is  to  be  run  using  the  steady-state  sub- 
routines to  initialize  the  system,  this  load  factor  must  be 
the  initial  hourly  data  point  (i.e.,  the  factor  at  time 
t = IBHR)  for  the  period  being  studied.  The  factor  is  dimen- 
sionless and  merely  represents  the  ratio  of  the  instantaneous 
system  electrical  load  to  the  peak  load  for  the  period 
under  investigation.  (Inputting  such  a combination  of  a 
peak  load  and  load  factors  allows  the  simple  retention  of  a 
given  load  profile  while  varying  the  size  of  the  community 
to  be  supplied) . 

Variables : 

ELF(25)  - electric  load  factor  (ratio  of  instantaneous 
to  peak  non-space-conditioning  electric 
load)  at  time  t = IBHR 
Format:  c.c.  1-5  F5.2  ELF (25) 

33 . Initial  Building  Room  Temperature  Data  Card 
This  card  inputs  the  room  temperatures,  by  building 
type,  to  be  used  in  determining  the  building  space  condition- 
ing loads  for  the  steady-state  system  calculations.  If  a 
dynamic  simulation  is  to  be  run  using  the  steady-state  sub- 
routines to  initialize  the  system  parameters,  these  tempera- 
tures must  be  the  initial  hourly  data  points  (i.e.,  the 
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temperatures  at  time  t = IBHR)  for  the  period  being  studied. 
Each  building  type  is  assumed  to  have  a characteristic  room 
temperature  throughout  all  the  rooms  of  that  building  type 
in  the  community.  These  temperatures  are  allowed  to  vary 
throughout  the  simulation  to  correspond  to  the  changing 
requirements  of  the  people  using  the  buildings.  All  tempera- 
tures must  be  specified  in  °F.  One  temperature  3s  required 
for  each  of  the  building  types  specified  on  the  Building 
Type  Number  Card.  A maximum  of  20  building  types  is  allowed. 
If  more  than  10  building  types  are  specified,  two  cards  must 
be  used  to  input  the  temperatures. 

Variables : 

TSB(25,I)  - desired  room  temperature  for  building  type 
Format : 


I at  time 

t = IBHF 

(it 

z °F)  (maximum 

(Card  1) 

c.c.  1-7 

F7.2 

TSB( 25 , 1 ) 

c.c.  8-14 

F7.2 

TSB( 25 , 2 ) 

c.c.  64-70 

F7.2 

TSB( 25,10) 

(Card  2)  (not 

required 

if 

NTYPE  < 10) 

c.c.  1-7 

F7.2 

TSB (25,11) 

c.c.  8-14 

F7.2 

TSB (25,12) 

F7.2 


TSB(25,NTYPE) 
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34 . Initial  Building  Use  Factor  Data  Card 
On  this  card  are  Input  the  occupancy-related 
building  use  factors  which  multiply  the  peak  ventilation  air 
flowrate  (from  the  Building  Miscellaneous  Data  Cards)  and 
the  peak  non-space-conditioning  equipment  heat  generation 
rate  (from  the  Building  Solar  Data,  Internal  Heating,  and 
Domestic  Hot  Water  Data  Cards)  to  obtain  the  instantaneous 
values  of  these  parameters  for  each  building  type  at  time 
t = IBHR.  It  is  assumed  by  the  code  that,  in  general,  venti- 
lation requirements  and  internal  heat  generation  are  func- 
tions of  building  occupancy  and  intensity  of  usage.  There- 
fore, these  use  factors  provide  a measure  of  the  relative 
occupancy  of  each  building  type  and  should  be  given  as  the 
ratio  of  the  instantaneous  ventilation  or  internal  heating 
rates  to  the  peaks  observed  during  the  simulation  period. 

If  more  than  10  building  types  are  specified,  two  cards 
must  be  used  to  input  the  factors. 

Variables : 

BUF(25,I)  - use  factor  for  building  type  I at  time 
t = IBHR  (decimal  fraction  of  peak  use) 
Format:  (Card  1) 


c.c.  1-7 

F7.2 

BUF (25,1) 

c.c.  8-14 

F7.2 

BUF (25,2) 

c.c.  64-70 

F7.2 

BUF (25, 10) 
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(card 

2)  (not 

requi 

red 

if  NTYPE  < 

10) 

c . c . 

1-7 

F7 

.2 

BUF (25 ; 

,11) 

c . c . 

8-114 

F7 

.2 

BUF (25, 

,12) 

F7.2  3UF( 25 ,NTYPE) 

3 5 • Initial  Domestic  Hot  Water  Use  Factor  Da":  Carl 
This  card  inputs  the  use  factors  which  multiply  the 
peak  domestic  hot  water  demands  from  the  Building  Solar  Data, 
Internal  Heating,  and  Domestic  Hot  Water  Data  Cards  to  obtain 
the  instantaneous  energy  demands  for  domestic  hot  water 
service  at  time  t = IBHR.  Although  it  is  correlated  with 
building  occupancy,  hot  water  consumption  generally  follows 
a different  schedule  from  that  specified  by  the  general  build- 
ing use  factor  and  is  therefore  allowed  to  vary  independently 
from  the  building  ventilation  and  internal  heat  generation 
rates  determined  by  BUF.  If  a building  has  no  domestic  hot 
water  consumption  (e.g.,  a warehouse  or  small  commercial 
building),  its  use  factor  should  simply  be  input  as  0 or 
left  blank  on  this  card.  If  more  than  10  building  types  are 
specified,  two  cards  must  be  used  to  input  the  factors. 

Variables : 

WUF(25,I)  - domestic  hot  water  use  factor  for  building 


type  I at  time  t = IBHR  (decimal  fraction 
of  peak  use) 


Format : 
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(card  1) 


c . c . 

1-7 

F7.2 

WUF (25,1) 

c . c . 

8-14 

F7 .2 

WUF (25,2) 

c . c . 

64-70 

F7.2 

WUF (25,10) 

(card 

2)  (not 

required 

if  NTYPE  < 10) 

c . c . 

1-7 

F7.2 

WUF  (25,  ID 

c . c . 

^r 
1 ) 

\ 

CO 

F7.2 

WUF (25, 12) 

• 

F7.2 

WUF (25, NTYPE) 

3 6 . Hourly  Data  Point  Number  Card 
This  card  inputs  the  number  of  hourly  data  points 
to  be  used  in  the  dynamic  system  simulation.  The  first 
hourly  data  point  is  assumed  to  be  given  at  time  t = IBHR+1. 
Therefore,  the  number  specified  on  this  card  must  not  include 
the  initial  data  values  used  by  the  steady-state  initializa- 
tion subroutines,  and  it  should  correspond  to  the  total  dura- 
tion of  the  simulation  in  hours.  For  example,  a 6-hour  simu- 
lation would  require  6 data  points  in  addition  to  the  initial- 
ization values,  and  "6"  would  be  specified  on  this  card. 

Since  the  maximum  single  simulation  period  is  24  hours,  a 
maximum  of  24  hourly  data  points  per  simulation  is  allowed. 

Variables : 

NTEMP  - number  of  hourly  data  points  in  simulation 
(maximum  of  24) 
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Format:  c.c.  1-3  13  NTEMP 

3 7 . Hourly  Air  Temperature  Data  Card 
This  card  inputs  the  hourly  external  air  tempera- 
ture values  to  be  used  in  determining  the  building  space 
conditioning  loads  during  the  dynamic  system  simulation. 

The  first  value  specified  should  be  the  air  temperature  at 
time  t = IBHR+1  (i.e.,  one  hour  into  the  simulation).  One 
temperature  value  must  be  given  for  each  of  the  data  points 
specified  on  the  Hourly  Data  Point  Number  Card.  The  program 
uses  a linear  interpolation  routine  to  determine  intermediate 
temperatures  for  each  simulation  time  step.  Therefore,  the 
specified  temperatures  should  accurately  reflect  the  actual 
air  temperature  profile  for  the  time  period  being  studied. 

The  temperatures  should  be  listed  consecutively  on  this  card 
according  to  the  given  input  format;  a maximum  of  12  values 
per  card  and  a maximum  of  2 cards  per  simulation  (i.e.,  a 
maximum  of  24  hourly  data  points)  is  allowed.  The  temperatures 
must  be  given  in  units  of  °F. 

Variables : 

TAIR(I)  - external  air  temperature  at  time  I hours 
into  the  simulation  (maximum  of  I = 24 
values ) . (in  °F) 


1 


m. 


i 
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Format : 

(card  1) 

c . c . 1-5 

F5.2 

TAIR(l) 

c.c.  6-10 
« 

F5.2 

TAIR ( 2 ) 
» 

• 

C.C.  56-60 

F5.2 

TAIR( 12 ) 

(card  2)  (not  required 

if  NTEMP  < 12) 

c.c.  1-5 

F5.2 

TAIR ( 13 ) 

c.c.  6-10 

F5.2 

TAIR(14) 

• 

F5.2 

TAIR (NTEMP) 

(Max  of  NTEMP 

38. 

Hourly  Wind 

Velocity 

Data  Card 

This  card  inputs  the  hourly  wind  velocities  for 
the  simulation  period.  The  first  value  specified  is  assumed 
' to  obtain  at  time  t = IBHR+1  (i.e.,  one  hour  into  the  simula- 

tion), and  the  total  number  of  data  points  must  correspond 
to  the  number  of  simulation  hours  specified  on  the  Hourly 
Data  Point  Number  Card.  Velocities  should  be  in  units  of 
miles/hr.  If  the  simulation  is  to  span  more  than  12  hours, 
two  data  cards  are  required  as  shown  below. 

Variables : 

WV ( I ) - wind  velocity  at  time  I hours  into  the  simula- 
tion (maximum  of  1=24  values)  (in  mph) 

kJ. _ . . _ _ . __ 


(card  1) 


For"'.',  : 


c.c.  1-5 

?5.2 

WV(1) 

c.c.  6-10 

F5-2 

V.'  1(2) 

c.c.  56-60 

F5.2 

WV ( 1 2 ) 

(card  2)  (not 

required  if 

NTEMP  < 12) 

c.c.  1-5 

F5.2 

WV ( 13 ) 

c.c.  6-10 

F5.2 

WV(llJ) 

F5.2  WV(NTEMP) 

39 • Hourly  Wind  Direction  Data  Card 
On  this  card  are  input  the  wind  directions  for  each 
hour  of  the  simulation  period.  The  first  hourly  data  point 
is  defined  at  time  t = IBRH+1  (i.e.,  one  hour  into  the  simu- 
lation), and  the  total  number  of  values  specified  should 
correspond  to  the  number  of  simulation  hours  input  on  the 
Hourly  Data  Point  Number  Card.  The  wind  incidence  directions 
should  be  given  in  degrees,  measured  clockwise  from  North. 

If  the  simulation  is  to  cover  more  than  12  hours,  two  data 
cards  are  required. 

Variables : 

WDIR(I)  - wind  incidence  direction  at  time  I hours 

into  the  simulation  (in  degrees  from  Nor'h) 
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r 


t 


(card 

1) 

c . c . 

1-6 

F6.2 

WDIR ( 1 ) 

c . c . 

7-12 

F6 . 2 

WDIR(2 ) 

c . c . 

67-72 

F6.2 

WDIR ( 12 ) 

(card 

2)  (not  required 

if  NTEMP  < 12) 

c . c . 

1-6 

F6.2 

WDIR(13) 

c . c . 

7-12 

F6.2 

WDIR (14) 

• 

F6.2 

WDIR (NTEMP) 

Hourly  Electric 

Load  Factor  Data  Card 

card 

inputs  the 

hourly 

electric  load  . 

which  multiply  the  peak  community  load  specified  on  the 
Power  Plant  Rating  and  Peak  Non-Space-Conditioning  Electrical 
Demand  Data  Card  to  determine  the  non-space  conditioning  elec- 
tric loads  during  the  dynamic  system  simulation.  The  first 
value  specified  should  be  the  load  factor  at  time  t = IBHR+1 
(i.e.,  one  hour  into  the  simulation).  One  value  must  be 
given  for  each  of  the  data  points  specified  on  the  Hourly 
Data  Point  Number  Card.  The  program  uses  a linear  interpola- 
tion routine  to  determine  intermediate  loads  for  each  simula- 
tion time  step.  Therefore,  the  specified  load  factors  should 
accurately  reflect  the  actual  community  electrical  load  pro- 
file for  the  time  period  being  studied.  The  factors  should 
be  listed  consecutively  on  this  card  according  to  the  given 


A 
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input  format;  a maximum  of  12  values  per  card  and  a maximum 
of  2 cards  per  simulation  (i.e.,  a maximum  of  2'4  hourly 
data  points)  is  allowed. 

Variables : 

ELF (I ) - electrical  load  factor  (ratio  of  instantane- 
ous total  non-space -conditioning  electric 
load  to  peak  for  period)  at  time  I hours  into 


Format : 

the 

(card 

simulation 

1) 

(maximum  of  1= 2 h values). 

c . c . 

1-5 

F5.2 

ELF(l) 

c . c . 

6-10 

F5.2 

ELF ( 2 ) 

c . c . 

56-60 

F5.2 

ELF (12 ) 

(card 

2)  (not  required  if  NTEMP  < 12) 

c . c . 

1-5 

F5.2 

ELF (13) 

c . c . 

6-10 

F5.2 

ELF(liJ) 

- 

• 

F5.2 

ELF (NTEMP) 

(Max.  of  NTEMP=24) 

Hi. 

Hourly  Building 

Room 

Temperature  Data  Cards 

This 

card 

inputs  the 

room 

temperatures,  by  building 

type,  to  be  used  in  determining  the  building  space  conditioning 
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loads  for  the  dynamic  system  simulation.  Each  building  type 
is  assumed  to  have  a characteristic  room  temperature  which 
must  be  maintained  by  the  space  heating/cooling  systems 
throughout  all  the  rooms  of  that  building  type  in  the  commu- 
nity. These  temperatures  are  allowed  to  vary  on  an  hourly 
basis  to  correspond  to  the  changing  requirements  of  th^ 
people  using  the  buildings.  (e.g.,  Business  Offices  need 
not  be  fully  heated  during  unoccupied  periods).  The  program 
uses  a linear  interpolation  routine  to  determine  intermediate 
temperatures  for  each  simulation  time  step.  Therefore,  the 
specified  temperatures  should  accurately  reflect  the  actual 
desired  room  temperature  variations  for  the  time  period 
being  studied.  All  temperatures  must  be  specified  in  °F. 

One  temperature  per  card  is  required  for  each  of  the  building 
types  specified  on  the  Building  Type  Number  Card;  a maximum 
of  20  building  types  is  allowed,  and  if  the  number  of  types 
specified  exceeds  10,.  two  data  cards  must  be  input  per  simu- 
lation hour.  One  card  (or  pair  of  cards,  if  NTYPE  > 10)  must 
be  supplied  for  each  of  the  hourly  data  points  specified  on 
the  Hourly  Data  Point  Number  Card,  with  the  first  card(s)  in 
the  set  corresponding  to  time  t = IBHR+1  (one  hour  into  the 
simulation) . A maximum  of  24  data  points  per  building  type 


is  allowed. 
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Variables : 

TSB(I,NN)  - internal  room  temperature  at  time  I hours 
into  the  simulation  for  building  type  W 
(maximum  of  1=24,  NN=20)  (in  °F) 

Format:  (card  1) 


c . c . 

1-7 

F7.2 

TSB (1,1) 

c . c . 

8-14 

F7 .2 

TSB (1 , 2 ) 

c ,c . 

64-70 

F7 .2 

TSB (1,10) 

(card 

2,  if 

NTYPE  > 10) 

c . c . 

1-7 

F?  . 2 

TSB (1,11) 

c . c . 

8-14 

F7 .2 

TSB (1,12) 

• 

F7.2 

TSB  ( 1 , NTY, 

(card 

2,  if 

NTYPE  < 10) 

(card  3,  if 

c . c . 

1-7 

F7.2 

TSB (2,1) 

c . c . 

8-14 

F7.2 

TSB(2 , 2 ) 

c . c . 

64-70 

F7-2 

TSB (2,10) 

To  NTEMP  (maximum  of  24)  cards  or  pairs  of  cards, 


42 . Hourly  Building  Use  Factor  Data  Cards 
This  card  inputs,  for  each  building  type  specified 
or.  the  Building  Type  Number  Card,  the  building  use  factor 
which  multiplies  both  the  peak  ventilation  and  the  peak 
internal  heat  generation  rates  to  determine  the  instantaneous 
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hourly  values  of  these  parameters.  If  more  than  10  building 
types  are  specified,  two  cards  must  be  input  per  simulation 
hour.  The  first  card,  or  pair  of  cards,  in  the  data  set 
should  contain  data  observed  at  time  t = IBHR+1  (i.e.,  one 
.hour  into  the  simulation),  and  the  total  number  of  data 
points  per  building  type  should  correspond  to  the  number  of 
hours  specified  on  the  Hourly  Data  Point  Number  Card.  If 
the  observed  ventilation  and  internal  heating  schedules  for 
a given  building  type  differ  significantly,  attempts  should 
be  made  to  adjust  the  specified  peak  values  of  each  so  that 
the  net  effects  of  the  single  schedule  will  approximate  the 
measured  behavior. 

Variables : 

BUV (I,NN)  - building  use  factor  at  time  I hours  into 

the  simulation  for  building  type  NN  (maxi- 
mum of  1=24,  NN=20)  (decimal  fraction  of 
peak  ventilation  and  internal  heating  rates) 


(card 

1) 

c . c . 

1-7 

F7 . 2 

BUF (1,1) 

c .c . 

8-14 

F7.2 

BUF (1,2) 

c . c . 

64-70 

F7.2 

BUF (1,10) 

-r-r: 


(card 

2,  if 

NTYPE  > 10) 

c . c . 

1-7 

F7.2 

3UF (1,11) 

c . c . 

8-14 

F7.2 

BUF(1,12) 

• 

F7.2 

BUF(1, NTYPE) 

(card 

2,  if 

NTYPE  < 10) 

(card  3,  if  NTYPE 

c . c . 

1-7 

F7.2 

BUF (2,1) 

c . c . 

8-14 

F7.2 

BUF (2,2) 

c . c . 

64-70 

F7.2 

BUF (2,10) 

To  NTEMP  (maximum  of  24)  cards  or  pairs  of  cards. 

43 . Hourly  Domestic  Hot  Water  Use  Factor  Data  Cards 
On  this  card  are  input,  for  each  building  type  speci- 
fied on  the  Building  Type  Number  Card,  the  domestic  hot 
water  use  factors  which  multiply  the  peak  water  consumption 
rates  on  the  Building  Solar  Data,  Internal  Heating,  and  Domes- 
tic Hot  Water  Data  Cards  to  obtain  the  hourly  hot  water  ser- 
vice energy  demands.  If  the  number  of  specified  building 
types  exceeds  10,  two  cards  must  be  supplied  per  simulation 
hour.  The  first  set  of  data  points  should  be  specified  at 
time  t = IBHR+1  (one  hour  into  the  simulation),  and  a set  of 
use  factors  should  be  input  for  each  hour  specified  on  the 
Hourly  Data  Point  Number  Card. 


Variables : 


WUF(I,NN)  - domestic  hot  water  use  factor  at  time  I 
hours  into  the  simulation  for  building 
type  NN  (maximum  of  1=24,  NN=20)  (decimal 
fraction  of  peak  use  rate) 

Format:  (card  1) 


c . c . 

1-7 

F7.2 

WUF (1,1) 

c . c . 

8-14 

F7.2 

WUF (1,2) 

c . c . 

64-70 

F7.2 

WUF (1,10) 

(card 

2,  if 

NTYPE  > 10) 

c . c . 

1-7 

F7.2 

WUF (1,11) 

c . c . 

8-14 

F7.2 

WUF(1 ,12) 

• 

F7.2 

WUF (1, NTYPE) 

(card 

2,  if 

NTYPE  < 10) 

(card  3 , if  NTYPE  > 10 

c . c . 

1-7 

F7.2 

WUF (2,1) 

c .c . 

8-14 

F7 . 2 

WUF (2, 2) 

c . c . 

64-70 

F7 .2 

WUF (2,10) 

To  NTEMP  (maximum  of  24)  cards  or  pairs  of  cards. 


172. 


: 

t 


A . 5 . 2 Sample  24-Hour  Input;  Data  Listing 

To  familiarize  the  user  with  the  composition  of  the 
daza  deck,  a complete  listing  of  the  input  data  required  for 
a sample  2d -hour  TES  simulation  is  given  below. 

The  specific  input  data  presented  models  a thermal 
energy  distribution  system  consisting  of  a single  primary 
HTV.'  loop  and  three  secondary  distribution  loops  - the  first 
serving  three,  the  second  serving  four,  and  the  third  serving 
five  load  center  heat  exchangers.  The  primary  loop  contains 
approximately  10  miles  of  nominal  18"  O.D.  pipe,  while  the 
secondary  loops  each  have  12"  O.D.  pipe.  All  heat  exchangers 
contain  1"  O.D.  tubes  having  1"  I.D.  interstitial  flow  chan- 
nels. All  circulation  pumps  are  located  in  the  return  Un- 
to the  loop  supply  heat  exchangers.  The  system  is  supplied 
from  the  TUS  supply  heat  exchanger  at  a constant  water  *.  • r 

ature  of  380  °P;  consumer  end-use  equipment  is  supplied  at 
nominal  temperatures  of  220  °F  and  2d0  °F,  depending  upon  thr 
type  of  equipment,  the  location,  and  the  load  to  be  served. 

The  consumer  sectors  contain  a total  of  eleven  different 
building  types,  representing  a community  with  a resident  popu- 
lation of  approximately  50,000.  The  buildings  are  distributed 
throughout  the  community  in  load  centers  based  upon  their 
general  heating  and  cooling  requirements.  Also  taken  into 
consideration  in  this  distribution  are  community  power  plant 
siting  restrictions,  with  residential  areas  being  located 
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furthest  from  the  plant  site  and  low  constant  population 
industrial  sectors  situated  at  a minimum  radius  of  approxi- 
mately three  miles  from  the  station.  End-use  equipment  is 
distributed  throughout  the  community  following  practices 
expected  to  yield  an  optimum  thermal/electrical  energy 
demand  ratio  affording  a minimum  required  power  plant  in- 
stalled capacity.  Approximately  75#  of  the  total  energy 
demands  of  the  community  are  supplied  through  the  TUS. 

The  given  air  temperature  data  corresponds  to  a rela- 
tively severe  winter  day  at  Fayetteville,  North  Carolina. 

The  non-space-conditioning  electrical  load  schedule  was 
taken  from  demand  data  for  a minimum  space  conditioning  day 
at  Fort  Bragg,  North  Carolina,  and  the  magnitude  of  the 
average  load  was  scaled  to  match  the  model  community  size. 

To  illustrate  the  effects  of  weather  variations  during  the 
period,  building  room  temperature  requirements  are  held 
constant  for  each  building  type.  General  building  usage 
and  domestic  hot  water  consumption  schedules  follow  observed 
conditions  throughout  typical  24-hour  building  occupancy 
cycles . 

The  "card  type"  listing  shown  below  refers  to  the  numeri- 
cal designations  of  the  various  input  data  card  descriptions 
presented  in  Section  A.5-1.  These  numbers  do  not  appear  on 
any  of  the  actual  data  cards  and  are  given  here  merely  as  a 
reference  for  the  user  to  facilitate  identification  of  the 


required 

card  formats  and 

input  variables 

. AS  S. 

reference 

for  data 

position  on  each 

card,  the  s 

ho  Kn  on 

the  first 

numerical 

data  card  appear 

■s  in  column  3 • 

Ati  -3 1-  ■ 

u.  _ ■ J C. 

CL  ul  kV  “DO 

and  including  the  first  Primary  Program  Flow  Control  Code 
Card  (type  27)  is  required  for  every  run.  However,  data 
following  this  card,  and  certain  parameter  values  specified 
in  the  preceding  data  sets,  may  be  modified  according  to  th 
specific  program  option  to  be  executed  as  summarized  in  the 
following  sub-sections.  It  is  recommended  that  the  data 
format  presented  here  (i.e.,  use  of  the  steady-state  sub- 
routines to  initialize  the  dynamic  system  calculations  - a 
sequence  of  K0DE=0,  IDES=0,  K0DE=1)  be  used  in  all  dynamic 
simulations . 


$-4  P^*-h  tH  OsJ  vt\\o  vt\vO  'AvO  'Avo  *AvO  ^vO  vt\vo  vt\vo  vrwo  ltwo  lp>\q  va\p>  vo, 

cS  >> 


•— 

r- 

p 

r*- 

rvj 

r- 

>J 

vC 

r- 

rr» 

p 

P 

IX'« 

to 

tO 

P 

rr\ 

>* 

in 

o 

m 

nJ 

X 

>0 

O 

< 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

oO 

u. 

o 

C 

O 

o 

C 

c 

C 

O 

o 

o 

U. 

nr* 

V- 

K- 

>fr 

c 

<\ 

P 

p 

P 

■o 

P 

p 

P 

<r 

P 

G 

C 

o 

o 

G 

G 

o 

o 

G 

c 

o 

p 

*-> 

G 

r* 

/-) 

u. 

r~ 

i.  ’ 

-\; 

c. 

( 

v_' 

OJ 

l! 

— 

u 

rO 

u 

Li. 

•nj 

{_ 

n- 

r- 

G 

tO 

c 

'j 

P 

r\j 

c 

P 

G 

P 

S' 

•p 

r- 

— 4 

Kj 

<\! 

m 

P 

p 

_4 

p 

r : 

<\ 

• 

• 

r> 

• 

• 

r-' 

• 

rK 

• 

-f 

• 

• 

G 

• 

p 

• 

♦ 

• 

• 

v* 

X 

• 

P 

vf 

CO 

P 

e 

p 

<r 

h- 

X 

V' 

c 

p 

-- 

jn 

CT' 

c* 

c 

r 

P 

i 

> 

• 

IT- 

■ • 

IP 

• 

p 

• 

p 

• 

p 

• 

« 

>c 

• 

X 

• 

p 

• 

x 

• 

p 

• 

>C 

• 

<1 

f-4 

» 4 

r\ 

>*N. 

p 

P 

o 

O' 

X 

p 

f\ 

c\- 

G 

r- 

■o 

p 

'Ll 

p 

p 

o 

o 

o 

G 

X 

LL 

a 

>c 

P 

P 

r\ 

■■—4 

o 

o 

p 

— - 

— 

_-4 

cr 

, 

Lu 

>* 

c 

«— • 

a 

<\j 

o 

— < 

o 

m 

c 

G 

'r 

o 

r° 

a 

o 

^-4 

G 

G 

p 

— « 

a 

h- 

o 

• 

c 

• 

G 

• 

g 

• 

o 

• 

g 

• 

G 

• 

r_j 

• 

G 

• 

c: 

• 

o 

• 

G 

• 

G 

2 

u. 

•—* 

u. 

Li 

4J_ 

— 

X 

—4 

a_ 

—* 

Li. 

_4 

LL 

•—4 

LL 

i_ 

—4 

-L 

y_ 

G 

< 

•-* 

>r 

o 

c 

G 

r-« 

o 

X 

P 

c 

c 

o 

nT 

G 

oO 

c 

o 

c 

c 

P 

o 

rv 

o 

o 

c 

f\j 

C 

m 

c 

X 

c 

r 

o 

X 

c 

- 

a 

fV 

O' 

c 

X 

e 

r ' 

«x 

3 

tT' 

=•'• 

>o 

o 

p 

c 

o 

sC 

p 

c 

o 

o 

P 

P 

c 

• 

r- 

• 

c 

• 

• 

c 

• 

p- 

• 

G 

• 

O' 

t 

p*» 

• 

r*~ 

• 

• 

rr 

• 

X 

>r 

c 

o 

>r 

c 

C 

O 

vf 

c 

o 

c 

c 

O' 

c 

o 

c 

C 

cr 

o 

G 

%/ 

C 

a 

• 

X 

• 

• 

• 

• 

t 

• 

• 

• 

• 

• 

• 

• 

— 1 

-r 

c 

p— 

g 

— 

c 

— 

rr 

c 

c 

G 

c 

c 

c 

c 

3 

3 

3 

c 

c 

3 

N 

C 

'C 

>c 

p 

p 

p 

0 

p 

p 

p 

>t 

P 

c 

c 

c 

c 

c 

c 

c 

c 

c 

o 

C 

c 

c 

Lj 

-- 

li. 

c 

u. 

•JJ 

o 

U 

m 

L-J 

rr 

LL 

c 

L. 

p 

a1 

p 

a 

>0 

U 

o 

X 

p 

rr 

rr 

X 

X 

X 

c 

p 

p 

vt 

G 

ir 

p 

ir> 

p 

X 

p 

p 

c 

•** 

3 

v<»- 

G 

c 

p 

G 

o 

c 

0 

p 

rr 

G 

o 

G 

c 

x G 

r- 

-r 

— - 

p 

—4 

p 

G 

G 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

—4 

—4 

•— 

y»4 

— 

r»~ 

— • 

0 

X 

r "* 

f*- 

— * 

\ 

r*“ 

p 

P 

K- 

«* 

n. 

p 

p 

>C 

rO  CD 

U PxO  'A\0  'AvO  v%\0  VT\\c  ^A\0  ^VO  'AvO  UAvO  C^COOGCG(X)OOCOCOCOCOCOCOCOOOCO 

Ctf  >3 

O Eh 


r- 

r- 

f^- 

sO 

r» 

sT 

>7 

>r 

nT 

O' 

cr 

cr 

0* 

X 

O' 

3 

O 

o 

a 

o 

c 

o 

o 

o 

O 

^ ‘4 

—4 

r4 

c 

o 

•0 

V.C 

nD 

£ 

n 

o 

X 

o 

G 

c 

3 

o 

G 

G 

G 

G 

G 

u. 

^“4 

IL. 

—4 

—4 

, J 

cr 

• i 

__ 

u. 

G 

G 

L 

a 

l- 

G 

L 

G 

a 

X 

o 

o 

X 

c 

f^- 

g 

r- 

C 

nT 

<N 

fX 

nr 

rv 

-o 

X 

XI 

• 

e 

• 

X 

• 

**■ 

• 

<x 

• 

G 

« 

\ 

« 

*> 

♦ 

O' 

• 

<\ 

• 

3 

• 

r- 

X 

—4 

rn 

—4 

X 

— 

*n 

r\i 

x 

— 

rC\ 

■ -4 

X 

— 

cr 

m 

— < 

X 

o 

• 

>c 

t 

X 

t 

o 

« 

C 

• 

• 

o 

« 

o 

• 

O 

« 

G 

• 

4L 

~4 

g 

* 

o 

G 

G 

o 

c 

c 

T. 

c 

O' 

O' 

O' 

O' 

O' 

o 

o 

G 

o 

cr 

o 

O' 

O' 

O' 

O' 

O' 

o 

•—• 

3 

~4 

c 

—4 

o 

vj 

cr 

—> 

o 

^4 

cr 

~4 

O' 

*^4 

O' 

«— i 

o 

•—4 

o 

• 

C5 

• 

G 

• 

G 

« 

rj 

• 

o 

• 

G 

• 

3 

• 

3 

« 

3 

• 

G 

t 

u 

Uj 

—4 

(j- 

—4 

LL 

G 

u 

G 

a. 

3 

u 

o 

ul. 

o 

o 

o 

o 

o 

cr 

g 

C 

a 

O 

G 

C 

c 

Q 

c 

in 

c 

C 

o 

c 

O 

O 

o 

X 

a 

X 

c 

X 

c 

j 

o 

~ 

o 

O 

X 

G 

X 

c 

~ 

0 

X 

G 

v*ina'^fsi\XsO 

—4  ^ oc3-cccv^oor'->r:oo 


<?  r-  x — * r\  -G  ^ r\j  m *r  — - ^ m vT 
^h'n;,*^— * — ^HfNjryj'vrv,  -n  m <g  ro 
<\'\\"Crr^ff)fnfOPOrCr«f^,,fi 


3 O O O'  ^ ro  ro  G r<* 

« X.  • *’■'•  • S • O • CZ  • "!C  • X # X • CL  • X • 

COO^CCCCTCCOCCGCCOCCC^C  r\sf  so 

• • • • • • rr  • x • x • a.  • c: 

cccocccoc  o 


C C C sj  c c c 

IT  IT  't  in  'S'  If* 

o c r c c c 

a,  •'O  LLJ  >?■  Li.  iP  u-  O u.  G l:  L.  \ U. 

*A  X -**-*  — rs,  x c — < — 

O rr  X -*  sj- 

*-  — -K-—  t— -- 
• ••••• 

^ X .%  a O' 

>r  it  £ r-  x — « <\. 


C c C c r\.'Jsu<\.Gx*<\iro.  ^incvjrn^m 

• x • c: 

C C 

^^^r\jPNkf\j<\jf^r\jfs.rg(\fNkfv. 

't  li.  IT  LL 

O m lp  — 4 rv  ^ -o 


1200 


90.0  2?noij. 

0.40  0.9S  0.9r>  0.9S  0.9S  1.00  0.91  0.  HO  007. Ot>  I . 1 j OF  06 


Card 


CM  CM  CM  CM  CM  CM  CM  CM  CM  C\1  CM  CM  CM  CM  CM  CM  C'YM'  vr^yr,  ocDOOrHCV  crA  frA.--v  Vf\  VT'  fx.vf}  r\_ 

>>  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  C'\  CM  C ACr'\  CM  C*S*?S  C*A  c~\  CM'  C.  C'  . tr 

Eh 


O rft 

CM  rO 


o lO 

IP  — 1 

O • 


o 

o 

o 

*— 

cc 

IP 

o 

c 

3 

o 

• 

« 

c 

c 

4- 

— * 

• 

• 

• 

IP 

"M 

• 

>c 

—4 

c 

• 

c 

CO 

t 

c 

c 

LP 

c 

• 

r- 

Cy 

o 

• 

1 

IP 

o 

c 

— 

c 

— * 

• 

• 

• 

— - 

IP 

• -4  p* 

CM  CC  — < 

N • » 


• • — - C 

• SC  n cc  f\  rr 

• • C M • • • • 

csj  ip  ^ l r 


C ard 


C.  C T 

vC  v-  ' 


? x r? 


:ar 


189. 


A.  5- 3 Input  Data  Modifications 


A. 5. 3.1  Building  Energy  Demands 
When  performing  only  the  single  time  instant 
building  load  analysis  (KODE=0,  IDES=1,  K0DE=2),  it  is  possi- 
ble to  eliminate  many  of  the  parameter  specifications  on 
the  cards  describing  the  TUS (types  5 and  6),  since  the  only 
information  required  is  the  general  TUS  piping  configuration 
and  the  locations  of  the  load  center  heat  exchangers.  Thus, 
all  fluid  flowrates,  fluid  parameters,  and  nodal  dimensions 
may  be  omitted,  although  it  may  be  desirable  to  input  non- 
zero dummy  values  for  these  parameters  to  eliminate  the  possi- 
bility of  computational  overflow  errors  occurring  during  the 
transition  from  the  single  time  analysis  to  an  hourly  load 
analysis.  The  consumer  fluid  supply  and  return  temperatures 
and  the  nodal  fluid  specific  heats  must  be  specified  if  a 
complete  analysis  is  desired,  since  the  initial  load  center 
heat  exchanger  fluid  mass  flowrates  are  computed  at  the  same 
time  the  load  distribution  is  printed.  All  TUS  nodal  cards 
must  be  input,  along  with  their  location  and  type  codes,  to 
facilitate  heat  exchanger  identification  by  the  load  aggrega- 
tion section  of  the  program.  The  number  of  tubes  in  each  heat 
exchanger  need  not  be  specified  (card  type  8),  but  its 
supply-and  demand-side  node  location  codes  must  be  supplied. 
Initial  conditions  at  the  reservoir  and  on  the  power  plant 
side  of  the  TJS  supply  heat  exchanger  are  unnecessary,  but 


lummy  values  should  be  input  to  prevent  possible  ov 
during  the  initialization  calculations.  -.11  tuildi 


ing  para- 


meter specifications  ana  all  building  distribution  ia-a  nust 
e supplied  to  every  program  run  option.  Finally,  for 
■ingle  time  analyses,  all  initial  system  data  shown  in  the 
sample  listing  between  IDES=0  (card  type  28)  and  K0DE=1 
(card  type  27)  must  be  input  following  ir'nS=~  to  supply  the 
load  calculation  subroutines  with  appropriate  weather  and 
building  requirements  specifications.  (Of  course,  KCDE=2 
ould  then  replace  KODE=i  following  card  type  35,  and  no 
further  data  would  be  required). 

Extension  of  the  single  time  instant  building  load  ana- 
lysis to  an  hourly  analysis  without  performing  any  Tl'S  cal- 
culations requires  non-zero  (and  non-equal)  values  to  be 
supplied  for  the  consumer  fluid  supply  and  return  temperatures 
and  non-zero  values  for  the  load  center  heat  exchanger  fluid 
specific  heats  to  avoid  overflows  during  the  calculation 
period  between  the  tv.’O  option  executions,  following  speci-t- 
cation  of  K0DE=1  after  the  initial  weather  data  and  building 
conditions  are  input,  data  must  be  provided  for  each  hour 
of  the  analysis  desired.  (These  hourly  data  values  need  not 
be  consistent  with  continuous  functions  for  a specified  time 
- eriod,  since  the  demands  are  computed  independently  at  each 
hourly  interval). 
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A. 5- 3- 2 Steady-State  System  Analysis 
Since  the  steady-state  TUS  analysis  performs 
essentially  the  same  computations  as  does  the  dynamic  system 
analysis  (but  under  steady-state  energy  and  momentum  trans- 
fer conditions),  the  input  data  requirements  for  the  steady- 
state  analysis  program  execution  option  are  exactly  the  same 
as  those  shown  for  the  dynamic  system  analysis,  with  the 
single  exception  that  K0DE=2  is  input  for  the  second  value 
of  KODE  (card  type  27)  and  no  hourly  data  is  required.  Esti- 
mates of  the  secondary  loop  fluid  flowrates  should  be  approxi- 
mately equal  to  the  expected  converged  system  values  to  mini- 
mize the  number  of  steady-state  iterations  required.  As  is 
demonstrated  by  the  sample  data,  all  secondary  loop  tempera- 
tures, all  consumer-side  load  center  heat  exchanger  fluid 
flowrates,  and  all  primary  HTW  loop  temperatures  except  for 
the  fixed  supply  temperature  at  the  loop  inlet  may  be  left 
unspecified  with  no  adverse  effects  upon  system  convergence 
or  the  overall  cost  of  the  calculations.  Since  typical  TUS 
configurations  may  be  expected  to  converge  well  within  1 CPU- 
minute  of  computing  time  if  the  estimated  initial  conditions 
are  at  all  reasonable  (see  Table  A. 10),  use  of  the  steady- 
state  system  analysis  with  a specified  upper  limit  of  1 
minute  for  the  CPU  time  can  provide  a means  for  verifying  the 
validity  of  all  but  the  hourly  input  data  directly  associated 
with  the  dynamic  system  analysis.  Any  mis-punched  data  cards 
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Whenever  IDES=1  in  input,  the  first  system  analysis  data 


to  be  output  is  a summary  of  the  individual  building  type 


space  conditioning  loads  calculated  for  the  specified  initial 


weather  and  building  usage  conditions  (i.e.,  the  loads  at 


time  t = IBHR) . These  individual  unit  demands  are  then 


aggregated  according  to  the  specified  distributions  o: 


buildings  and  end-use  equipment,  and  the  resulting  load  dis- 


tribution is  printed  by  building  type  for  each  load  center 


heat  exchanger  and  for  those  buildings  m the  "external 


system"  not  supplied  by  the  TUS.  (It  should  be  noted  that 


the  printed  loads  are  the  net  heat  losses/gains  for  the 


buildings  prior  to  being  scaled  by  the  end-use  equipment  COPs 


and  that  they  include  domestic  hot  water  demands).  For  TUS 


heat  exchanger  design  and  system  initialization  purposes,  the 


total  thermal  energy  demands  at  each  load  center  heat  exchan- 


ger (including  the  effects  of  the  consumer  equipment  COPs, 


but  not  of  heat  losses  in  the  piping)  are  printed.  If  an 


hourly  building  analysis  is  desired,  the  individual  building 


type  space  conditioning  load  summary  is  printed  for  every 


set  of  hourly  data  points  provided,  but  no  further  distribu- 


tion data  is  supplied.  Shown  below  are  the  outputs  of  the 


building  space  heating  demands,  the  building  load  distribution 


and  the  load  center  heat  exchanger  thermal  demands  correspon- 


ding to  the  initial  conditions  specified  In  the  sample  input 


data 
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EXCHANGER 

4 


5 

6 

7 

8 
9 

10 

1 1 - 

12 

13 

14 
--  15 


LOAD  i BTU/HR  J 
0.4C9LF+06 
C.  2E52F+C8 
0. 2378E+C3 
C.  1653F+08 
C.  22C5E*08 

0.35A1E+C8  

0.4C14E+08 

0.  1515E  + 0e 

0. 44 80E+08 
0. 4514F  + C8 
0.  1374E+08 
C. 3422F+G8--- 


**TERVlNATION  CODE  READ** 
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A . 6 . 3 Steady-State  System  Analysis 

When  the  steady-state  system  analysis  is  executed,  the 
individual  building  type  space  conditioning  loads  are  com- 
puted and  printed  as  in  the  building  load  analysis  option 
described  above.  The  demands  are  then  aggregated,  scaled 
by  the  equipment  COPs,  and  are  applied  to  the  load  center 
heat  exchangers.  When  TUS  convergence  is  achieved,  a nodal 
parameter  summary  is  printed  which  includes  fluid  temperatures 
fluid  mass  flowrates,  the  steady-state  nodal  energy  flows, 
loop  circulation  pump  pressure  settings,  fluid  conditions 
on  both  sides  of  the  TUS  supply  heat  exchanger,  and  speci- 
fications of  the  heat  stored  in  each  section  of  the  thermal 
reservoir.  (It  should  be  remembered  that  the  values  printed 
for  the  reservoir  stored  energy  are  only  relative  indicators 
of  the  actual  energy  storage  trends,  and  that  the  wat.r  volume 
data  provides  the  absolute  heat  storage  information).  Follow- 
ing this  TUS  summary,  the  building  load  distribution  on  a 
per-heat-exchanger  basis  is  printed  for  the  steady-state 
demand  conditions,  and  overall  TUS  and  power  plant  energy 
supply/demand  data  is  summarized.  The  building  loads  and 
the  demand  distribution  data  output  for  this  option  are 
identical  to  those  illustrated  in  Section  A. 6. 2.  The  con- 
verged steady-state  TUS  nodal  summary  corresponding  to  the 
sample  initial  conditions  is  shown  below. 


, c.-  Dynan. 


System  Analys 


When  the  lynamic  system  analysis  Is  initialized  inter- 
lly  from  the  converged  steady-state  TUS  parameters,  the 
first  block  of  analysis  output  data  corresponds  to  the 
printouts  summarized  in  Sections  k.6.2  and  A. 6. 3 for  the 
building  load  analysis  and  steady-state  TUS  analysis  program 
execution  options.  (In  this  case,  the  hourly  building  loaa 
summary  and  its  associated  input  data  are  printed  following 
the  steady-state  system  output,  since  the  hourly  data  is  not 
read  into  the  program  until  the  TUS  parameters  have  con- 
verged . ) 

The  output  from  the  dynamic  system  simulation  is  schedul- 
ed by  the  user  and  can  be  supplied  as  frequently  as  every 
time  step  or  as  infrequently  as  once  each  hour.  Printed  out- 
put includes  the  instantaneous  TUS  nodal  summary  at  the  end 
of  the  given  time  step,  the  instantaneous  distribution  of 
building  loads  at  each  load  center  heat  exchanger,  and  a 
summary  of  the  thermal  energy  supply/demand  mismatches  at 
each  heat  exchanger  in  the  system.  The  nodal  summary  and 
building  load' distribution  are  identical  in  form  to  the 
output  illustrated  for  the  steady-state  analysis;  a sample 
heat  exchanger  energy  balance  output  listing  is  shown  below. 


Also  output  from  the  dynamic  simulation  is  a punched 
card  deck  containing  specifications  of  the  instantaneous 
TJS  total  thermal  load  (including  all  losses),  the  instan- 
taneous TES  total  electrical  demand,  and  the  volume  of  water 
contained  in  the  hot  section  of  the  thermal  energy  storage 
reservoir  at  The  end  of  each  output  interval  time  step. 

These  overall  system  parameters  are  output  in  a format  easily 
adapted  tc  integrated  energy  consumption  analyses  (see,  for 
example,  Reference  14)  or  to  simple  hourly  system  summary 
data  plotting  routines.  The  composition  and  format  of  this 
data  deck  is: 

1.  Header  card: 

c.c.  1-7  F7 . 2 power  plant  electrical 

rating  in  MWe 

c.c.  8-12  F5-2  power  plant  thermal/elec' 

trical  energy  conversion 
efficiency  (decimal 
fraction) 

2.  System  data  cards  (one  punched  for  initial  conditions 

. and  at  each  specified  output  interval  during  dynamic 

analysis ) 

c.c.  1-5  F8.2  total  TUS  thermal  load 

in  MWt 

c.c.  12-19  F8 . 2 total  TES  electrical 


load  in  MWe 
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f; 


water  volume  in  reservoir 
hot  water  section  in 
102*  ft3* 


*Por  most  moderately-sized  systems,  this  unit  of  storage 
produces  punched  values  for  the  water  volume  of  magnitudes 
comparable  to  those  for  the  thermal  and  electrical  energy 
demands  and,  therefore,  easily  adapted  to  a single  set  of 
plotting  axes. 


I 


TDIST  Variable  Definitions 


Variable 
A ( J ) 

AC 

ACC (NP ) 

AEC 

AF  ( I ) 

AINFL(I) 

ALAT 

ALEAK 

ALPL 

AMAS(K,L,N 

AMS(L) 

AR(I) 

ASC 

ATD ( K , L , N ) 


Definition 

Coefficient  in  window  solar  heat  gain  formula 
for  absorption  of  incident  radiation 
Nodal  cross-sectional  flow  channel  area  (ft2) 
Fluid  acceleration  through  pump  NP  during 
one  time  step  (Ibm/hr  ) 

Atmospheric  extinction  coefficient  (dimension 
less ) 

Area  of  Building  Type  I in  contact  with  the 
earth  (ft^) 

Maximum  rate  of  forced-air  ventilation  for 
Building  Type  I (CFM) 

Latitude  of  site  (degrees) 

Total  air  exchanger  between  interior  and 
exterior  of  building  ( CFH  or  CFM) 

Heat  exchanger  fluid  temperature  decay  con- 
stant (dimensionless) 

Desired  nodal  fluid  mass  flowrate  (lbm/hr) 
Desired  fluid  mass  flowrate  in  secondary  loop 
L (lbm/hr) 

p 

Roof  area  of  Building  Type  I (ft  ) 

Apparent  solar  constant  (BTU/hr-ft2) 

Linear  average  difference  between  nodal  fluid 
temperature  and  that  of  surroundings  (°F) 


cro 


odal  fluid  outlet  temperature  (°P) 


-■■dal  heat  transfer  coefficient  (BTU/hr-°F) 


;round- contact  material  conduction 


coefficient  (BTU/hr-°F) 


Tota1  wall  and  window  conduction  heat  transfer 


icient  (BTU/hr-°F) 

exposed  wall  area  of  Building  Type  I (ft^) 
window  area  of  Building  Type  I (ft^) 


( 3TU/hr-°F ) 


conditioner 


formance 


ipment  coefficient  of  performance 


3 pmen 


ture-ad justed  compressive  air  conditioner 


er formance 


lding  Type  I (ft) 


_ 


217. 


BHPM 


BMAS (1,1,N) 
BPDIS 
BUF ( I , NN ) 


BUF (25,1) 


CAP ( K , L, N ) 
CAPG 

CITY 

COND(K,L,N) 

CONV(K,L,N) 

CS 

CSFAC 


Temperature-adjusted  heat  pump  coefficient 
of  performance 

Desired  fluid  mass  flowrate  in  bypass  around 
primary  loop  heat  exchanger  node  N (lbm/hr) 
Inverse  sum  of  parallel  node  pressure  drop 
factors  (ft'Vhr-lbf) 

Lse  factor  of  Building  Type  NN  at  hour  I of 
simulation  (decimal  fraction) 

0 = Building  unoccupied,  no  equipment  use 

1 = Full  occupancy,  all  lights,  appliances 

ventilation  in  use 

Use  factor  of  Building  Type  I at  beginning  of 
simulation  (decimal  fraction) 

0 = Building  unoccupied,  no  equipment  use 

1 = Full  occupancy,  all  lights,  appli-nces 

ventilation  in  use 

Nodal  fluid  specific  heat  (BTU/lbm-°F) 
Specific  heat  of  power  plant  coolant  in  pre- 
cooler (BTU/lbm-°F) 

Alphanumeric  run  title 

Nodal  heat  conduction  rate  (BTU/hr) 

Nodal  heat  convection  rate  (BTU/hr) 

Solar  incidence  direction  cosine 

Case  input  mass  flowrate  scale  factor  (1.0) 
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CTH(I) 


cvo 

CVOI 

cw 

cz 

DCL(I) 

DECL 

DENOM 

DPSI(NN) 

DFSIR 

DHW(I) 

DiA ( K , L , N ) 
DINC(I) 

DINN(I) 


'osine  cf  incident  solar  radiation  relative 


,o  norma. 


suriace  _ 


I = 1-4  * walls,  clockwise  from  north 
I = 5 = roof 

Volume  of  cold  water  section  of  reservoir  ( ft 3 ) 

Volume  of  cold  water  section  of  reservoir  (ft3) 

Solar  incidence  direction  cosine 

Solar  incidence  direction  cosine 

Total  length  of  cracks  around  all  exterior  doors 


in  Building  Type  I (ft) 

Solar  declination  (degrees) 

Heat  exchanger  heat  balance  coefficient  (dimen- 
sionless) 

Diffuse  solar  radiation  intensity  incident  on 
wall  NN  (BTU/hr-ft2) 

Diffuse  solar  radiation  intensity  incident  on 
roof  (BTU/hr-ft2) 

Maximum  energy  demand  for  domestic  hot  water  use 
in  Building  Type  I (BTU/hr) 

Nodal  flow  channel  diameter  (ft) 

Infiltration  air  flow  coefficient  C for  doors  in 
Building  Type  I (CFM/ft-in.  H2O) 

Infiltration  air  flow  coefficient  N f^r  doors  in 
Building  Type  1 (dimensionless) 
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DL  Distance  fluid  flows  during  one  time  step  (ft) 

DLK  Total  air  infiltration  through  cracks  around 

doors  (CFM) 

DM  Loop  fluid  mass  flowrate  at  end  of  time  step 

( lbm/hr ) (FLOWS) 

DM  Nodal  fluid  mass  flowrate  (lbm/hr)  (PRESS) 

DMAS ( 1 , 1 , N ) Desired  fluid  mass  flowrace  in  primary  loop 


DNSI 


DPM(L) 


heat  exchanger  node  N (lbm/hr) 

Direct  normal  solar  intensity  (BTU/hr-ft^) 

P 

Nodal  fluid  frictional  pressure  loss  (lbf/ft  ) 
Maximum  free-flow  parallel  node  pressure  loss 
for  loop  L (lbf/ft^) 


DQ ( ME ) 


Thermal  demand  at  heat  exchanger  NE  (BT'U/hr) 


DSI(NN) 


Direct  solar  radiation  intensity  incident  on 


DSIR 


wall  NN  ( BTU/hr-f t ; 

Direct  solar  radiation  intensity  incident  on 
roof  (BTU/hr-ft2) 

Nodal  outlet  increment  fluid  temperature  loss 
due  to  heat  conduction  to  the  surroundings  (°F) 


DTI  ME 


Size  of  one  time  step  (decimal  fraction  of  an 


hour ) 


DTLM 


Nodal  log-mean  temperature  difference  (°F) 


DTI. MX 


Log-mean  temperature  difference  over  heat 


exchanger  outlet  fluid  flow  increment  length  (°F) 
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EAES  .Nominal  coefficient  of  performance  of  absorp- 

tive air  conditioners 

ECMP  Nominal  coefficient  of  performance  of  compres- 

sive air  conditioners 

EPF  Power  plant  thermal/electrical  conversion 

efficiency  (%) 

EHPM  Nominal  coefficient  of  performance  of  heat  pump 

EHWS  Nominal  efficiency  of  hot  water  heating 

ELBE  Peak  non-space-conditioning  electrical  demand 

(MWe) 

ELF(I)  Non-space-conditioning  electric  load  factor  at 

hour  I of  simulation  (decimal  fraction  of  peak 
load ) 

ELF(25)  Non-space-conditioning  electric  load  factor  at 

beginning  of  simulation  (decimal  fraction  of 
peak  load) 

ERE3  Nominal  efficiency  of  resistance  heaters 

ESCL  Total  system  electric  space  conditioning  load 

(BTU/hr ) 

ETHL(IHR)  Total  system  electric  space  conditioning  load 
at  simulation  hour  IHR  (BTU/hr) 

F Mooiy-Veisbach  friction  factor  (dimensionless ) 

FAC1-FAC4  Coefficients  in  expansion  for  window  solar  heat 
gain  formula 
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PCTR 

FF(IX) 


FHPEX(I) 


FP 


FRCN 


FREXT(I) 


FV(IX) 


FWLIT( I ,NN ) 


G 

H 


Conversion  factor  ( lbm-f t^/lbf-hr^ ) 

Friction  factor  for  heat  exchanger  IX  exclu- 
sive of  valve  resistance  (dimensionless) 

Fraction  of  units  of  Building  Type  I not 
served  by  the  Thermal  Utility  System  which  are 
supplied  by  heat  pumps 

Nodal  fluid  friction  factor  including  valve 
resistance  (dimensionless) 

Incremental/total  heat  exchanger  heat  flow 
ratio  (dimensionless) 

Fraction  of  units  of  Building  Type  I not  served 
by  the  Thermal  Utility  System  which  are  supplied 
by  resistance  heaters 

Valve  friction  factor  for  heat  exchanger  IX 
(dimensionless ) 

Fraction  of  surface  NN  lit  by  direct  sunlight 
for  Building  Type  I (decimal  fraction) 

NN:  1-4  = walls,  with  1 defined  by  normal 

direction  WA,  numbered  clockwise  from 
north 
5 = roof 

Conversion  factor  ( lbm-ft/lbf-hr^ ) 

Angular  position  of  solar  disk  relative  to  noon 
normal  (radians  from  normal) 


i 
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HRF 


HVO 


KVOI 


IBEG 


IBHR 


IDES 


IEND 

IERR 


IEX(K,L,N) 


IPLO 


I HR 

IHRL 

IIN 


Elapsed  fraction  cf  hourly  interval  being 
examined 

Volume  of  hot  water  section  of  reservoir  (ft3) 
Volume  of  hot  water  section  of  reservoir  (ft^) 
Beginning  hour  of  simulation  period 
Beginning  hour  of  the  simulation  period  (0-23) 

0 = midnight 
12  = noon 

Design  option  code 

C = perform  standard  utility  system  analysis 

1 = perform  building  load  analysis  only 
Endpoint  of  simulation  period 

Data  error  code 

0 = no  error 

1 = error 

Nodal  heat  exchanger  code 
0 = no  heat  exchanger 
1-30  = heat  exchanger 
Fluid  flow  calculation  code 

0 = pump  reset  or  initialization 

1 = normal  time  step  sequence 

End  point  of  hourly  Interval  being  examined 
Beginning  of  hourly  interval  being  examined 
Computer  system  card  reader  input  device  number 


INIT 


IOUT 

IPMP(K,L,N) 


IPNCH 

IPRE 

IPRM 

IPRX 

ISER(K,L,N) 


ITIM 

JPM 


Kl(NEX) 
K2 (NEX) 
KODE 


Pump  pressure  initialization  code 

0 = compute  loop  pressure  losses 

1 = reset  flow  valves  and  pump  settings 
Computer  system  line  printer  output  device 
number 

Nodal  pump  code 

0 = no  pump 

1 = pump 

Computer  system  card  punch  output  device  number 
Unused  variable 
Unused  variable 
Unused  variable 

Nodal  series/parallel  flow  code 

0 = series 

1 = parallel 

Hour  of  simulation  being  studied 
Thermostat  limit  indicator 

0 = limits  not  exceeded 

1 = limits  exceeded,  pump  reset  to  take  place 
Heat  exchanger  NEX  supply-side  group  index 
Heat  exchanger  NEX  demand-side  group  index 
Simulation  control  code 

0 = steady-state  analysis 

1 -*  transient  analysis 
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KODP 


Ll(NEX) 

L2(NEX) 

LONG(K,L,N 

LS 

M 

MASS (K,L,N 
MGAS 

MTR(I) 

MTW(I) 


mu(k,l,n) 

Ni(NEX) 

N2(NEX) 

NA 

NABS ( IX , I ) 


Consumer  summary  print  code 
C = no  print 
1 = print 

Heat  exchanger  HEX  supply-side  loop  index 
Heat  exchanger  NEX  demand-side  loop  index 
) Modal  flow  channel  length  (ft) 

Number  of  secondary  distribution  loops 
Primary  loop  node  indicator 
) Modal  fluid  mass  flowrate  (lbm/hr) 

Fower  plant  coolant  mass  flowrate  through 
precooler  (lbm/hr) 

Roof  surface  material  type  for  Building  Type  I 
(surface  roughness  classes  identical  to  those 
for  MTW(I)) 

Wall  exterior  surface  material  type  for  Building 
Type  I 

1 = stucco  3 = concrete  5 = plaster 

2 = brick  4 = pine  6 = glass 

Nodal  fluid  viscosity  (lbm/hr-ft) 

Heat  exchanger  NEX  supply-side  node  index 
Heat  exchanger  NEX  demand-side  node  index 
Nodal  indicator  number 

Number  of  units  of  Building  Type  I served  by 
heat  exchanger  IX  which  have  absorptive  air 


conditioning 


NB(L) 

N30I 

NCMP ( IX , I ) 

ND(  L) 

NDOI 

NDUM 

NEX 

NEXCH 

NH 

NHPM( IX , I ) 
NHR 

nhws(ix,i) 

NITR 

NMK 

NNODE 

NP 
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Number  of  consumer  load  heat  exchangers  served 
by  loop  L 

Dummy  variable  for  number  of  consumer  heat 
exchangers  served  by  a secondary  loop 
Number  of  units  of  Building  Type  I served  by 
heat  exchanger  IX  which  have  compressive  air 
conditioning 

Number  of  nodes  in  secondary  loop  L 
Dummy  variable  for  number  of  nodes  in  a secon- 
dary loop 

Load  center  heat  exchanger  supply  node  indicator 
Heat  exchanger  number 

Total  number  of  heat  exchangers  in  utility 
system 

Number  of  nodes  in  primary  HTW  loop 

Number  of  units  of  Building  Type  I served  by 

heat  exchanger  IX  which  have  heat  pumps 

Endpoint  of  hourly  interval  being  examined 

Number  of  units  of  Building  Type  I served  by 

heat  exchanger  IX  which  have  hot  water  heat 

Number  of  print  intervals  per  hour 

Load  center  heat  exchanger  supply  node  indicator 

Dummy  variable,  unused 

Secondary  loop  node  indicator 


| 
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NFMP 

NPRKT 

;-:res(ix,i) 


NSTEP 
NT3(NEX) 
NTEMP 
NTYPE 
NU ( IX , I ) 

NUEX(I) 

PC 

PC(I) 

PCI 

PEX(I) 

PMF 

PI^FI 

POL(I) 

PP(I) 


Number  of  pumps  in  utility  system 

Number  of  time  steps  between  print  intervals 

Number  of  units  of  Building  Type  I served  by 

heat  exchanger  IX  which  have  resistance  heaters 

Number  of  time  steps  per  hour 

Number  of  tubes  in  heat  exchanger  NEX 

Number  of  hours  in  simulation  period 

Number  of  building  types  in  community 

Number  of  units  of  Building  Type  I served  by 

heat  exchanger  IX 

Number  of  units  of  Building  Type  I not  served 
by  the  Thermal  Utility  System 
Stack-effect  pressure  on  building  (in.  E^O) 
Percent  change  in  fluid  mass  flowrate  through 
pump  I during  one  time  step 

Percent  increase  in  loop  fluid  mass  flowrate 
during  one  time  step 

Compensated  pressure  setting  for  pump  I (Ibf/ft^) 
Precooler  fluid  mass  flowrate  (Ibm/hr) 

Precooler  fluid  mass  flowrate  (lbm/hr) 

Maximum  rate  of  heat  generation  due  to  lighting, 
appliances,  people,  etc.  for  Building  Type  I (KW) 
Total  fluid  frictional  pressure  head  faced  by 
pump  i (lbf/ft‘~) 


PREF(I) 

PS(I) 

PSENS 

PT 

PTKN 

PTKO 

PV 

PW 

Q 

QABS 

QAEXT 


QCEXT 


QCI 


QCMP 


Reference  pressure  setting  for  pump  I (lbf/ft^) 

O 

Pressure  setting  for  pump  I (lbf/ft  ) 

Pump  response  sensitivity  factor 

Total  pressure  drop  across  wall  (in.  HpO) 

Normal  air  flow  pressure  drop  coefficient 
(dimensionless) 

Oblique  air  flow  pressure  drop  coefficient 
(dimensionless ) 

Normal  wind-induced  pressure  loading  on  wall 
(in.  P^O) 

Wind-induced  pressure  loading  on  wall  (in.  H20) 
Total  heat  loss  from  buildings  of  one.  type  served 
by  a given  load  center  heat  exchanger  (BTU/hr) 
Energy  demand  from  absorptive  air  conditioners 
(BTU/hr) 

Energy  demand  from  absorptive  air  conditioners 
in  buildings  not  served  by  Thermal  Utility  System 
(BTU/hr) 

Energy  demand  from  compressive  air  conditioners 
in  buildings  not  served  by  Thermal  Utility  System 
(BTU/hr) 

Heat  stored  in  reservoir  cold  water  section 
(BTU  referenced  to  0 °P) 

Energy  demand  from  compressive  air  conditioners 
( BTU/hr ) 
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QCO 

QCOKD(I) 

QCW 

QD 

QDEXT 

QDHW 

^DIF 

QDO(NX) 

QEXT 

QHEXT 

QHI 

QHO 

QHPM 


.Heat  stored  in  reservoir  cold  water  section 
(BTU  referenced  tc  0 °F) 

Total  heat  loss  due  to  conduction  from  3uilding 
Type  I (BTU/hr ) 

Heat  stored  in  reservoir  cold  water  section 
(BTU  referenced  to  C °F) 

Thermal  demand  at  heat  exchanger  (BTU/hr) 

Total  domestic  hot  water  energy  demand  for 
buildings  of  a given  type  not  served  by  Thermal 
Utility  System  (BTU/hr) 

Total  domestic  hot  water  energy  demand  for  build- 
ings of  one  type  served  by  a given  load  center 
heat  exchanger  (BTU/hr) 

Supply /demand  imbalance  at  heat  exchanger  (ETU/hr) 
Thermal  demand  at  heat  exchanger  NX  (BTU/hr) 

Total  heat  loss  from  all  buildings  of  a given 
type  not  served  by  Thermal  Utility  System  (BTU/hr) 
Energy  demand  from  heat  pumps  in  buildings  not 
served  by  Thermal  Utility  System  (BTU/hr) 

Heat  stored  in  reservoir  hot  water  section 
(BTU  referenced  to  0 °F) 

Heat  stored  in  reservoir  hot  water  section 
(BTU  referenced  to  0 CF) 

Energy  demand  from  heat  pumps  (BTU/hr) 


A 


kd 


Thermal  energy  available  to  thermal  utility 


system  after  electrical  generation  (BTU/hr) 


Heat  stored  m reservoir  hot  water  section 


(BTU  referenced  to  0 F) 


Energy  demand  from  hot  water  heating  in  build 


QHWEX 


ngs  not  served  by  Thermal  Utility  System  (BTU/hr 


Energy  demand  from  hot  water  heat  (BTU/hr) 


QHWS 


QL(IHR,I) 


Total  heat  loss  from  Building  Type  I at  simula- 


tion hour  IHR  (BTU/hr) 


QLEAK(I) 


Heat  loss  rate  from  Building  Type  I due  to  infil 


tration/ventilation  (BTU/hr) 


QLIG(I) 


Internal  heat  generation  rate  for  Building  Tyoe 


due  to  non-space  conditioning  equipment  and 


people  (BTU/hr) 


QRES 


Energy  demand  from  resistance  heaters  (BTU/hr) 


QREXT 


Energy  demand  from  resistance  heaters  in  build- 


ings not  served  by  Thermal  Utility  System  (BTU/hr) 


Thermal  supply  at  heat  exchanger  (BTU/hr) 


QS(L) 


Total  thermal  demand  at  load  center  heat  exchan 


gers  in  loop  L (BTU/hr) 


QTLD 


Total  thermal  demand  at  load  center  heat  exchan- 


gers (BTU/hr) 


QTLLD 


Total  thermal  demand  at  load  center  heat  exchan- 


gers (BTU/hr) 


numbe 


thermal  power  ratin 


Reference  primary  loop  fluid  mass  fxowrate 


( lbm/hr ) 


REFT 


Reference  primary  loop  return  temperature  (CF) 


Roof 


urface  convective  film  heat  transfer 


BTU/'h: 


RHO(K,L,N)  Nodal  fluid  density  (lbm/ft-3) 


RMAS 


Desired  fluid  mass  flowrate  in  primary  loop 


( lbm/hr ) 


RPE(L) 


( lbf-hr/f t" ) 


Power  plant  rating  to  meet  total  system  electri 


cal  demand  at 


ed  generation  efficiency 


Power  plant  electrical  generation  capacity  (KWe) 


RSAB(I) 


Roof  material  solar  absorptivity  for  Building 


Power  plant  base  - load  thermal  output  (MWt) 


RTBL 


Building  Type  I ground-contact  material  thermal 
resistance  (hr-ft2-°F/BTU) 


RTHRM 

RTR(I) 

RTW(I) 

RTWD(I) 


SDF 

SHGF 

STLD 

STOR 

SUMM 

SUMT 

SWA 

T(J) 

TAIR(I) 

TAIR ( 25 ) 

TAM3 


Power  plant  base-load  thermal  output  (BTU/hr) 

Roof  thermal  resistance  of  Building  Type  I 
(hr-ft2-°F/BTU) 

Wall  thermal  resistance  of  Building  Type  I 
(hr-ft2-°F/BTU) 

Window  thermal  resistance  of  Building  Type  7 
(hr-ft2-°F/BTU) 

Sky  diffusion  factor  (dimensionless) 

Window  solar  heat  gain  factor  (ETU/hr-ft2) 

Total  thermal  utility  system  load  including 
transmission  losses  (MWt) 

, -5  m 

Reservoir  hot  water  section  volume  (ft'xlO  ) 

Sum  of  load  center  heat  exchanger  fluid  mass 
flowrates  (lbm/hr) 

Sum  of  mass-flow-weighted  load  center  heat  exchan- 
ger outlet  temperatures  (lbm-°F/hr) 

Building  orientation  relative  to  south  (radians) 
Coefficient  in  window  solar  heat  gain  formula 
for  transmission  of  incident  radiation 
Ambient  air  temperature  at  hour  I of  simulation 
(°F) 

Ambient  air  temperature  at  beginning  of  simula- 
tion (°F) 

Ambient  earth  temperature  (°F) 


232. 

TC3  Power  plant  coolant  temperature  at  outlet  of 

precooler  (°F)  -i 

TCW  Reservoir  cold  water  section  temperature  (°F) 

TCWI  Reservoir  cold  water  section  temperature  (°?) 

* 

TD  Nodal  outlet  temperature  difference  (°F) 

TDX  Temperature  difference  across  nodal  boundary 

one  incremental  flow  length  from  outlet  (°F) 

TELD  Total  electrical  demand  from  consumer's  space 

conditioning  and  domestic  hot  water  equipment 
(BTU/hr,  or  MWe) 

TELLD  Total  system  electrical  load  (MWe) 

TEO  Secondary  loop  fluid  reference  temperature  (°F) 

THG  Power  plant  coolant  temperature  at  inlet  to  pre- 

cocler  (°F) 

THI  Reservoir  hot  water  section  temperature  (°F) 

THS  Primary  loop  fluid  reference  temperature  (°F) 

TKW  HTW  temperature  setting  at  inlet  to  primary 

loop  ( °F) 

THWI  Reservoir  hot  water  section  temperature  (°F) 

TT(K,L,N)  Nodal  fluid  inlet  temperature  (°F) 

7IAB  Absolute  temperature  to  be  maintained  in  build- 

ing ( °K ) 

ft 

TIB(NM)  Temperature  to  be  maintained  in  Building  Type 


NN  ( °F ) 


233. 


. < 
r! 
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TIBM 

TIME 

TIP 

TITLE 

TLD 

TLH(IX) 

TL? (IX) 

TM 

TMAS 

TN 

TO ( K , L , N ) 
TOAB 
TOB 
TODP 


TOP ( K , L ,N) 


TPO 

TQ(IX) 


f 

f 

I 

» 

f 

Minimum  of  temperatures  to  be  maintained  in  all 
building  types  (°F) 

Elapsed  time  from  beginning  of  simulation  (hr) 

Load  center  heat  exchanger  fluid  inlet  tempera- 
ture at  end  of  time  step  (°F) 

Alphanumeric  run  title 

Total  community  space  conditioning  and  domestic 
hot  water  energy  demand  (BTU/hr,  or  MWt ) 

Lower  temperature  limit  setting  for  supply  water 
to  consumers  served  by  heat  exchanger  IX  (°F) 

Lower  temperature  limit  setting  for  water  at 
outlet  of  primary/secondary  heat  exchanger  IX  (°F) 

Heat  exchanger  temperature  difference  ratio 
Total  mass  of  fluid  contained  in  a node  (lbm) 

Nodal  inlet  temperature  difference  (°F) 

Nodal  fluid  outlet  temperature  (°F) 

Absolute  ambient  air  temperature  (°K) 

Ambient  air  temperature  (°F) 

Fluid  temperature  at  one  incremental  flow 
length  from  nodal  outlet  (°F) 

Nodal  fluid  outlet  temperature  at  end  of  time 
step  ( °F ) 

Total  pressure  loss  around  loop  (lbf/ft  ) 

Thermal  demand  at  load  center  heat  exchanger 
IX  (BTU/hr) 

J 


' ■ " 


TQL(IHR,IX)  .Thermal  demand  at  load  center  heat  exchanger 


IR(L) 


?RP ( IX ) 


TSAW 


TSB(I,NN) 


TSB( 25  j I ) 


TSET 


TTLD 


TUH(IX) 


TUL(IX) 


IX  at  simulation  hour  IHR  (ETU/hr) 
Hass-flowrate-averaged  fluid  temperature  at 
load  center  heat  exchanger  return  header  for 
loop  L (°P) 

Design  temperature  setting  of  water  at  outlet 
of  primary/secondary  neat  exchanger  IX  (°F) 
Sol-air  temperature  at  roof  surface  (°F) 
Sol-air  temperature  at  exterior  wall  surface 


Temperature  to  be  maintained  in  building  type 
NN  at  hour  I of  simulation  (°F) 

Temperature  to  be  maintained  in  Building  Type 
I at  beginning  of  simulation  (°F) 

Heat  exchanger  steady-state  reference  tempera- 
ture (°F) 

Total  thermal  demand  at  a load  center  heat 
exchanger  (BTU/hr) 

Upper  temperature  limit  setting  for  supply 
water  to  consumers  served  by  heat  exchanger 
IX  (°F) 

Upper  temperature  limit  setting  for  return 
water  from  consumers  served  by  heat  exchanger 
IX  (°F) 


.1.. - 


Upper  temperature  limit  setting  for  water  at 
outlet  of  primary/secondary  heat  exchanger  IX 


VLEAK 


Total  air  exchange  between  interior  and  exterior 


of  building  due  to  ventilation  (CFH) 


Orientation  of  Building  Type  I with  respect  to 


North  (0°<WA<90°) 


WDCL(I) 


Total  length  of  cracks  around  all  windows  in 


Building  Type  I (ft) 


WDINC(I) 


Infiltration  air  flow  coefficient  C for  windows 


in  Building  Type  I (CFM/ft-in.  H~0) 


WDINN(I) 


Infiltration  air  flow  coefficient  N for  windows 


in  Building  Type  I (dimensionless) 


WDIR(I) 


Wind  direction  at  hour  I of  simulation  (degrees 
clockv/ise  from  north) 


WDIR ( 25 ) 


Wind  direction  at  beginning  of  simulation 
(degrees  clockwise  from  north) 


Total  air  infiltration  through  cracks  around 


windows  (CFM) 


WDSC(I) 


Window  interior  shading  coefficient  for  Building 
Type  I (decimal  fraction;  0=no  solar  transmission 


l=no  shading) 


Exterior  wall  convective  film  heat  transfe 


coefficient  (BTU/hr-ft  -°F) 


Surface  roughness  coefficients  for  Building 


Inf iltr 


Infiltration  air 


low  coefficient  N for  walls 


In  Building  Type  I (dimensionless) 


Total  air  infiltration  through  cracks  in  walls 


WSAB(I) 


Wall  material  solar  absorptivity  for  Building 


Aggregate  window  solar  heat  gain  factor 


( BTU/hr-ft^ ) 


Use  factor  for  domestic  hot  water  in  building 
type  NN  at  hour  I of  simulation  (decimal  fraction) 


ak  demand  rate 


domestic  hot  water  in  Building 


imulation  (decimal 


1 = use  at  peak  demand  rate 


?37. 


WV(I) 

wv(?5) 


XHAT 


Wind  velocity  at  hour  I of  simulation  (mph) 
Wind  velocity  at  beginning  of  simulation  (mph) 
Heat  exchanger  energy  supply  scale  factor 
( hr-°F/BTU ) 

Relative  direction  of  wind  relative  to  normal 
to  wall  I (radians) 

Ratio  of  fluid  mass  flowrate  through  primary 
heat  exchanger  to  primary  loop  mass  flowrate 


Diffuse  solar  radiation  incidence  coefficient 


REAL  LONG, MASS, MtQMGAS  MMN0001 
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